%使用電影方式制作動畫
n=20;
for i=1:n
 x=0:0.1:i;
 y=1-1/sqrt(1-0.3^2)*exp(-0.3*x).*sin(sqrt(1-0.3^2)*x+acos(0.3));
 plot(x,y)
 axis([0,20,0,1.5]);  %固定坐標軸
 M(i)=getframe; %抓取畫面
end
movie(M,3)  %播放3次