% script for animation of wave traveling
for ii=1:100,
    plot(Xfin,Y2_save*(am_cnv'.*cos(list*0.5*ii)));
    axis([0 1.01 -1.5e2 1.5e2]);
    M(ii) = getframe;
end
