var('x')
y=function('y',x)
print desolve(diff(y,x)+3*y/x-(1/x),[y,x])