var('x')
y=function('y',x)
print desolve(diff(y,x)-y*(1+3*(x^(-1)))-x-2,[y,x],ics=(1,exp(1)-1))