Romberg积分算法代码.c

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总是指示待计算元素的前一个元素(同一行)