function I=Trapezoidal(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)+2*sum(y(2:n)))/(2*n);