用复合梯形求积公式求积实验报告.docxVIP

用复合梯形求积公式求积 实验报告 学院：计算机学院 班级：信计 0902 班 姓名：高天宇 学号： 0908060220 1、 目的 1） 察复化梯形公式和复化辛普森公式随区 数  n 增加各自 差的减少 律；研究广 分的数 算如何将其 化 普通 分，再由已有数 分方法 行 算； 2）利用复化梯形公式和复化辛普森公式 算定 分， 程 。 2、 内容 分 用 复 化 梯 形 公 式 和 复 化 辛 普 森 公 式 算 定 分 ， 研 究 随 着 n 增 加 各 自 差 的 减 少 律 2 1 ex dx ，取 n=2,4,8,16 ，精确 I=4.006994 。 I 0 3、 原理 原理：将区 [a,b] 等分成 N 个子区 [x(k),x(k+1)] （k=0,1, ?..,N -1 ）， h=（ b-a ） /N ，在每个子区 [x(k),x(k+1)] 上使用体形公式，可求得 果。 4、 思想 1 x ln xdx ，然后再将 分区 用上述的的公式来 算 I [0 ，1] 行 n 等分，得到 n 个小区 [x i ,x i+1 ] 0 （ i=0 ， 1 ， 2 ， ? ， n-1 ） , 每个小区的度 h=(b-a)/n, 其 中 xi =a+i*h 。 按 b h n 1 h [ f (a) n 1 T f (x)dx [ f ( xi ) f (xi 1 )] f (b)] h f ( xi ) 行 算。 a 2 i 0 2 i 0 5、程序框 6、 程序 #includeiostream.h #includemath.h double fun1(double a,double h,int n) { int i; double sum=0,m; for(i=1;i=n-1;i++) { m=a+i*h; if(m==0)continue; sum=sum+sqrt(m)*log(m); } return sum; } double get(double a) { if(a==0)a=0.00001; return sqrt(a)*log(a); } void main() { double a=0,b=1; int n; while(1) { cout 请输入要等分的子区间个数 (n0):; cinn; double h; double T,f0,f1,f2; h=(b-a)/n; f0=get(a); f1=fun1(a,h,n); f2=get(b); T=h/2*(f0+2*f1+f2); cout 步长 h 为 :hendl; cout 复化梯形公式计算结果 T=Tendl; } } 7、实验结果 8、实验体会 通过本次实验我熟悉了用复化梯形公式求数值积分的全过程，并更加熟悉的掌握了复化梯形公式原理。 使我对计算方法这门课有了进一步的理解，为以后的学习打下了良好的基础。