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