(PARI){数字我的数字=20000; 默认值(realprecision,数字我的数字+80); s=1.0;n=5*数字我的数字; j=4*n+1;si=-1.0;对于(i=3,j-2,s+=si/i^2;si=-si;i++;);s+=0.5/j^2;ttk=4.0;d=4.0*j^3;xk=2.0;xkp=xk;对于(k=2100000000,项=(ttk-1)*ttk*xkp;xk++;xkp*=xk;如果(k>2,项*=xk;xk++;xkp*=xk;);术语*=bernreal(k)/d;sn=s+项;如果(sn==s,断裂);s=sn;ttk*=4.0;d*=(k+1)*(k+2)*j^2;k++;);x=10*s;对于(n=0,数字, 我的数字, d=地板(x);x=(x-d)*10;写(“b006752.txt”,n,“”,d);}/*贝塔系数(2)=1-1/3^2+1/5^2-…-1/(J-2)^2+1/(2*J^2)+2*Bern(0)/(2*J^3)-2*3*4*Bern(2)/J^5+, */