function I = Simpson(x, y)
    n = length(x);
    A = 0;
    for i = 1 : n - 2 : 2
        A =  A + (x(i+2) - x(i)) * (y(i) + y(i+2) + 4 * y(i+1));
    end
    I = A/6;
end