trapez2.m

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));