#include <iostream>
#include <cmath>
 
using namespace std;
 
// Funkcja f(x) – NIEUJEMNA na [a, b]
double f(double x)
{
    return x * x;   // przykład
}
 
int main()
{
    // Dane wbudowane
    double a = 0.0;
    double b = 1.0;
    int n = 1000;
 
    double h = (b - a) / n;
    double s = 0.0;
 
    for (int i = 0; i < n; i++)
    {
        double x = a + (i + 0.5) * h; // środek przedziału
        s += f(x);
    }
 
    s *= h;
 
    cout << "Przyblizona wartosc pola: " << s << endl;
 
    return 0;
}
 

I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKLy8gRnVua2NqYSBmKHgpIOKAkyBOSUVVSkVNTkEgbmEgW2EsIGJdCmRvdWJsZSBmKGRvdWJsZSB4KQp7CiAgICByZXR1cm4geCAqIHg7ICAgLy8gcHJ6eWvFgmFkCn0KCmludCBtYWluKCkKewogICAgLy8gRGFuZSB3YnVkb3dhbmUKICAgIGRvdWJsZSBhID0gMC4wOwogICAgZG91YmxlIGIgPSAxLjA7CiAgICBpbnQgbiA9IDEwMDA7CgogICAgZG91YmxlIGggPSAoYiAtIGEpIC8gbjsKICAgIGRvdWJsZSBzID0gMC4wOwoKICAgIGZvciAoaW50IGkgPSAwOyBpIDwgbjsgaSsrKQogICAgewogICAgICAgIGRvdWJsZSB4ID0gYSArIChpICsgMC41KSAqIGg7IC8vIMWbcm9kZWsgcHJ6ZWR6aWHFgnUKICAgICAgICBzICs9IGYoeCk7CiAgICB9CgogICAgcyAqPSBoOwoKICAgIGNvdXQgPDwgIlByenlibGl6b25hIHdhcnRvc2MgcG9sYTogIiA8PCBzIDw8IGVuZGw7CgogICAgcmV0dXJuIDA7Cn0K