10   point 10
   20   dim Bern(100)
   30   Bern(2)=1//6
   40   Bern(4)=1//30
   50   Bern(6)=1//42
   60   Bern(8)=1//30
   70   Bern(10)=5//66
   80   Bern(12)=691//2730
   90   Bern(14)=7//6
  100   Bern(16)=3617//510
  110   Bern(18)=43867//798
  120   Bern(20)=174611//330
  130   Bern(22)=854513//138
  140   Bern(24)=236364091//2730
  150   Bern(26)=8553103//6
  160   Bern(28)=23749461029//870
  170   Bern(30)=8615841276005//14322
  180   Bern(32)=7709321041217//510
  190   Bern(34)=2577687858367//6
  640   A=10:J=34:X#=0
  650   for S=2 to 36 step 2
  660     X#=0:Y#=0
  670     for N=1 to A-1
  680       X#=X#+1/(N^S)
  690     next N
  700     Y#=1/(S-1)/(A^(S-1))+1/2/(A^S)
  710     for N=2 to J step 2
  720       Y#=Y#+(-1)^(N\2+1)*Bern(N)*!(S+N-2)/!(S-1)/!(N)/(A^(S+N-1))
  750     next N
  760     D=15
  770     U#=int((X#+Y#)*10^D)/10^D:V#=int((X#+Y#)*(1-1/2^S)*10^D)/10^D
  780     Z#=int((X#+Y#)*(1-1/3^S)*10^D)/10^D:W#=int((X#+Y#)*(1-1/6^S)*10^D)/10^D
  790     A#=int((1-1/2^S)*10^D)/10^D
  800     B#=int((1-1/3^S)*10^D)/10^D:C#=int((1-1/6^S)*10^D)/10^D
  802     D#=int((1/6^S)*10^D)/10^D
  803     E#=int((1/25^S+1/35^S+1/49^S+1/55^S+1/65^S+1/77^S+1/91^S+1/121^S+1/163^S+1/169^S+1/289^S)*10^(D+5))/10^(D+5)
  810     print S,V#+Z#-W#-A#-B#+C#+D#-E#
  820   next S
  860   end