What did i do wrong here?
显示 更早的评论
L=1;
alpha= 10^-3;
T1=800;
T2=200;
x=0:0.01:L;
N=50;
t=0:5:25;
A_0=(T1+T2)/2;
Teaxct=zeros(length(x),length(t));
for m=1:length(x)
for k=1:5:length(t)
for n=1:(N)
A_n(n)=(2*(T1-T2)/pi*(n))*(sin((pi/2)*(n)))*(cos((pi*n*m)/L))*(exp(-1*alpha*((pi*(n))/L)^2*k));
Texact(m,k)=A_0+A_n(n)
end
end
end
2 个评论
KSSV
2020-8-25
Why three loops? Type error here in the variable Teaxct.
You can achieve what you want with meshgrid.
Muhammad Umar Khan
2020-8-25
回答（1 个）
Cris LaPierre
2020-8-26
0 个投票
In your third loop, the values of m and k do not change. This means the results of each loop are written to the same element of Texact, overwriting the previous value. The result is you only capture the final value in your third loop.
Your code also runs really slow because you are printing out the results of every single interation.
You also have a typo I suspect. You initialize Teaxct but then use Texact in your loops.
2 个评论
Muhammad Umar Khan
2020-8-26
Cris LaPierre
2020-8-26
Fix your typo, add a semicolon to the end of your Texact line, and then describe what your expected output size is. If you want to capture every value, you need to add a third dimension to Texact
Texact(m,k,n)= ...
类别
在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息
2020-8-25
2020-8-26
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
