% Implementering av trapesmetoden for numerisk integrasjon.
% Integrasjonsgrensene a og b og integranden funk
% blir gitt heilt i toppen av skriptet.
% For å gi oppdelinga n brukar vi input-funksjonen

% Integrasjonsgrenser
a=-1;
b=1;

% Integranden
funk=@(x) x^6;

% Oppdeling
n=input('Gi oppdelinga n: ');
h=(b-a)/n;

% Bidrag fra endane
T=h/2*(funk(a)+funk(b));

% Resten av bidraga
for i=1:(n-1)
    xi=a+i*h;
    T=T+h*funk(xi);
end

% Skriv svaret til skjerm
T