#include <iostream>
#include <cmath>
using namespace std;
double trapezy(double a, double b, int n)
{
double dx = (b - a) / n;
double s = 0.0;
double x = a;
for (int i = 1; i <= n; i++)
{
s = s + dx * sin(x);
x = x + dx;
}
return s;
}
int main()
{
double a1 = 0.0;
double b1 = 3.14;
int n1 = 10;
double a2 = 0.0;
double b2 = 3.14;
int n2 = 100;
double s1 = trapezy(a1, b1, n1);
double s2 = trapezy(a2, b2, n2);
cout << "s1 (n = 10) = " << s1 << endl;
cout << "s2 (n = 100) = " << s2 << endl;
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKZG91YmxlIHRyYXBlenkoZG91YmxlIGEsIGRvdWJsZSBiLCBpbnQgbikKewogICAgZG91YmxlIGR4ID0gKGIgLSBhKSAvIG47CiAgICBkb3VibGUgcyA9IDAuMDsKICAgIGRvdWJsZSB4ID0gYTsKCiAgICBmb3IgKGludCBpID0gMTsgaSA8PSBuOyBpKyspCiAgICB7CiAgICAgICAgcyA9IHMgKyBkeCAqIHNpbih4KTsKICAgICAgICB4ID0geCArIGR4OwogICAgfQoKICAgIHJldHVybiBzOwp9CgppbnQgbWFpbigpCnsKICAgIGRvdWJsZSBhMSA9IDAuMDsKICAgIGRvdWJsZSBiMSA9IDMuMTQ7CiAgICBpbnQgbjEgPSAxMDsKCiAgICBkb3VibGUgYTIgPSAwLjA7CiAgICBkb3VibGUgYjIgPSAzLjE0OwogICAgaW50IG4yID0gMTAwOwoKICAgIGRvdWJsZSBzMSA9IHRyYXBlenkoYTEsIGIxLCBuMSk7CiAgICBkb3VibGUgczIgPSB0cmFwZXp5KGEyLCBiMiwgbjIpOwoKICAgIGNvdXQgPDwgInMxIChuID0gMTApICA9ICIgPDwgczEgPDwgZW5kbDsKICAgIGNvdXQgPDwgInMyIChuID0gMTAwKSA9ICIgPDwgczIgPDwgZW5kbDsKCiAgICByZXR1cm4gMDsKfQ==