Найти значение интеграла 5*x*Sin(x) на промежутке [0,π] с точностью ε=0,001, используя метод средних прямоугольников

}while(fabs(f) > eps && iter < 20000);
  cout << iter << " iterations\n\n";

#include <stdio.h> 
#include <math.h> 
 
double f(double x)
{
    return 5 * x * sin(x);
}
 
int main() 
{
    double h = 0.019;
    const double pi = 4 * atan(1.0);
    int n = (int)(pi / h);
    n = n + 1;
    h = pi / n;
    double ApprResult = 0;
    int i;
    double xi;
    for (i = 0; i < n; i++)
    {
        xi = i * h;
        ApprResult = ApprResult + h * (f(xi) + f(xi + h)) / 2;
    }
    printf("Approach of Result = %f\n", ApprResult);
 
    double exactResult = 5 * sin(pi) - 5 * pi * cos(pi);
    printf("Exact Result = %f\n", exactResult);
    double eRR = fabs(exactResult - ApprResult);
    printf("Error = %f", eRR);
    return 0;
}

double func(double x)
{
    return 5*x*sin(x);
}
double integral(double x)
{
    return 5*sin(x)-5*x*cos(x);
}
void main()
{
    double xa=0,xb=M_PI;
    double Sint=integral(xb)-integral(xa);
    cout<<"Sintg="<<Sint<<endl;
    double eps=0.001,Srct;
    int n=0;
    do
    {
        double dx=(xb-xa)/++n;
        Srct=0;
        double x=xa+dx/2;
        for(int i=0; i<n; x+=dx,i++)
            Srct+=func(x);
        Srct*=dx;
    }while(fabs(Sint-Srct)>=eps);
    cout<<"Srect="<<Srct<<"  n="<<n<<endl;
    system("pause");
}