function z = fun(t,y)
z = y/t-(y/t)^2;