龙贝格+实验报告.docx

学院：数学科学学院　　　专业：实验名称成绩姓名年级学号实验三 龙贝格求积分实验目的 用算法实现龙贝积分运算，求解I=∫01X3/2dx的积分实验原理 算法设计与程序流程图算法设计：(先假定所求积分二分最大次数次数为20)先求T[k][0] 再由公式求T[i][j]在求出的同时比较T[k][k]与T[k-1][k-1]的大小，如果二者之差的绝对值小于1e-5，就停止求T[k][k]；此时的k就是所求的二分次数，而此时的T[k][k]就是最终的结果打印出所有的T[i][j]； 程序流程图： m=n N Yk=N N m=k Nfabs(T[k][k]-T[k-1][k-1])1e-5 N Y 源程序 #includestdio.h#includestdlib.h#includemath.h#define a 0 //区间下限#define b 1 //区间上限#define N 20 //最大加速次数double f(double x){return(pow(x,1.5));}void main(){double h=b-a;int i,j,k,m,n;double T[N][N],s;T[0][0]=h*((*f)(a)+(*f)(b))/2.0; for(k = 1;kN;k++) {h = (b-a) * pow(0.5,k); n = (int)pow(2,k-1);s=0.0; for(m=1;m=n;m++) { s+=(*f)(a+(2*m-1)*h); } T[k][0] = 0.5*T[k-1][0] + h*s; }for(k=1;kN;k++){for(m=1;m=k;m++){T[k][m]=pow(4,m)/(pow(4,m)-1)*T[k][m-1]-1/(pow(4,m)-1)*T[k-1][m-1];}if(fabs(T[k][k]-T[k-1][k-1])1e-5){break;}}printf(R=%.6lf\n,T[k][k]);printf(k=%d \n,k);printf(T:\n);for(i = 0;i=k;i++) { for(j = 0;j=i;j++) printf(%.6lf ,T[i][j]); printf(\n); }}程序运行结果分析如上所示的结果与课本中求得的结果完全一样，表明程序编写正确，且符合要求，事实上，只要再将所求值的精度设置得更小，则所求的结果将更加准确，最终将无限接近于标准值，由上表也可以看出用龙贝格积分法求函数的积分值在精度比较低的情况下就能求到很准确的值！实验体会通过这次实验上机，使我更进一步了解了龙贝格法的计算思想，其在精度上很有保证，收敛较快，是解积分问题的有效方法。