2010数值积分方法.docVIP

实验二 数值积分法（4课时） 实验目的： 掌握复化辛卜生求积公式，理解Romberg方法及其求法，能写出其基本算法及程序。 实验要求： 按照题目要求完成实验内容 写出相应的C语言程序． 给出实验结果． 对实验结果进行分析讨论． 写出相应的实验报告． 实验内容： 分别用复化辛卜生求积公式及龙贝格方法上机计算以下积分结果（取eps=5e-6) 程序框图 龙贝格积分算法 程序框图 龙贝格实验程序算法： float Romberg(float a,float b,float(*f)(float),float eps) { int n=1,k; float h=b-a,x,temp; float T1,T2,S1,S2,C1,C2,R1,R2; T1=(b-a)/2*((*f)(a)+(*f)(b)); while(1) { temp=0; for(k=0;k=n-1;k++) { x=a+k*h+h/2; temp+=(*f)(x); } T2=(T1+temp*h)/2; if(fabs(T2-T1)eps) return T2; S2=T2+(T2-T1)/3; if(n==1) {T1=T2;S1=S2;h/=2;n*=2;continue;} C2=S2+(S2-S1)/15; if(n==2) {C1=C2;T1=T2;S1=S2;h/=2;n*=2;continue;} R2=C2+(C2-C1)/63; if(n==4) {R1=R2;C1=C2;T1=T2;S1=S2;h/=2;n*=2;continue;} if(fabs(R2-R1)eps) return R2;R1=R2;C1=C2;T1=T2;S1=S2;h/=2;n*=2; } } 用补偿法求新龙贝格积分R2 用补偿法求新柯特斯积分C2 用补偿法求新辛卜生积分S2 求新的梯形面积T2 新增加结点函数值 判断积分近似值是否达到精确度要求 为求新积分准备二分区间 输出 求积分值 数据类型说明定义 函数输入 使用龙贝格法求函数的积分 用补偿法求新龙贝格积分R2 用补偿法求新柯特斯积分C2 用补偿法求新辛卜生积分S2 求新的梯形面积T2 新增加结点函数值 判断积分近似值是否达到精确度要求 为求新积分准备二分区间 输出 求积分值 数据类型说明定义 函数输入 使用龙贝格法求函数的积分 用辛卜生积分公式求S值并返回 用辛卜生积分公式分解算法求S1，S1 新增加结点函数值 为求新积分准备n分区间， 输出 求积分值 数据类型说明定义 被积函数的定义 使用复合辛卜生法求函数的积分