数值逼近实习报告八(龙贝格).docVIP

数 值 逼 近 实 习 报 告 班级 姓名 指导教师 实习题目 龙贝格积分算法 计算公式： 目的意义 掌握龙贝格积分公式的计算方法 数学模型 如果将区间等分，则复化梯形积分公式为 现变步长，将区间等分，则有 令 算法描述 对于给定的和初始 Step1：用复化梯形公式 求得外推出 Step2: 对于已经计算出来的 计算出 Step3: 用外推法一次计算 Step4： 输出 （流程图）程序 #include iostream.h #include math.h const int k=4; //次数 const double a=0; //起点 const double b=1; //终点 double h=b-a; //步长 double (*data)[k]=new double[k][k]; double f(const double x); //求积函数 double sum(const int n,const double h); //多项式求和 int main() { cout ///////////////////////////////////////endl; cout // Ronberg积分公式 //endl; cout ///////////////////////////////////////endl; cout利用Ronberg积分公式计算圆周率πendl; cout步长为(b-a)的时候积分结果为：endl; data[0][0]=h/2*(f(a)+f(b)); coutT0=; coutdata[0][0]endl; cout变步长之后的计算结果为：endl; for (int i=1;ik;++i) { data[i][0]=data[i-1][0]/2+h/2*sum(i,h); //复化梯形积分计算公式 h=h/2; } for(i=1;ik;i++) for (int j=i;jk;j++) data[j][i]=(pow(4,i)*data[j][i-1]-data[j-1][i-1])/(pow(4,i)-1); for (i=0;ik;i++) { for (int j=0;j=i;j++) coutTj+1=data[i][j] ; coutendl; } } double f(const double x) { return 4/(1+x*x); } double sum(const int n,const double h) { int t=pow(2,n-1); double s=0; for (int i=1;i=t;i++) { s=s+f(a+(2*i-1)*(h/2)); } return s; } 数值算例 对计算结果进行分析评价 参考文献 [1] ***，数值逼近，***：****大学***学院，20**.