龙贝格+实验报告.docxVIP

学院：数学科学学院　　　专业：实验名称成绩姓名年级学号实验三 龙贝格求积分实验目的 用算法实现龙贝积分运算，求解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 = (