#include <iostream>

#include <cstdlib>

using namespace std;

double Funkcja(double x)

{

return x*x-x+1;

}

int main()

{

double eps=0.000001,k,l,przedzial,pole;

cout<<"Podaj poczatek przedzialu calkowania"<<endl;

cin>>k;

cout<<"Podaj koniec przedzialu calkowania"<<endl;

cin>>l;

if(k>l)

cout<<"pierwszy przedzial musi byc mniejszy od drugiego!"<<endl;

przedzial = l-k;

int iloscTrapezow=0;

iloscTrapezow=(int)(przedzial/eps);

int i;

for(i=0; i<iloscTrapezow;++i)

pole +=((Funkcja(k+i*eps)+Funkcja(k+(i+1)*eps))/2.0)*eps;

cout<<"Metoda trapezow: "<<pole<<endl;

system("pause");

return 0;

}