最简单求解一个微分方程数值解得方法:Euler法 function [x,y]=Euler_method(dufun,span,h,x0,y0) %EuLer格式, %求解方程y'=dufun(x,y);其中x \in[a,b];y0为初始值:n为自变量的离散个数:y为求解结果 x=span(1):h:span(2); n=length(x); y=zeros(1,n);%存放数值的解 x(1)=x0; y(1)=y0; for i=1:n-1 y(i+1)=y(i)+h.*feval(dufu…