# D (derived independently) == R U R on 144 monomials: PASS
k_2 = 4  == direct
k_4 = 8  == direct
k_6 = 80/3  == direct
k_8 = 304/3  == direct
k_10 = 6112/15  == direct
k_12 = 25328/15  == direct
k_14 = 448864/63  == direct
k_16 = 212924/7  == direct
k_18 = 53017256/405  == direct
k_20 = 892403651/1575  == direct
# monomials in B^42 1: 3784;  time 0.6s
# all cumulants k_2..k_44 from the model compared with the direct ones (available up to k_44): ALL EQUAL
# B^1 1 = (1)*x1^0*x3^1*s^0 + (1)*x1^1*x3^0*s^0
#   <B^1 1> = 0
# B^2 1 = (1)*x1^0*x3^0*s^0 + (1)*x1^0*x3^2*s^0 + (2)*x1^1*x3^1*s^0 + (1)*x1^2*x3^0*s^0
#   <B^2 1> = 2
# B^3 1 = (1)*x1^0*x3^1*s^0 + (1)*x1^0*x3^1*s^1 + (1)*x1^0*x3^3*s^0 + (2)*x1^1*x3^0*s^0 + (-1)*x1^1*x3^0*s^1 + (3)*x1^1*x3^2*s^0 + (3)*x1^2*x3^1*s^0 + (1)*x1^3*x3^0*s^0
#   <B^3 1> = 0
#   <B^4 1> = 20/3
