Нахождение интеграла методом Симпсона

double Simpson(double (*Fx)(double))
{
        double I;
        double h;
        double x;
        double a,b;
        int m;
        
 
        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++;
        }
 
        return I*h/3;
}

double sinus(double x)
{
       /* формула от какой найти интеграл */
       return x;
}

double sinus(double x)
{
       /* формула от какой найти интеграл */
       return x;
}

double f(double x)
{
    double r = sin(x);
    r = 1. - r*r/4.;
    return r;
}

double I2 = 0, I4 = 0;
for(int k = 2; k<m;k+=2)
{
    I2 += Fx(a+k*h);
    I4 += Fx(a+(k-1)*h);
}
I = Fx(a)+Fx(b)+4*I4+2*I2;
I *= h/3;

I4 = Fx(a+h);
for(int k = 2; k < m; k += 2 )
{
    I4 += Fx(a+(k+1)*h);
    I2 += Fx(a+k*h);
}

double func(double x )
{
    // function code
}
double simpson( double (*Fx)(double), double a, double b, int m);
 
int main()
{
    // ввод параметров
    // ....................
 
    // вычисление интеграла
 
    I = simpson( &func, a, b, m );
 
    // вывод результата
    // ................
}
 
double simpson( double (*Fx)(double), double a, double b, int m)
{
    // здесь код только вычислений
}

double somefunc(double x)
{
    return x*x + 3*x;
}

//---------------------------------------------------------------------------
 
#include <vcl.h>
#include <math.h>
#include <iostream.h>
#include <conio.h>
#pragma hdrstop
 
using namespace std;
 
//Ïîäûíòåãðàëüíàÿ ôóíêöèÿ
double func(double x)
{
    double r = sin(x);
    r = 1. - r*r/4.;
    return r;
}
double simpson( double (*Fx)(double), double a, double b, int m);
 
#pragma argsused
int main(int argc, char* argv[])
{
   double I;
   double a,b;
   int m;
 
   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;
 
   // âû÷èñëåíèå èíòåãðàëà
 
   I = simpson( &func, a, b, m );
 
   // âûâîä ðåçóëüòàòà
   cout << endl << "      I = " << I;
   getch();
}
 
//---------------------------------------------------------------------------
// Ìåòîä ñèìïñîíà
double simpson( double (*Fx)(double), double a, double b, int m)
{
 
   double h;
   h=(b-a)/m;
 
   double I, I2 = 0, I4 = 0;
   I4 = Fx(a+h);
   for(int k = 2; k < m; k += 2 )
   {
       I4 += Fx(a+(k+1)*h);
       I2 += Fx(a+k*h);
   }
   I = Fx(a)+Fx(b)+4*I4+2*I2;
   I *= h/3;
 
   return I;
}
//---------------------------------------------------------------------------

std::cin.get();

system("pause");

double Fx(double x)                 //  функція
{ double s; 
s=sin(x);
return s;                                //   повертає значення "х"
}
void __fastcall TForm1::ToolButton1Click(TObject *Sender)
{
double a,b,h,n,s=0,s1=0,s2=0,x; // змінні
int i;
a=StrToFloat(Edit1->Text);        //  нижня межа інтеграла
b=StrToFloat(Edit2->Text);        //  верхня межа інтеграла
n=StrToInt(CSpinEdit1->Text);   //  кількість розбиттів
h=(b-a)/n;                              //  величина кроку
if(a==b)
ShowMessage("Невірні значення a i b");
{for(i=0;i<=n-1;i++)
{x=a+i*h;
if (i%2==0)
s2=s2+Fx(x);
else s1=s1+Fx(x);}
s=h/3*(Fx(a)+Fx(b)+4*s1+2*s2);
Form1->Memo1->Lines->Text="Інтеграл від функції X на проміжку["+FloatToStr(a)+";"+FloatToStr(b)+"] дорівнює "+FloatToStr(s);
}
}