y=function('y',x)
deq=x*diff(y,x,1)-(y-x)^3-y==0
desolve(deq,[y,x],contrib_ode=True)