function z=f3(t,x)
% f(t,x)=-x+2 x^2 exp(-t)
% x'=f, x(0)=1/3, T=1
% Soluzione: x=1/(2 exp(t) + exp(-t))
% Questa funzione serve anche a confrontare numero di passi 
% tra RK a passo fisso e ode45
z= - x + 2* x^2 *exp(-t);

