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

% Integrasjonsgrenser
a=1;
b=3;

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

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

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

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

% Skriver svaret til skjerm
T