Метод Симпсона три восьмых

#include <stdlib.h>
#include <conio.h>
#include <stdio.h>
#include <iostream>
#include <locale.h>
#include <math.h>
using namespace std;
float Y(double p);

int main()
{
        setlocale(LC_CTYPE,"Russian");
       
            double x, a, b, h, s,i;
            int n;
            cout<<"Отрезок интегрирования [a,b] --> (a) =\n";
            cin>>a;
           cout<<"Отрезок интегрирования [a,b] --> (b) =\n";
            cin>>b;
                       cout<<"число разбиений --> (n) =\n";
                       cin>>n;

           

            h = (b - a)/ (2*n);
            cout<<"H="<<h;
            s = 0; x = a + h;
            while (x < b)        
            {
                s = s + 2 * Y(x);
                x = x + h;
                s = s + 3 * Y(x);
                x = x + h;
            }
            s = (3*h)/ 8 * ( s + Y(a) + Y(b));
            cout<<"Интеграл = {0}"<<s;
           
}
float Y(double p)
{
return 1/(1+p*p);
}