Вычислить интеграл от функции f (x) на промежутке [a, b] методом прямоугольников

Доброго времени суток, помогите, пожалуйста. Не понимаю в чем ошибка.
Вычислить определенный интеграл от функции f (x) на промежутке [a, b] методом прямоугольников.
f(x)=(x/(sqrt(x)+4))*atan(x)+1 a=-2, b=2
Число разбиений n = 10; 40; 160; 640.
Мой код:
#include <iostream>
#include <math.h>
#include <stdlib.h>
using namespace std;

int main()
{
int i, j, n, m;
//создаем массив
cout << "Number of equations: ";
cin >> n;
cout << "Number of variables: ";
cin >> m;
m += 1;
float **matrix = new float *[n];
for (i = 0; i<n; i++)
matrix[i] = new float[m];

//инициализируем

for (i = 0; i<n; i++)

for (j = 0; j<m; j++)
{
cout << " Element " << "[" << i + 1 << " , " << j + 1 << "]: ";

cin >> matrix[i][j];
}

//выводим массив
cout << "matrix: " << endl;
for (i = 0; i<n; i++)
{
for (j = 0; j<m; j++)
cout << matrix[i][j] << " ";
cout << endl;
}
cout << endl;

//Метод Гаусса
//Прямой ход, приведение к верхнетреугольному виду
float tmp;
int k;
float *xx = new float[m];

for (i = 0; i<n; i++)
{
tmp = matrix[i][i];
for (j = n; j >= i; j--)
matrix[i][j] /= tmp;
for (j = i + 1; j<n; j++)
{
tmp = matrix[j][i];
for (k = n; k >= i; k--)
matrix[j][k] -= tmp*matrix[i][k];
}
}
/*обратный ход*/
xx[n - 1] = matrix[n - 1][n];
for (i = n - 2; i >= 0; i--)
{
xx[i] = matrix[i][n];
for (j = i + 1; j<n; j++) xx[i] -= matrix[i][j] * xx[j];
}

//Выводим решения
for (i = 0; i<n; i++)
cout << xx[i] << " ";
cout << endl;

delete[] matrix;
system("pause");
return 0;
}
ответы: -3.628, 6.709, 4.317, -0.0402, -0.266

0

Сообщение от Temerest a=-2, b=2

брал промежуток a = 0, b = 2, ибо x не может принимать отрицательных значений из-за квадратного корня (если подразумеваем получить вещественное число)

1