%ex0808 用ode23 ode45 ode113解多階微分方程
clear,clc
[x23,y23]=ode23('myfun03',[1,10],[1 10 30]);
[x45,y45]=ode45('myfun03',[1,10],[1 10 30]);
[x113,y113]=ode113('myfun03',[1,10],[1 10 30]);
figure(1) %第一幅圖
plot(x23,y23(:,1),'*r',x45,y45(:,1),'ob',x113,y113(:,1),'+g') %作出各種函數(shù)所得結(jié)果
legend('ode23解','ode45解','ode113解')
title('ODE函數(shù)求解結(jié)果')
figure(2)
plot(x45,y45)  %以ode45為例作出函數(shù)以及其各階導(dǎo)數(shù)圖
legend('y','y一階導(dǎo)數(shù)','y兩階導(dǎo)數(shù)')
title('y,y一階導(dǎo)數(shù),y二階導(dǎo)數(shù)函數(shù)圖')