ex3-13-3.m

matlab程序

%  建立函数文件xprime.m

%  解微分方程

t0=0;  tf=20;

x0=[0,0.25];

[t,x]=ode23('xprime',t0,tf,x0);

plot(t,x(:,1),':b',t,x(:,2),'-r')

legend('速度','位移')