function y=f5(t,x)
% f(t,x)=-5 t x^2 + 5/t - 1/t^2
% La soluzione di x'=f(t,x), x(1)=1 e'
% x(t)=1/t (molto "smooth")
% Integrare su [1,25] con LMM Adams, con vari valori di k, e calcolare
% il massimo errore globale dell'approssimazione
%
% Esempio
% h   err(k=1)   err(k=2)  err(k=4)
% 0.2    0.03     inf        inf
% 0.1    0.009    0.009      inf
% 0.05   0.003    0.002      0.8
% 0.01   6e-5     1.e-6      1.e-6
% 0.005  3e-5     2.e-7      2.e-7
% 0.002  1e-5     4.e-8      4.e-8

% Si osserva che se k e' alto l'approssimazione e' migliore, ammesso
% che h sia abbastanza piccolo, da rimanere nella regione di stabilita'
% assoluta. Al contrario di RK, la regione di stabilita' assoluta si
% restringe se k cresce.
% Prova anche a calcolare la soluzione con ode113, ode15s (Matlab) o odebdi
% (Octave)
y = -5*t*x^2 + 5/t - 1/t^2;
