Интеграл

double func(const double x0)
{
    return Math::Sin(x0)*Math::Sin(x0);
}
 
double Simpson(const double a, const double b, const int n)
{
    double h = (b - a)/n;
    double s = 0;
    double x = a;
    for (int i = 1; i <= n; i++)
    {
        s += func(x) + 4*func(x + h/2) + func(x + h);
        x += h;
    }
    s *= h/6;
    return s;
}

#include "iostream.h"
#include "math.h"
 
double func(const double x0)
{
        return sin(x0)*sin(x0);
}
 
double Simpson(const double a, const double b, const int n)
{
        double h = (b - a)/n;
        double s = 0;
        double x = a;
        for (int i = 1; i <= n; i++)
        {
                s += func(x) + 4*func(x + h/2) + func(x + h);
                x += h;
        }
        s *= h/6;
        return s;
}
 
void main(void)
{
     double a = 0;
     double b = 3.14/2;
     int n = 100;//количество отрезков на которое разбивается интервал, чем больше тем выше точность
     double s = Simpson(a, b, n);
     std::cout << "S = "<< s;
     std::cin.get();
}