function T=litenTrapets(funktion,x0,xend,steg)
delta=xend-x0;
h=delta/steg;
x=x0:h:xend;
f=funktion(x);
T=h*(sum(f)-(f(1)+f(end))/2);