数值分析上机(九).docxVIP

PS：均来自数值分析第五版 李庆扬，王能超，易大义编出版社：清华 大学出版社微分方程问题的初值解法 欧拉及四阶龙格库塔 问题局部： 给定初值问题： V-—7, 1 V x ( 2, XXy⑴=1 y(2) =- 50P + 50/ + 2x, 0xly(o) = 1 y (2) 要求：(a)用改进欧拉法(h=0.05)及经典四阶R-K法(h=0.1)求⑴的数值解，并打印x=l+0.1i 的值 (b)用经典四阶R-K方法解(2),步长分别取h=0.l,0.025,0.0L计算并打印x=0.1i(i=0...10) 各点的值，与准确解y(x)=l/3exp(-50*(x))+x八2;代码局部: /**R_K.m**/function z = R_K(h,a,b,yO) %a,b输入参数：边界输入参数:步长 %y0——_输入参数：初值%y———输出:函数值 x=a:h:b; N=(b-a)/h+l;y(i)=yO； for i=l:N-lKl=fun(x(i),y(i)); K2二fun(x ⑴+h/2,y ⑴+h/2*K 1);K3=fun(x ⑴+h/2,y ⑴+h/2*K2); K4=fun(x(i)+h,y(i)+h*K3);y(i+l)=y ⑴+h/6*(Kl+2*K2+2*K3+K4); end z=y;end 函数局部 /***fun.m**/ function z