n=1;?%?阻尼系數
w=6;?%?圓頻率
ww=sqrt（w^2-n^2）;?%?等效的圓頻率
B1=80;?%?系數
B2=-20;?%?系數
axis（[-105105-105105]）;hold?on;?%?設置坐標軸范圍
set（gcf‘doublebuffer‘‘on‘）;?%?設置圖形窗口的渲染效果
G=plot（1000‘r‘）;?%?軌跡曲線
x=100;?%?初始位置的橫坐標
y=0;?%?初始位置的縱坐標
r=2;?%?小球半徑
t1=r*sin（0:.1:2*pi）;?%?計算小球的邊界數據
t2=r*cos（0:.1:2*pi）;?%?計算小球的邊界數據
X=x+r*t1;?%?得到小球的橫坐標數據
Y=y+r*t2;?%?得到小球的縱坐標數據
Q=fill（XYrand（13））;?%?畫出小球
t=0;?%?時間的初始數值
plot（[-100100][00]）;?%?畫出格線
plot（[00][-100100]）;?%?畫出格線
k=1;?%?控制循環(huán)的參數
xlabel（‘Please?press?“space“?key?and?stop?this?program!‘...
????‘fontsize‘14‘color‘‘r‘）;?%?添加提示字符
title（[‘《damp?turbination》?‘‘B1=