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)