f033b4fd0c1d4d751e08c65f60e7fe31  a36_c19a_check.py
X case2_identity: True
X case3_crit: True
X case3_at_3/2: True
X case3_decreasing: True
X case4_deriv: True
X case4_at_5/2: True
X c_gt_32/57: True
X P1: True
X P2: True
X case3_formula: True
X arb: c = [0.606530659712633 +/- 4.26e-16]  c - 32/57 = [0.045127150940704 +/- 7.04e-16]
PART X PASS
D closed form vs definition, 2,000 points of [-12, 4]: max diff 1.7106e-49
D max of D(-g)/g^2 + c/2 over 300000 grid points of (0, 30]: 0.0 at g = 1.0
D values at g = 0.999, 1.001: ['-3.03873e-7', '-3.0266e-7']
PART D PASS
K Richardson derivative vs D(gamma)/(sqrt2 gamma^2) at 24 points: max diff 2.293e-19
K largest derivative over the negative points: -0.214440971240177 at x = -1.4142136; -sqrt(2/e)/4 = -0.214440971240177
PART K PASS

real	1m28.826s
user	1m28.770s
sys	0m0.032s
