function [I]=trapznew(x,y)
nx=length(x);
ny=length(y);
if nx~=ny
 disp('x and y must be of the same length!')
 return
else
 n=nx;
end
avg_y=y(1:n-1)+diff(y)/2;
I=sum(diff(x).*avg_y);