#include <iostream>
#include <cmath>

using namespace std;

// Definicja funkcji
double f(double x)
{
    return sin(x);
}

// Funkcja obliczająca całkę metodą trapezów
double trapezy(double a, double b, int n)
{
    double h = (b - a) / n;
    double suma = 0.0;

    suma += (f(a) + f(b)) / 2.0;

    for (int i = 1; i < n; i++)
    {
        suma += f(a + i * h);
    }

    return suma * h;
}

int main()
{
    double a = 0.0;
    double b = 3.14;

    int n1 = 10;
    int n2 = 100;

    double s1 = trapezy(a, b, n1);
    double s2 = trapezy(a, b, n2);

    cout << "s1 (n = 10)  = " << s1 << endl;
    cout << "s2 (n = 100) = " << s2 << endl;

    return 0;
}
