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