% Implementering av Simpsons metode.
% 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=[*];

% Integranden
funk=@(x) [**];

% Oppdeling (kontrollerer at n er et partall)
n=input('Gi oppdelinga n: ');
if round(n/2) ~= n/2
    disp('n må være et partall')
    return
end
h=[**];                  % Skrittlengda

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

% Oddetallsbidrag:
for i=1:2:[**];
    xi=a+i*h;
    S=S+h/3*4*funk(xi);
end

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

% Skriver svaret til skjerm
S