numerator of d gamma / dY: -M**2*(4*M**3 + 6*M**2*Y - 3*M**2 - 6*M*Y - 2*M - 2*Y**3 - 3*Y**2 + 1)
symbolic checks: LQ, T, (T) identities, Dg+, gamma closed form, gamma(M/2), derivative numerator 216M^2-96M+96,
  (g3) parametrisation, max of s/(c+s)^3 = 4/(27c^2) at s = c/2, constants 3840/27636 and 1024/25947,
  first integral and equilibrium curve of (5) (all equilibria on q = 0), variant: first integral and
  equilibrium curve a = y, b = x, xy = lambda(x+y)
TOTAL FAILS 0
nice -n 10 python3 -I r2sym.py   [a Python 3 with sympy; interpreter path omitted]  0.43s user 0.03s system 99% cpu 0.463 total
