function I=Simpson38(x,y)
N = length(x); % Tamaño del vector x
n = N-1;  % Número de segmentos
a = x(1);  % Limite inferior
b = x(N);  % Limite superior
s3 = 0; sn3 = 0; % inicializando sumatorias
for i=2:n
    if (i-1)/3 == ceil((i-1)/3)
        s3 = s3 + y(i); % indices multiplos de 3
    else
        sn3 = sn3 + y(i);% indices no multiplos de 3
    end
end
I=3*(b-a)*(y(1) + y(N) + 3*sn3 + 2*s3)/(8*n);