python3 -c 'from decimal import * getcontext().prec=80 D=Decimal s=D(3).sqrt();b=D(2).sqrt() d11=-4-4*b-s;d13=-4-4*b;d15=-4-4*b+s;d31=-4-2*b-3*s-2*b*s;d35=-4-2*b-s-2*b*s;d41=-12-s-4*b*s;d43=-12-4*b*s;d45=-12+s-4*b*s S1=[16*s*((3*s-5)/(24*d11)+(1-s)/(12*d13)+(s+1)/(48*d15)),16*s*(-(s-1)**2/(48*d11)+1/(3*d13)-(s+1)**2/(48*d15)),16*s*(-(s-1)/(48*d11)+(1+s)/(12*d13)-(3*s+5)/(24*d15))] S2=[(2*s-3)/4,(s-3)/2,-s/4] S3=[s*(3*s-5)/(3*d31),(2*s/3)*((2-s)/d31+(2+s)/d35),-s*(5+3*s)/(3*d35)] S4=[(s-2)/12,-(s-1)/6,D(1)/12] S5=[(2*s/3)*(5-3*s)/d41+(4*s/3)*(s-1)/d43+(s/3)*(1+s)/d45,-(s/3)*(s-1)**2/d41+(16*s/3)/d43-(s/3)*(s+1)**2/d45,(s/3)*(1-s)/d41-(4*s/3)*(s+1)/d43+(2*s/3)*(5+3*s)/d45] dir=[-(3-s)/192,-s/48,(3+s)/192] nu=2+b+2*s G=[dir[i]+(3*S1[i]+6*S2[i]+12*S3[i]+6*S4[i]+3*S5[i])/192 for i in range(3)] for i,den in enumerate([2-s,D(2),2+s]): g=G[i]/(den-nu**2) print(g);print(g*D(100000))'