tt = linspace(0,12,41);
x = [0:0.01:6];

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

for i =1:nn
f(:,i) = fwave(N,tt(i));
plot(x,f(:,i))
axis([0,6,-1.5,1.5]);
text(5,1.1,['t = ',num2str(tt(i))])
M(:,i) = getframe;
end