function  T = trapezrule(f,a,b,m)
% Integración trapecio compuesto
 x = linspace(a,b,m+1);             % subintervalos
 T = (f(a)+f(b))/2;                 % comienzo y final
 for  i = 1:m-1
    T = T + f(x(i+1));              % interior intervalo
 end
 T = T*(b-a)/m;                     % multiplica por h
end

function f = intfun_1(x)
% Función que se integra
 f = 1./(1+x);
end
