文件名称:Romberg积分算法代码.c
文件大小:2KB
文件格式:C
更新时间:2023-05-16 12:29:52
Romberg 积分算法 代码
double Romberg(double aa, double bb) { //aa,bb 积分上下限 int m, n;//m控制迭代次数, 而n控制复化梯形积分的分点数. n=2^m double h, x; double s, q; double ep; //精度要求 double *y = new double[MAXREPT];//为节省空间,只需一维数组 //每次循环依次存储Romberg计算表的每行元素,以供计算下一行,算完后更新 double p;//p总是指示待计算元素的前一个元素(同一行)