I have this code and want to eliminate the inner most loop and use sum
How can I do that?
thetao=1;alfa=1;l=1;theta2(1:7,1:21)=0;T=[0.01, 0.1, 0.2, 0.5, 1, 2, 5];x=0:0.05:l;n=0:50;for t=1:7 for j=1:length(x) for k=1:length(n) theta1(t,j)=(sin((2*n(k)+1)*pi*x(j)/(2*l))*exp(-(2*n(k)+1).^2*pi^2*alfa*T(t)/(4*l^2))/((2*n(k)+1)*pi)); theta2(t,j)=theta2(t,j)+theta1(t,j); theta(t,j)=thetao-4*thetao*theta2(t,j); end endend
Best Answer