python3 -c 'import math s=math.sqrt(3);b=math.sqrt(2) def f(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,t):return sum(f(j,r)*f(j,p)*f(j,q)*f(j,t) for j in range(6)) print(C(1,1,2,2),C(1,5,2,2),C(1,3,4,4),C(5,3,4,4),C(2,3,4,4),C(4,2,2,4))'