实验四常微分方初程值问题数值解法.docVIP

实验四 常微分方程初值问题数值解法 1、改进Euler公式 function [x,y]=euler_r(ydot_fun, x0, y0, h, N) % 改进Euler公式，其中 % ydot_fun --- 一阶微分方程的函数 % x0, y0 --- 初始条件 % h --- 区间步长 % N --- 区间的个数 % x --- Xn 构成的向量 % y --- Yn 构成的向量 x=zeros(1,N+1); y=zeros(1,N+1); x(1)=x0; y(1)=y0; for n=1:N x(n+1)=x(n)+h; ybar=y(n)+h*feval(ydot_fun, x(n), y(n)); y(n+1)=y(n)+h/2*(feval(ydot_fun, x(n), y(n))+feval(ydot_fun, x(n+1), ybar)); end 书本第九章习题1解答 ydot_fun=inline(1./x^2-y./x,x,y); [x,y]=Euler_r(ydot_fun,1,2,0.05,20) x = Columns 1 through 13 1.0000 1.0500 1.