#include <iostream>
#include <cmath>
 
using namespace std;
 
double f(double x) {
    return sqrt(1 - x * x);
}
 
int main() {
    int n = 100000;
    double a = -1.0;
    double b = 1.0;
    double h = (b - a) / n;
 
    double suma = 0.5 * (f(a) + f(b));
 
    for(int i = 1; i < n; i++) {
        double x = a + i * h;
        suma += f(x);
    }
 
    double pi = 2 * h * suma;
 
    cout << "Przyblizona wartosc pi (metoda trapezow): " << pi << endl;
 
    return 0;
}

I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKZG91YmxlIGYoZG91YmxlIHgpIHsKICAgIHJldHVybiBzcXJ0KDEgLSB4ICogeCk7Cn0KCmludCBtYWluKCkgewogICAgaW50IG4gPSAxMDAwMDA7CiAgICBkb3VibGUgYSA9IC0xLjA7CiAgICBkb3VibGUgYiA9IDEuMDsKICAgIGRvdWJsZSBoID0gKGIgLSBhKSAvIG47CgogICAgZG91YmxlIHN1bWEgPSAwLjUgKiAoZihhKSArIGYoYikpOwoKICAgIGZvcihpbnQgaSA9IDE7IGkgPCBuOyBpKyspIHsKICAgICAgICBkb3VibGUgeCA9IGEgKyBpICogaDsKICAgICAgICBzdW1hICs9IGYoeCk7CiAgICB9CgogICAgZG91YmxlIHBpID0gMiAqIGggKiBzdW1hOwoKICAgIGNvdXQgPDwgIlByenlibGl6b25hIHdhcnRvc2MgcGkgKG1ldG9kYSB0cmFwZXpvdyk6ICIgPDwgcGkgPDwgZW5kbDsKCiAgICByZXR1cm4gMDsKfQ==