figure(1)title('Exact Steady State Temperature Profile')hold onplot(r_profile,Up_exact,'r-');xlabel('r')ylabel('T')for i = 1:ceil(tau/dt) Mplot(:,i) = U; U = Atinv*(D*U+dt*B); [Uedge, redge] = GetBHboundary(theMesh,3,U); figure(2) plot(redge, Uedge,'k') axis([1 10 0 0.4]) title('FE Temperature Profile') drawnow pause(0.05)end
I have the following code which plots my steady state temperature profile in figure 1 and an animated transient temperature profile in figure 2. Is there a way i can have them both in one plot showing how the transient solution approaches the steady state?
Thanks 🙂
Best Answer