trapez1.m

% Integration mit der Trapezregel
%   anpassen: 
%       Grenzen a, b
%       Funktion (statt sin)
%       Zahl der Intervalle N
a = 0;
b = pi/2;
N = 100;

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