Решение интеграла методом симпсона

#include "stdafx.h"
#include <iostream>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
 
using namespace std;
float f (float x)
{
   
   float result;
   result=(x*x)/(8+(x*x*x));
   return result;
}
float integralsimpson (float a, float b, float e,int &k);
 
int main (void)
{
    setlocale(LC_ALL,"rus");
   float a, b, e,result;
   int K=0;
    cout << "Введите нижнее значение интеграла:";
    cin >> a;
    cout << "Введите верхнее значение интеграла:";
    cin >> b;
    cout << "Введите точность:";
    cin >> e;
    result = integralsimpson(a,b,e,K);
    cout<<endl<<"Интеграл равен: "<<result<<endl;  
    cout<<"Количество шагов: "<<K<<endl;
    system("pause");
   return 0;
}
 
 
float integralsimpson (float a, float b, float e,int &k)
{
   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
     {
       k++;
       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; 
   return result;
}

result=(x*x)/(8+(x*x*x));

#include<iostream>
#include<cmath>
 
using namespace std;
 
double f(double x)
{
 return (x*x)/(8+(x*x*x));
}
 
double simpson(double a, double b, double e)
 {
 double y,x,h,I,I1,I2,I3;
 
 I2 = 1;
 
 h = b-a;
 I = f(a) + f(b);
 
 do
 {
  I3 = I2;
  h = h/2;
  I1 = 0;
  x = a + h;
  do
  {
   I1 = I1 + 2*f(x);
   x = x + 2*h;
  }
  while(x < b);
 
  I = I + I1;
  I2 = (I + I1)*h/3;
  x = fabs(I3 - I2)/15;
 }
 while(x > e);
 
 y = I2;
 return y;
}
 
int main()
{
 double a;
 double b;
 double e;
 
cout<<"vvedite a = ";
cin>>a;
cout << endl;
 
cout<<"vvedite b = ";
cin>>b;
cout << endl;
 
cout<<"vvedite pogreshnost = ";
cin>>e;
cout << endl;
 
cout<<"Integral Summa ravna  = " << simpson(a,b,e) <<endl;
return 0;
}