Вычисление интеграла

void integral()
{
    float a, b;
    float dx;
    float s;
    int n;
    float x;
    float y;
    int i;
    SetColor(12, 0); printf("Приближенное вычисление интеграла (f(x)=(log(1 + x) / (1 + x*x))*dx");
    printf("\n\n");
    SetColor(15, 0); printf("Используем метод прямоугольников\n");
    printf(" \n\n");
    printf("Для вычисления введите следующие значения:\n");
 
    printf("Нижняя граница интервала:\n");
    scanf_s("%f", &a);
    printf("Верхняя граница интервала:\n");
    scanf_s("%f", &b);
    printf("\nПриращение аргумента:\n");
    scanf_s("%f", &dx);
 
    n = (b - a) / dx + 1;
    x = a;
    s = 0;
 
    for (i = 1; i <= n; i++)
    {
        y = (log(1 + x) / (1 + x*x))*dx;
        s += y*dx;
        x += dx;
    }
    printf("\n\nЗначение интеграла: %6.3f", s);
    _getch();
    
}

#include <iostream>
 
using namespace std;
double F(double x)
{
    return 1 / (2 + x*x);
}
 
double methodTrap(double a, double b, int n, double(*f)(double))
{
    double h = (b - a) / n, sum = 0;
    for (int i = 0; i < n; i++)
    {
        sum += 0.5 * h * (f(a + i * h) + f(a + (i + 1) * h));
    }
    return sum;
}
 
int main()
{
    int n;
    cout << "n=";
    cin >> n;
    double a, b;
    cout << "a=";
    cin >> a;
    cout << "b=";
    cin >> b;
 
    cout << methodTrap(a, b, n, F) << "\n";
    system("pause");
}