Интеграл - C++ - Ответ 193842

cout << Simpson(f) << endl;

#include <iostream.h>
#include <math.h>
 
// интегрируемая функция
double f(double x)
{
    return 2*asin(x/2)-1/2*sin(4*asin(x/2));
}
 
// нахождение интегралла методом Симпсона
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;
}
 
int main()
{
    cout << Simpson(&f) << endl;
    cin.get(); // pause
    return 0;
}