I have to plot these equations about the heat in a aluminium bar:
But for some reason, my graph of u(x,t)(in this case u on the code) is showing something different (i think the correct would be a exponential). Can someone review my code and see if is there any mistakes? Obs: n is the value where sin=0, for example sin(n*pi)=0 like in pi, 2pi, 3pi...
clc
clear all
syms x k t n
for n=1:1:20
cn=(100/(n*pi))*(cos(pi*x/4)-cos(3*pi*x/4));
for t=0:5
num=(cos(pi*x/4)-cos(3*pi*x/4));
t1=exp((-(k^2)*(pi^2)*t)/1600);
t2=sin(k*pi*x/40);
S1=symsum((num*t1*t2),k,1,inf);
u=((100/pi)*S1);
end
end
fplot(cn)
figure
fplot(u)
