#include <iostream>
#include <cmath>
using namespace std;
int main()
{
double a = 0.0;
double b = 1.57;
int n = 10;
double h = (b - a) / n;
double S = 0.0;
S += sin(a) + sin(b);
for (int i = 1; i < n; i++)
{
double x = a + i * h;
S += 2 * sin(x);
}
S *= h / 2;
cout << "Przyblizona wartosc pola (metoda trapezow) S = " << S << endl;
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKaW50IG1haW4oKQp7CiAgICBkb3VibGUgYSA9IDAuMDsKICAgIGRvdWJsZSBiID0gMS41NzsKICAgIGludCBuID0gMTA7CgogICAgZG91YmxlIGggPSAoYiAtIGEpIC8gbjsKICAgIGRvdWJsZSBTID0gMC4wOwoKICAgIFMgKz0gc2luKGEpICsgc2luKGIpOwoKICAgIGZvciAoaW50IGkgPSAxOyBpIDwgbjsgaSsrKQogICAgewogICAgICAgIGRvdWJsZSB4ID0gYSArIGkgKiBoOwogICAgICAgIFMgKz0gMiAqIHNpbih4KTsKICAgIH0KCiAgICBTICo9IGggLyAyOwoKICAgIGNvdXQgPDwgIlByenlibGl6b25hIHdhcnRvc2MgcG9sYSAobWV0b2RhIHRyYXBlem93KSBTID0gIiA8PCBTIDw8IGVuZGw7CgogICAgcmV0dXJuIDA7Cn0=