显示欧拉方法与梯形方法.docVIP

《数值分析》实验报告 姓名 学号 日期 2012．11.20 实验室 设备编号 实验题目 用显示欧拉方法和梯形方法求y/=-y+x+1 一 实验目的 1.掌握显示欧拉格式和梯形格式。 2．用显示欧拉方法和梯形方法求y/=-y+x+1,y(0)=1,0≤x≤0.5 二 实验内容 算法：设在点xn列出方程y’(xn)=f(xn,y(xn))并用商差y(xn+1)- y(xn)/h替代其中的导数项y’(xn),则有y(xn+1)≈yn+hf(xn,y(xn))若用y(xn)的近似值yn代入上式右端，并记所得结果为yn+1得欧拉公式: yn+1= yn+ hf(xn,yn)，n=0,1,2,… 设方程y’=f(x,y)的两端从xn到xn+1求积分，即得梯形求积公式y(xn+1)=y(xn)+∫f(x,y(x))dx，梯形格式为yn+1= yn+h/2[f(xn,yn)+ f(xn+1,yn+1)]。 源程序： 1)显示欧拉方法 #includeiostream #includemath.h #includeiomanip using namespace std; class euler { private: int i,n;//i步计数器，n等分数 float f,h,x,y,yb;//f表示F(x,y),h步长，x自变量 public: void