function result = simpsonsMetode(f, start, stopp, delintervall)
    xs = []; % liste for grovparten av Simpsons metode
    dx = (stopp - start)/delintervall;
    for i=(start + dx):2*dx:(stopp - dx) % konstruksjon av dx/3 ( ... 4f1 + 2f2 + 4f4 ... 4fn-1)
        xs = [xs, 4*f(i), 2*f(i + dx)];
    end
    xs = [f(start), xs, f(stopp)]; % (f(a) + 4f1 + 2f2 + ... + 2fn-1 + f(b)
    result = sum(dx/3 * xs); % dx/3 (f(a) + 4f1 + 2f2 + ... + 2fn-1 + f(b)
end