ex3-13-3.m

来自「matlab程序」· M 代码 · 共 15 行

%  建立函数文件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('速度','位移')