integrate 1/(8*pi**2)*ln(21-4*(cos(x)+cos(y)+cos(x+y)))/ln(9/2) for x = 0..2pi and y = 0.. 2pi