【求助】【急】求高手知道Matlab编程，利用函数分析法画阻尼振动曲线

哎，老了，matlab生疏了。
哥们啊，我一个心血来潮给你认真了一个下午才弄出来的，要求追加分数！
程序代码如下：
beta=0.1;
w0=1;
delta_t=25/1000;
t(1)=0;
for i=1:1000
t(i+1)=delta_t*i;
end

threshold=0.9;
x(1)=100; %为了求解方便，最后再更改比例尺
temp_v1=0;
n=1; %x的角标
for i=1:1000
for j=-100:0.001:100
temp_x=j;
temp_v2=(temp_x-x(n))/delta_t;
temp_a=(temp_v2-temp_v1)/delta_t;
if abs(temp_a+2*beta*temp_v2+w0^2*temp_x) n=n+1;
x(n)=temp_x;
temp_v1=temp_v2;
break;
end
end
end
for i=1:1001 %回归真正的比例
x(i)=x(i)/100;
end
plot(t,x);
xlabel('Time');
ylabel('X-distance');
title('弹簧振子阻尼振动曲线');

直接运行便可以得到弹簧振子阻尼振动曲线。