q=1/5 n=1: max_f sum Delta^2 = 1 = max_k k(1-q)^(2k-2)
q=1/5 n=2: max_f sum Delta^2 = 32/25 = max_k k(1-q)^(2k-2)
q=1/5 n=3: max_f sum Delta^2 = 32/25 = max_k k(1-q)^(2k-2)
q=1/3 n=1: max_f sum Delta^2 = 1 = max_k k(1-q)^(2k-2)
q=1/3 n=2: max_f sum Delta^2 = 1 = max_k k(1-q)^(2k-2)
q=1/3 n=3: max_f sum Delta^2 = 1 = max_k k(1-q)^(2k-2)
q=2/7 n=1: max_f sum Delta^2 = 1 = max_k k(1-q)^(2k-2)
q=2/7 n=2: max_f sum Delta^2 = 50/49 = max_k k(1-q)^(2k-2)
q=2/7 n=3: max_f sum Delta^2 = 50/49 = max_k k(1-q)^(2k-2)
q=9/20 n=1: max_f sum Delta^2 = 1 = max_k k(1-q)^(2k-2)
q=9/20 n=2: max_f sum Delta^2 = 1 = max_k k(1-q)^(2k-2)
q=9/20 n=3: max_f sum Delta^2 = 1 = max_k k(1-q)^(2k-2)
C2: identity exact on 3000 random cases (11257 endpoint coincidences exercised)
C3: {L<0} = {no success}, E|Y-mu| = 2 q (sum alpha)(1-q)^n and the bound, exact on 679 cases
C4: martingale coupling X = E[N|X] exists (e^-q >= 1-q; residual law on {X=1} has mass 1, mean 1)
C4: E|N_1/2 - 1/2| = e^-1/2 exactly, so R_Pois >= sqrt(2/e) = 0.857763884961
C4: convex-order inequality E|sum a(X-q)| <= E|sum a(N-q)| on 300 random cases (numerical evidence): max (lhs-rhs)/lhs = -5.03e-05
ALL C1-C4 CHECKS PASSED

real	0m5.901s
user	0m5.863s
sys	0m0.036s
