常微分方程初值问题_Matlab代码.pdfVIP

常微分方程的初值问题 1.DEEuler 用欧拉法求一阶常微分方程的数值解 function y = DEEuler(f, h,a,b,y0,varvec) format long; N = (b-a)/h; y = zeros(N+1,1); x = a:h:b; y(1) = y0; for i=2:N+1 y(i) = y(i-1)+h*Funval(f,varvec,[x(i-1), y(i-1)]); end format short; 2.DEimpEuler 用隐式欧拉法求一阶常微分方程的数值解 function y = DEimpEuler(f, h,a,b,y0,varvec) format long; N = (b-a)/h; y = zeros(N+1,1); y(1) = y0; x = a:h:b; var = findsym(f); for i=2:N+1 fx = Funval(f,var(1),x(i)); gx = y(i-1)+h*fx - varvec(2); y(i) = NewtonRoot(gx,-10,10,eps); end format short; 3.DEMod