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=36:X#=0 240 for S=2 to 48 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)) 330 next N 340 print S,int((X#+Y#)*(1-1/2^S)*10^23)/10^23 350 'print S,(int((X#+Y#)*10^23)-int((X#+Y#)*(1-1/2^S)*10^23))/10^23 360 next S 370 end