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