function S=SimpsonFunk(funk,a,b,n)

% Implementering av Simpsons metode.
% Integranden funk, grensene a og b, og oppdelinga n 
% blir gitt som input.

% Skrittlengda
h=(b-a)/n;

% Bidrag fra endane
S=h/3*(funk(a)+funk(b));

% Oddetalsbidrag:
for i=1:2:(n-1);
    xi=a+i*h;
    S=S+h/3*4*funk(xi);
end

% Partalsbidrag
for i=2:2:(n-2)
    xi=a+i*h;
    S=S+h/3*2*funk(xi);
end
