显式Euler、梯形方法、预校Euler方法的比较.docVIP

微分方程数值解 姓名 丁建伟 学号 200708020211 日期 2010．10．6 实验项目 显式Euler、梯形方法、预校Euler方法的比较考虑一阶常微分方程初值问题dy/dx = -y^2, 0x=1, y(0)=1，其精确解为y = 1/(1+x),使用三种方法求初值问题数值解。给出步长h=1/16,1/32,1/64,1/128,1/256时的x=1点处的整体截断误差E(h)=|yn - y(1)|，并从误差、收敛阶、计算量方面比较三种方法。 目的与要求：掌握三种方法的程序实现掌握比较算法优缺点的方法原理：显式Euler方法： 梯形方法： 预校Euler方法： 算法描述：由已知条件和根的存在区间，经过变换，三种方法的迭代公式分别为： 显式Euler方法： 梯形方法： 预校Euler方法： 主要程序代码或命令： #includestdio.h #includemath.h //需要用到绝对值abs（）和求平方根sqrt（）函数 void main() { int n; int i; printf(请输入n的值：); scanf(%d,n); float y=1.0,y0=1.0,y1=1.0; //赋初始值。 float e,e0,e1; //绝度误差 float h; h=1/float(n);