一般地，凡表示未知函数，未知函数的导数与自变量之间的关系的方程叫做微分方程未知函数是一元函数的，叫常微分方程;未知函数是多元函数的，叫做偏微分方程,常微分方程的数值解,如,数值解的龙格库塔方法,龙格库塔方法的一般形式,其中,t,x=ode23(f,ts,x0,options,p1,p2,) t,x=ode45(f,ts,x0,options,p1,p2,)function dx=f(t,x) dx=f1;f2;fn,Matlab实现,options=odeset(reltol, rt, abstol , at) 默认时rt=10-3,at=10-6,Matlab实现,食饵捕食者模型,表示时刻 食饵的密度， 表示捕食者的密度；表示食饵独立生存时的增长率；表示捕食者独立生存时的死亡率；表示捕食者的存在对食饵增长的影响系数，反映捕食者对食饵的捕获能力；表示食饵的存在对捕食者增长的促进系数，反映食饵对捕食者的喂养能力,高阶常微分方程的解法,高阶常微分方程,可以将以上高阶微分方程化为如下一阶常微分方程组,令,Bessel方程,令 时求解,解：首先将其化为方程组,令,首先编译如下函数M-文件,然后编译如下M-文件，用求解常微分方程的命令ode45调用以上函数文件，画图比较所求结果和精确解,运行结果,
