char4_6.m

matlab 文件

t=[0,50];
y0=[0,0.25];
[t,y]=ode45('xp',t,y0);
plot(t,y(:,1),':b',t,y(:,2),'-r')
legend('速度','位移')
pause
clf;
axis([0,50,-4,4]);
hold on
ode23s('xp',t,y0)
legend('速度','位移')