[tt,xx] = ode45(@finducts,[0:1:5],2.1);
format long
[tt,xx]

function xp = finducts(t,x)
%Nonlinear function
xp = 20*x^2/(16+x^2) - 2*x;
end