[tt,xx] = ode45(@finduct,[0:1:5],3.1);
[tt,xx]

function xp = finduct(t,x)
%Nonlinear function
xp = 32.5*x^2/(30+x^2) - 2.5*x;
end