Решить методом трапеций или Симпсона определенный интеграл

                           int N = 0;
        
                           double h = DBL_MAX;
        double eps = 0.01;
                           double a, b; 
                            a = 1;
                           b = 2;   
        double S1 = DBL_MAX;
        double S2 = DBL_MAX;
        do
        {
            
            S2 = S1;
            S1 = 0;
            N+=2;
            h = (a - b)/N;
            double x = a + h;
            for(int n = 1; n< N; n++)
            {
                n++;
                double x1 = (x-h);
                double x2 = (x+h);
                S1 += FormulaValue(x1) + 4*FormulaValue(x)+FormulaValue(x2);
                x +=2.*h;
            }
            S1 = S1*h/3.;
        }
        while((2*abs(S1-S2)/abs(S1+S2)) > eps);
        
             dblVal = (S1+S2)/2.;
    
    return dblVal;