数值计算第3次大作业答案宋宝瑞.docxVIP

PAGE \* Arabic 4 / NUMPAGES \* Arabic 4 数值计算第三次大作业 Euler法 function y=euler(x) %实现euler法 %输入值x代表x=x0时，y0的值 x1=x(1); y1=x(2); h=1E-2; n=1; u(1)=y1; for x=x1+h:h:2*pi n=n+1; y2=y1+h*fx(x); u(n)=y2; y1=y2; end y=y2; %绘图 x=x1:h:2*pi; plot(x,u); grid on xlabel(角度/ °); ylabel(函数值); title(Euler法); 结果为：y=5.0000 梯形法 function u=tixingfa(x) %实现euler法 %输入值x代表x=x0时，y0的值 x1=x(1); y1=x(2); h=1E-2; n=1; y(1)=y1; for x=x1+h:h:2*pi y(n+1)=y(n)+h/2*(fx(x)+fx(x+h)); n=n+1; end u=y(n); x=x1:h:2*pi; plot(x,y); grid on xlabel(角度/ °); ylabel(函数值); title(梯形法); 结果为：y=4.9999 龙格库塔法 functio