function I = trapez(func, a, b, N) % Integration mit der Trapezregel % Argumente % zu integrierende Funktion func % Grenzen a, b % Zahl der Intervalle N % Ergebnis % Integral von func von a bis b, % genaehert mit Trapezregel mit N Intervallen h = (b-a)/N; x = a:h:b; y = func(x); I = h*(0.5*y(1) + sum(y(2:N)) + 0.5*y(N+1));