function y = f2(t,x)

y = -x.**2./(1+3.*exp(-t).*x);

end