function y = trapecio(F,a,b,n)

//
// Integrales usando regla trapezoidal
//

y = 0;
x = a;
h = (b-a)/n;

for i = 1:n
   y = y + (h/2)*(feval(x,F)+feval(x+h,F));
   x = x + h;
end
