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