tt = [0:0.01:.1,.3:.2:5];
x = [0:0.005:pi];

nn = length(tt)
clf
M = moviein(nn);

for i =1:nn
f(:,i) = fheat(N,tt(i));
plot(x,f(:,i))
axis([0,pi,0,12]);
text(2.75,11,['t = ',num2str(tt(i))])
M(:,i) = getframe;
end