syms y(t)
y(t) = dsolve(diff(y,t) == -(1/t)*y + y^3);
y(t) = simplify(y)