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 comentarios

KSSV
KSSV el 25 de Ag. de 2020
Why three loops? Type error here in the variable Teaxct.
You can achieve what you want with meshgrid.
Muhammad Umar Khan
Muhammad Umar Khan el 25 de Ag. de 2020
One for each variable

Iniciar sesión para comentar.

Respuestas (1)

Cris LaPierre
Cris LaPierre el 26 de Ag. de 2020

0 votos

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 comentarios

Muhammad Umar Khan
Muhammad Umar Khan el 26 de Ag. de 2020
How to resolve it any idea
Cris LaPierre
Cris LaPierre el 26 de Ag. de 2020
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)= ...

Iniciar sesión para comentar.

Categorías

Más información sobre Loops and Conditional Statements en Centro de ayuda y File Exchange.

Etiquetas

Preguntada:

el 25 de Ag. de 2020

Comentada:

el 26 de Ag. de 2020

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by