n=1: max over p in (0,q] of W1 - (qp + beta p^2) = -2.776e-17
n=2: max over p in (0,q] of W1 - (qp + beta p^2) = +2.776e-17
n=3: max over p in (0,q] of W1 - (qp + beta p^2) = +5.551e-17
n=4: max over p in (0,q] of W1 - (qp + beta p^2) = +1.110e-16
n=5: max over p in (0,q] of W1 - (qp + beta p^2) = +1.665e-16
overall max excess +1.665e-16 at n=5, p=0.4000
