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