function [I]=comptrapz(a,b,F,N)

Dx=(b-a)/N;
I=0;
for i=1:N-1
    I=I+Dx*F(a+i*Dx);
end
I=I+(F(a)+F(b))*Dx/2;

end