function I = simpson_1(f,a,b,n)
% Integración de f en [a,b] con la regla de Simpson en n+1
% puntos igualmente espaciados
 h = (b-a)/n;
 xi = a:h:b;
 I = (h/3)*(f(xi(1))+2*sum(f(xi(3:2:end-2)))+4*sum(f(xi(2:2:end)))+f(xi(end)));
end

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