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