trapez.m

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