(1) n=3: N_3 == -512 (Q1-Q2): True
    M_1 == 16*((d_2-d_3)^2 + (d_2+d_3)/2 - 3/16): True
    M_2 == 16*((d_1-d_3)^2 + (d_1+d_3)/2 - 3/16): True
    M_3 == 16*((d_1-d_2)^2 + (d_1+d_2)/2 - 3/16): True
(2) n=4: N_4 has 424 monomials (degree 8), M_i has 35 monomials (degree 4); written to N4_poly.txt
    random rational points: 1499, agreement 1499; counts {'first': 667, 'second': 540, 'outside': 292}
(3) n=5: 3000 random points, agreement 3000
(3) n=6: 3000 random points, agreement 3000
(2b) n=4: Q1-Q2 proportional to N_4? False ; Q1-Q2 at (1/4,1/4,1/4,1/4) = 1/16  vs N_4 there = 65536
(4) r_j = 9/20 for all j (d_j = 319/1600): polygon 2r <= 2 holds (interior of G(B)); N_4 = 5.036818e+03
(4) r_j = 11/20 for all j (d_j = 279/1600): polygon 2r <= 2 holds (interior of G(B)); N_4 = -3.257437e+03
DONE
