function erg=mittel(a,b,n,f)
h=(b-a)/n;
tmp=0;
for j=1:n
    mop=a+(j-0.5)*h;
    tmp=tmp+f(mop);
end
    erg=h*tmp;
end