tt = linspace(0,2*pi,31);
x = [0:0.005:pi];

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

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