a=1;
b=4;
n=5;
h=(b-a)/(n-1);
p=a;
k=a+h;
calka=(h/2)*(f(p)+f(k));
for i=2:n
p=k;
k=k+h;
z=f(p);
x=f(k);
calka=calka+(h/2)*(z+x);
end
poch2=2;
reszta=(-1/12)*n*(h^3)*poch2;

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