Вычисление определенного интеграла методом трапеций

#include <cstdlib> // для system
#include <iostream> 
#include <string>
 
using std::cout;
using std::cin;
using std::endl;
        int n,i;
        double a,b,h,x,s;
        double fun(double x);
int main()
{    
        setlocale(0, "");
        cout<<"Введите начало интервала a=";
        cin>>a;
        cout<<endl;
        cout<<"Введите конец интервала b=";
        cin>>b;
        cout<<endl;
        cout<<"Введите количество разбиений интервала n=";
        cin>>n;
        cout<<endl;
        //n=100;
        //a=0;
        //b=4;
        h=(b-a)/n;
        x=a;
        s=0;
        for (i = 0; i<n; i++) 
        {
            s=s+0.5*(fun(x)+fun(x+h))*h;
            x=x+h;
        }
        cout << "Интеграл="<< s <<endl;
        system("pause"); // Только для тех, у кого MS Visual Studio
}
 
double fun (double x)
{
    return (2*x+1)/sqrt(3+(pow(x,3)));
}