trapez1.m

% Integration mit der Trapezregel
%   anpassen: 
%       Grenzen a, b
%       Funktion (statt sin)
%       Zahl der Intervalle N

N = 100;
a = 0;
b = pi/2;
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))