#include <iostream>
#include <cmath>
 
using namespace std;
 
double pole_trapezow(double a, double b, int n) {
    double h = (b - a) / n;
    double suma = sin(a) + sin(b);
 
    for (int i = 1; i < n; i++) {
        double x = a + i * h;
        suma += 2 * sin(x);
    }
 
    return (h / 2) * suma;
}
 
int main() {
    double a = 0.0;
    double b = 3.14;
 
    int n1 = 10;
    int n2 = 100;
 
    cout << "Pole metoda trapezow dla n = 10:  "
         << pole_trapezow(a, b, n1) << endl;
 
    cout << "Pole metoda trapezow dla n = 100: "
         << pole_trapezow(a, b, n2) << endl;
 
    return 0;
}
 

I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKZG91YmxlIHBvbGVfdHJhcGV6b3coZG91YmxlIGEsIGRvdWJsZSBiLCBpbnQgbikgewogICAgZG91YmxlIGggPSAoYiAtIGEpIC8gbjsKICAgIGRvdWJsZSBzdW1hID0gc2luKGEpICsgc2luKGIpOwoKICAgIGZvciAoaW50IGkgPSAxOyBpIDwgbjsgaSsrKSB7CiAgICAgICAgZG91YmxlIHggPSBhICsgaSAqIGg7CiAgICAgICAgc3VtYSArPSAyICogc2luKHgpOwogICAgfQoKICAgIHJldHVybiAoaCAvIDIpICogc3VtYTsKfQoKaW50IG1haW4oKSB7CiAgICBkb3VibGUgYSA9IDAuMDsKICAgIGRvdWJsZSBiID0gMy4xNDsKCiAgICBpbnQgbjEgPSAxMDsKICAgIGludCBuMiA9IDEwMDsKCiAgICBjb3V0IDw8ICJQb2xlIG1ldG9kYSB0cmFwZXpvdyBkbGEgbiA9IDEwOiAgIgogICAgICAgICA8PCBwb2xlX3RyYXBlem93KGEsIGIsIG4xKSA8PCBlbmRsOwoKICAgIGNvdXQgPDwgIlBvbGUgbWV0b2RhIHRyYXBlem93IGRsYSBuID0gMTAwOiAiCiAgICAgICAgIDw8IHBvbGVfdHJhcGV6b3coYSwgYiwgbjIpIDw8IGVuZGw7CgogICAgcmV0dXJuIDA7Cn0K