Форум программистов [Powered by Invision Power Board]

Код

#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;
}