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