换手机号微信还能用吗:romberg积分 c++程序

里面有一个错误 望高手更正
#include <iostream>
#include <cmath>

using namespace std;

const double E = 2.718281828;
const double PI = 3.1415926;
const double EPS = 1e-5;
const int ROW_COL = 100;

double matrix[ROW_COL][ROW_COL];

void InitialMatrix()
{
int i;
for (i=0; i<ROW_COL; i++)
{
int j;
for (j=0; j<ROW_COL; j++) matrix[i][j] = 0.0;
}
return;
}

inline double function(double x) //此处修改被积函数
{
// return 2 / sqrt(PI) * (pow(E, -x));
return x * sqrt(1 + pow(x, 2));
}

void PrintMatrix(int x)
{
int i;
for (i=0; i<x; i++)
{
int j;
for (j=0; j<x; j++) cout << matrix[i][j] << " ";
cout << endl;
}
return;
}

int Romberg(double a, double b) // a为积分下界, b为上界
{
double h = b - a;
matrix[0][0] = (h / 2) * (function(b) + function(a));

int k = 0;
int divideNum = 2;
int i;
do
{
k++;
double sum = 0;

double stepLen = h / (divideNum);
double temp = a;
for (i=0; i<divideNum - 1; i++)
{
temp += stepLen;
sum += function(temp);
}
divideNum *= 2;

matrix[k][0] = stepLen / 2 * (function(a) + function(b) + 2 * sum);

for (i=0; i<k; i++)
{
double reserve = pow(4, (i + 1));
matrix[k][i + 1] = (reserve / (reserve - 1)) * matrix[k][i] - (1 / (reserve - 1)) * matrix[k - 1][i];
}

} while (fabs(matrix[k][k] - matrix[k - 1][k - 1]) > EPS);
return k;
}

int main()
{
cout.setf(ios::fixed,ios::floatfield);
cout.precision(7);

InitialMatrix();
int x = Romberg(0, 3);
PrintMatrix(x + 1);
return 0;
}