实验四数值积分.docVIP

数学与计算科学学院 实 验 报 告 实验项目名称 数值积分 所属课程名称 数值计算 实 验 类 型 验证 实 验 日 期 2015年5月17日 班 级 桥梁1301班 学 号 201388250115 姓 名 殷兆镶 成 绩 实验四 数值积分 一、实验概述： 【实验目的】 掌握复合梯形公式和复合抛物线公式求积分的方法 【实验原理】 复合梯形公式和复合抛物线公式求积分的基本原理和公式 【实验环境】 二、实验内容： 【实验方案】 教材p280数值试验七 1.复合求积公式计算定积分 实验题目(1)(2)： 要求：① 将区间10等分，分别用复合梯形公式和复合抛物线公式求定积分； ② 将区间n等分，比较当n = 10、100①”的实验过程 【实验结论（结果）小结/*#includestdio.h void main() {printf(hello); }*/ #includestdio.h #includemath.h double f(double x) {return -2 ; } void main() { double h,T=0.0,s; int k,n=100;double a=2,b=3; h=(b-a)/n; for(k=1;kn;k++) {T+=f(a+k*h);} s=(h/2.0)*(f(a)+2*T+f(b)); printf(结果是:%lf\n,s); } 实验题目(2) /*#includestdio.h void main() {printf(hello); }*/ #includestdio.h #includemath.h double f(double x) {return ; } void main() { double h,T=0.0,s; int k,n=100;double a=1,b=2; h=(b-a)/n; for(k=1;kn;k++) {T+=f(a+k*h);} s=(h/2.0)*(f(a)+2*T+f(b)); printf(结果是:%lf\n,s); 复合抛物线公式: #includestdio.h #define Len 100 void main() { int i; double n,a,b; double y[Len]; double y1=0; double y2=0; double h,S,fa,fb; printf(请输入端点a,b的值和n:\n); scanf(%lf %lf %lf,a,b,n); h=(b-a)/(2*n); printf(%lf\n,h); printf(请输入端点的函数值fa,fb:\n); scanf(%lf %lf,fa,fb) for(i=1;i2*n;i++) { printf(请输入函数值y[%d]:\n,i); scanf(%lf,y[i]); } for(i=1;i=n-1;i++) { y1=y1+y[2*i]; } for(i=1;i=n;i++) { y2=y2+y[2*i-1]; } S=h/3*(fa+4*y2+2*y1+fb); printf(S=%lf\n,S); 1