a=1;
b=4;
h=b-a;
n=5;
z=f(a);
x=f(b);
calka=(h/2)*(z+x);

poch2=2;
reszta=(-1/12)*(h^3)*poch2;

function y=f(x)
y=(x*x)-x+1;
end