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

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <conio.h>
#include <iostream>
 
using namespace std;
 
double integralsimpson (float a, float b, float e);
 
int main (void)
{
   setlocale(LC_ALL, "RUS");
   float a1, b1, a=0, b=0, e=0;
   float epsilon;
   cout<<"Введите a "<<endl;
   cin>>a1;
   cout<<"Введите b "<<endl;
   cin>>b1;
   cout<<"Введите e "<<endl;
   cin>>epsilon;
   integralsimpson(a,b,e);
   return 0;
}
 
 
float f (float x)
{
   float result;
   result=log (x);
   return result;
}
 
double integralsimpson (float a, float b, float e)
{
   float result;
   float h;
   float s;
   float s1;
   float s2;
   float s3;
   float x;
   s2 = 1;
   h = b-a;
   s = f (a) +f (b);
   do
     {
       s3 = s2;
       h = h/2;
       s1 = 0;
       x = a+h;
       do
          {
            s1 = s1+2*f (x);
            x = x+2*h;
          }
       while (x<b);
       s = s+s1;
       s2 = (s+s1) *h/3;
       x = fabs (s3-s2) /15;
     }
   while (x>e);
   result = s2;
   cout<<endl<<"Интеграл равен: "<<result<<endl;
   return result;
}

double integralsimpson (float a, float b, float e);
 
int main (void)
{
   setlocale(LC_ALL, "RUS");
   float a=0, b=0, e=0;
   cout<<"Введите a "<<endl;
   cin>>a;
   cout<<"Введите b "<<endl;
   cin>>b;
   cout<<"Введите e "<<endl;
   cin>>e;
   integralsimpson(a,b,e);
   return 0;
}
 
 
float f (float x)
{
   float result;
   result=log (x);
   return result;
}
 
double integralsimpson (float a, float b, float e)
{
   float result;
   float h;
   float s;
   float s1;
   float s2;
   float s3;
   float x;
   s2 = 1;
   h = b-a;
   s = f (a) +f (b);
   do
     {
       s3 = s2;
       h = h/2;
       s1 = 0;
       x = a+h;
       do
          {
            s1 = s1+2*f (x);
            x = x+2*h;
          }
       while (x<b);
       s = s+s1;
       s2 = (s+s1) *h/3;
       x = fabs (s3-s2) /15;
     }
   while (x>e);
   result = s2;
   cout<<endl<<"Интеграл равен: "<<result<<endl;
   system("pause");
   return result;
}

float f (float x)
{
   float result;
   result=x/sin(x);
   return result;
}