figure('color',[1 1 1])
plot([-0.2;0.2],[0;0],'color','y','linestyle','-','linewidth',10);%畫(huà)懸線的橫梁
g=9.8;                                                            %重力加速度
l=1;                                                              %線長(zhǎng)
theta0=pi/6;                                                      %初始角度
x0=l*sin(theta0);                                                 %初始 x 值
y0=-l*cos(theta0);                                                %初始 y 值
axis([-0.75,0.75,-1.25,0]);                                       %畫(huà)坐標(biāo)范圍
axis('off');                                                      %關(guān)閉坐標(biāo)顯示
head=line(x0,y0,'color','r','linestyle','.','erasemode','xor','markersize',40);    %定義小球格式
%mmonth=line(x0,y0,'color','g','linestyle','.','erasemode','xor','markersize',30);
body=line([0;x0],[0;y0],'color','b','linestyle','-','erasemode','xor');            %定義懸線格式
t=0;                                                              %時(shí)間初值
dt=0.005;                                                          %時(shí)間增量
while 1                                                           %死循環(huán)(關(guān)閉窗口后中止循環(huán))
   t=t+dt;
   theta=theta0*cos(sqrt(g/l)*t);                                 %計(jì)算時(shí)刻 t 時(shí)的幅角                               
   x=l*sin(theta);                                                %計(jì)算時(shí)刻 t 時(shí)的 x 坐標(biāo)
   y=-l*cos(theta);                                               %計(jì)算時(shí)刻 t 時(shí)的 y 坐標(biāo)
   set(head,'xdata',x,'ydata',y);                                 %計(jì)算時(shí)刻 t 時(shí)小球的位置
   set(body,'xdata',[0;x],'ydata',[0;y]);                         %計(jì)算時(shí)刻 t 時(shí)懸線的位置
   drawnow;                                                       %重畫(huà)小球與懸線
end