#include <vcl>
#include <iostream>
#include <conio>
#include <math>
using namespace std;
 
double Simpson(double (*Fx)(double))
{
        double I;
        double h;
        double x;
        double a,b; // a- нижня грацница интегрирования, b - верхняя граница интегрирования
        int m; // точность интегирования ( кол-во промежутков разбиения )
 
 // Этот код - заполнение данных для консоли
 // вам возможно прийдется его преписать, если вы делаете через GUI
 
        cout << "\n Input integral atributes:\n" << endl;
        cout << "\t-> Enter a = ";
        cin >> a;
 
        cout << "\t-> Enter b = ";
        cin >> b;
 
        cout << "\t-> Enter m = ";
        cin >> m;
 
 // вычисление интеграла
        h=(b-a)/(m-1);
        x = a;

I = Fx(a) + Fx(b);
        int n=0;
 
        while (n < m-2)
        {
                x = x+h;
                if (n%2==0) I=I+4*Fx(x);
                        else  I = I+2*Fx(x);
                n++;
        }
 
}
 
main()
{
 return I*h/3;
}