Решение интеграла методом трапеции С++

#include<iostream>
#include<conio.h>
#include<math.h>
using namespace std;
 
double f(double x) {return 3*x*x+sin(x);}
void main()
{
   const double e=0.0001;
   double a, b, s=0, h, s1;
   int n=10, i;
   cout<< " a= "; cin>>a;
   cout<< " b= "; cin>>b;
 
   h=(b-a)/n;
   s1=0; s=0; i=1;
   
   for (int i=1; i<n; i++)
   {s=s+f(a+i*h);}
 
   s=h*(s+(f(a)+f(b))/2;
 
   while (fabs(s1-s)>3*e)
   {n=n*2; h=h/2; s1=s;}
 
   cout<<"  s= ";
   _getch();
}

while (fabs(s1-s)>3*e)
{n=n*2; h=h/2; s1=s;}

while (fabs(s1-s)>3*e){
   //блаблабла подбираем злачение эн больше прежнего
  for (int i=1; i<n; i++) {
     //блаблабла интегрируем
  }
}

#include <iostream>
#include<math.h>
using namespace std;
 
double f(double x) {return 3*x*x+sin(x);}
int main()
{
   const double e=0.0001;
   double a, b, s=0, h, s1;
   int n=10;
   cout<< " a= "; cin>>a;
   cout<< " b= "; cin>>b;
   do
   {
       h=(b-a)/n;
       s1=s;
       s=0;
       for (int i=1; i<n; i++)
           s=s+f(a+i*h); 
       s=h*(s+(f(a)+f(b))/2);
       n=n*2; 
   }
   while (fabs(s1-s)>3*e);
   cout<<"  s= "<<s;
return 0;
}