а це царський сімпсон

// ConsoleApplication49.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include <iostream>
#include <cmath>
double f(double x)  // Ліві частини рівнянь
{
    double y;
    y = x;//           тут потрібно ввести вашу функцію

    return y;
}
double sym( int m, double x, double w)
{
    double s;
    s = 0;
    for (int i = 1; i <= m; i++)
    {
        s = s + f(x);
        x = x + w;

    }
    return s;
}

int main()
{
    double a,b,eps,N;
    std::cout << "vvedit mezi a ta b\n";
    std::cin >> a >> b;
    std::cout << "vvedit tochnist eps\n";
    std::cin >> eps;
    std::cout << "vvedit pochatkovui podil\n";//ну типу наближення(кількість поділів)
    std::cin >> N;

    double h = (b - a) / (N);
    double int_new, int_old = 0;
    int_new = h*(f( a) + f( b) + 4 * sym( N, a + h / 2, h) + 2 * sym( N - 1, a + h, h)) / 6;
    while (abs(int_new - int_old) > eps)
    {
        N = N + 1;
        int_old = int_new;
        h = (b - a) / N;
        int_new = h*(f( a) + f( b) + 4 * sym( N, a + h / 2, h) + 2 * sym( N - 1, a + h, h)) / 6;
    }
    std::cout << "znachennja integralu = " << int_new;
    system("pause");
}