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