10 point 10 20 dim Bern(50) 30 dim Zr(50) 40 dim Br(50) 50 Bern(2)=1//6 60 Bern(4)=1//30 70 Bern(6)=1//42 80 Bern(8)=1//30 90 Bern(10)=5//66 100 Bern(12)=691//2730 110 Bern(14)=7//6 120 Bern(16)=3617//510 130 Bern(18)=43867//798 140 Bern(20)=174611//330 150 Bern(22)=854513//138 160 Bern(24)=236364091//2730 170 Bern(26)=8553103//6 180 Bern(28)=23749461029//870 190 Bern(30)=8615841276005//14322 200 Bern(32)=7709321041217//510 210 Bern(34)=2577687858367//6 220 Bern(36)=26315271553053477373//1919190 230 A=10:J=34:X#=0 240 for S=2 to 36 step 2 250 X#=0:Y#=0 260 for N=1 to A-1 270 X#=X#+1/(N^S) 280 next N 290 Y#=1/(S-1)/(A^(S-1))+1/2/(A^S) 300 for N=2 to J step 2 310 Y#=Y#+(-1)^(N\2+1)*Bern(N)*!(S+N-2)/!(S-1)/!(N)/(A^(S+N-1)) 320 next N 330 D=20 340 Zr(S)=int(log(X#+Y#)*10^D) 350 next S 360 Br(36)=14551:Br(34)=58207:Br(32)=232830:Br(30)=931326 '326 370 Br(28)=3725333:Br(26)=14901555:Br(24)=59608184:Br(22)=238450446 380 Br(20)=953961124:Br(18)=3817278702:Br(16)=15282026219 '1124 390 Br(14)=61244396725:Br(12)=246026470035 400 Br(10)=993603574437:Br(8)=4061405366518 410 Br(6)=17070086850639:Br(4)=76993139764252 420 for N=36 to 2 step -2:Br(N)=Br(N)*10^(D-15):next N 430 S0=16 440 for S=S0 to 2 step -2 450 K=2:M=K*S:W=0 460 repeat 470 if M<37 then W=W+Br(M)/K 480 if M>36 then W=W+fnPzeta(M)/K 490 K=K+1:M=K*S 500 until M>100 510 Br(S)=int((Zr(S)-W)/10^(D-15))*10^(D-15) 520 next S 530 for S=2 to S0 step 2:print S,Br(S)/10^D:next S 540 end 550 fnPzeta(S) 560 local X 570 X=1/2^S+1/3^S+1/5^S+1/7^S+1/11^S+1/13^S+1/17^S 580 return(X*10^D)