python3 -c 'import math
def b(j,m): return math.sqrt(1/3)*(math.sin((j+1)*m*math.pi/6)-math.sin(j*m*math.pi/6))
def C(r,p,q,s): return sum(b(j,r)*b(j,p)*b(j,q)*b(j,s) for j in range(6))
print(C(1,2,3,4))'