function erg2=trapez(a,b,n,f)
h=(b-a)/n;
tmp=0;

for j=1:n
    mop=a+(j-1)*h;
    tmp=tmp+f(mop);
end

erg2=(h/2)*(f(a)+2*tmp+f(b));

end