• Book	Advanced Engineering Mathematics
  Author	Erwin Kreyszig
  Edition	8th Edition

code_begins

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)

code_ends

• Solution by:

phanindra,student,nmrec