10 'E101 list PA: push F10 to continue 20 print "PA" '**PA** 30 point 10 40 dim Bern(50) 50 dim Zr(50) 60 dim Br(50) 70 Bern(2)=1//6 80 Bern(4)=1//30 90 Bern(6)=1//42 100 Bern(8)=1//30 110 Bern(10)=5//66 120 Bern(12)=691//2730 130 Bern(14)=7//6 140 Bern(16)=3617//510 150 Bern(18)=43867//798 160 Bern(20)=174611//330 170 Bern(22)=854513//138 180 Bern(24)=236364091//2730 190 Bern(26)=8553103//6 200 Bern(28)=23749461029//870 210 Bern(30)=8615841276005//14322 220 Bern(32)=7709321041217//510 230 Bern(34)=2577687858367//6 240 Bern(36)=26315271553053477373//1919190 250 A=10:J=36:X#=0 260 for S=2 to 44 step 2 270 X#=0:Y#=0 280 for N=1 to A-1 290 X#=X#+1/(N^S) 300 next N 310 Y#=1/(S-1)/(A^(S-1))+1/2/(A^S) 320 for N=2 to J step 2 330 Y#=Y#+(-1)^(N\2+1)*Bern(N)*!(S+N-2)/!(S-1)/!(N)/(A^(S+N-1)) 340 next N 350 print S,int((X#+Y#)*(1-1/2^S)*10^23)/10^23 360 next S 370 end