Код: Выделить всё
#include <vcl>
#include <iostream>
#include <conio>
using namespace std;
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()
{
system("pause");
}
1 -ый интеграл какой вычислить нужно (1/sqrt(0.5*pow(x,2)+1.5))
2 - ой интеграл (x+1)*cos(x*x)