DSolve[y[x]*(c + b*x^3*y[x]) + x*(c + a*x*y[x]^3)*Derivative[1][y][x] == 0,y[x],x]