DSolve[Derivative[1][y][x] == (4*a*x - y[x]^2)^3/(y[x]*(-1 + 4*a*x - y[x]^2)),y[x],x]