function?dy=myfun03(x,y)

dy=zeros(3,1)?%初始化變量dy

dy(1)=y(2);?%dy(1)表示y的一階導(dǎo)數(shù),其等于y的第二列值

dy(2)=5/x*y(3)-y(1);?%dy(2)表示y的二階導(dǎo)數(shù)

%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ù)圖')

總結(jié)

以上是生活随笔為你收集整理的matlab求微分数值,用MATLAB语言求微积分方程的数值解.(xd^2y)/dx^2-5dy/dx+y=0y(0)=0y'(0)=0...的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。