dsolve(x*y(x)*diff(diff(y(x),x),x)-x*diff(y(x),x)^2+a*y(x)*diff(y(x),x)+b*x*y(x)^3=0,y(x))