function I=Simpson13(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
I=(b-a)*(y(1)+y(N)+4*sum(y(2:2:n))+...
    2*sum(y(3:2:n-1)))/(3*n);