Найти приближенное значение интеграла

#include <math.h>
#include <iostream>
 
using namespace std;
 
// численно интегрируемая функция
double F(double x)
{
    int n=14;
    return sin(0.02*n +1.5*x)/(1.4+cos(1.2*n-0.3*x));
}
 
 
double Rectangular(double a, double b, int n)
{
    double sum = 0.0;
    double h = (b - a) / n ;
    for (int i = 0; i < n ; i++)
    {
        sum += h * F(a + i * h);
    }
    return sum;
}
 
double Trapezoidal(double a, double b, int n)
{
    double sum = 0.0;
    double h = (b - a) / n;
    for (int i = 0; i < n; i++)
    {
        sum += 0.5 * h * (F(a + i * h) + F(a + (i + 1) * h));
    }
    return sum;
}
 
double Simpson(double a, double b, int n)
{
    double sum = 0.0;
    double h = (b - a) / n;
    if (n % 2 != 0)
    {
        for (int i = 0; i < n - 1; i += 2)
        {
            sum += h * (F(a + i * h) + 4 * F(a + (i + 1) * h) + F(a + (i + 2) * h)) / 3;
        }
    }
    else
    {
        sum = 3 * h * (F(a) + 3 * F(a + h) + 3 * F(a + 2 * h) + F(a + 3 * h)) / 8;
        for (int i = 3; i < n - 1; i += 2)
        {
            sum += h * (F(a + i * h) + 4 * F(a + (i + 1) * h) + F(a + (i + 2) * h)) / 3;
        }
    }
    return sum;
}
 
int main()
{
    int n = 1000;
    double a = 0.1;
    double b = 1.7;
    
    cout << "Testing Rectangular method" << endl;
    cout << "Result = " << Rectangular(a, b, n) << endl;
 
    cout << "Testing Trapezoidal method" << endl;
    cout << "Result = " << Trapezoidal(a, b, n) << endl;
 
    cout << "Testing Simpson method" << endl;
    cout << "Result = " << Simpson(a, b, n) << endl;
 
    system("Pause");
    return 0;
}