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

#include "stdafx.h"
#include <iostream>
#include <conio.h>
using namespace std;
void glavnaja();
double summa(double);
double S1,v,S2,l=20,h,k,e=0.000006;
 
 
 
int main()
{
    
    glavnaja();
    cout << v << endl;
    _getch();
    return 0;
}
void glavnaja()
{
    do
    {
        cout << "ja rabotaju " << endl;
        v = summa(l);
        l = l * 2;
        S1 = v;
        
        v =summa(l);
        S2 = v;
        
 
    } while (abs(S1 - S2) >= e);
 
}
double summa(double p)
{
 
    double v,S=0, N, c = 0, x = 0.1, i, t, s,  a = 0.1, b = 1.2;
    N = p;
    h = (b - a) / N;
 
    
        s = 0.1;
        t = 1.4;
        
        while (x <  b ){
 
            S = (t * atan(pow(x, 0.2))) / x*x + s*sin(x);
            c = c + S;
            v = c / 2 * h;
            k = k + v;
            x = x + h;
        }
 
        
        
    
 
    return k;
 
}