数值分析-第九章计算实习题答案-昆工.docxVIP

⑻程序: cic;clear; a二l;b=2;%定义域h = 0.05;% 步长 n=(b-a)/h;y0=l;%初值 f= @(x,y) 1/x 八 2-y/x;% 微分函数Xn=linspace(a,b,n+1);%将定义域分为n等份 Yn 二 zeros(Ln);% 结果矩阵Yn⑴=y0;%赋初值 %以下根据改进欧拉公式求解for i=l:n xn=Xn(i); xnn=Xn(i + l); yn=Yn(i); yp=yn+h*f(xn,yn); yc=yn + h*f(xnn,yp); yn=(yp+yc)/2; Yn(i + l)=yn;end Xn=Yn;%以下根据经典四阶R-K法公式求解 for i=l:n xn = Xn(i); yn=Yn(i); kl=f(xn,yn); k2=f(xn+h/2,yn + h/2*kl); k3=f(xn+h/2,yn + h/2*k2); k4=f(xn+h,yn+h*k3); yn = yn+h/6*(kl+2*k2+2*k3+k4); Yn(i+l)=yn;end dispC改进欧拉法 四阶经典R-K法)disp([Xn Yn]) 结果如下：改进欧拉法10.99887 结果如下： 改进欧拉法 1 0.99887 0.99577 0.99114 0.98532 0.97857 四阶经典R-K法 1 0.99885 0.997