function s=integriere(f,a,b,n) x=0:(n-1); h=(b-a)/n; x=a+h*x; s=0; for i=1:n s=s+f(x(i))*h; end end