var('x') y = function('y', x) f(x) = desolve(diff(y, x) + y/x, y) show(f(x))

plot(f(x).subs(c=1), (x, 0.1, 1), detect_poles=True)

phanindra,student,nmrec