I got all the curves to plot in the top subplot by adding matrices that stored the various values of the vectors you plot in the loops. (Before, it was only plotting the last vector of ‘T1’ and ‘T2’ calculated.) The curve in the top subplot looks a bit strange, but I’ll let you sort that, since I don’t know what you’re doing or how it should look.
I’m only listing the loops here, since those are the only changes I made:
for i=1:nxp
depth(i)=i*dx;
for j=2:ntp
time(j)= j*dt;
Tx=Tsd+((Ti-Tsd)*(erf(depth(i)/(2*sqrt(aLf*(time(j)))))));
T1(j)=Tx;
end
subplot(223);
plot(t,T1); hold on; xlabel('time (t)'); ylabel('Temperature(Tx)');grid on
Tni(i)=Tx;
T1M(i,:) = T1; % Added: Store ‘T1’ Vectors In ‘T1M’
end
time2=time;
time2(1)=1;
for k=1:nxp
depth(k)=k*dx;
for r=1:ntp
Ty=Tsn+((Tni(k)-Tsn)*(erf(depth(k)/(2*sqrt(aLf*(time2(r)))))));
T2(r)=Ty;
end
Tnii(k)=Ty;
subplot(224);
plot(t,T2); hold on; xlabel('time (t)'); ylabel('Temperature(Tx)');grid on
T2M(i,:) = T2; % Added: Store ‘T2’ Vectors In ‘T2M’
end
T=[T1M T2M]; % Concatenate ‘T1M’ And ‘T2M’ To Form ‘T’
subplot(211);
plot(y,T); hold on; xlabel('time (t)'); ylabel('Temperature(Tx)');grid on
