Bernoulli numbers

[N. J. A. Sloane]

Dear SeqFans,   Michael Somos made an interesting discovery.
A001067 and A046968 are not the same!
In fact: (speaking baby Mathematica talk, which is
all that i know): 

a[n_] := Numerator[BernoulliB[2n]/(2n)]          (* A001067 *)
b[n_] := Numerator[BernoulliB[2n]/(2n(2n-1))]    (* A046968 *)

For[n=1, n <= 580, n++, 
If[ a[n] != b[n], Print[n, " ", a[n]/b[n]] ]
]

produces one line of output

                  574 37

In other words, the sequence of values of n such that
A001067(n) differs from A046968(n) starts
574, ...
and the associated ratios begin
37, ...

Could someone with more computing power than I have extend
these two sequences?  PARI gives the same result.  Maple dies.


Thanks

Neil


 Neil J. A. Sloane
 AT&T Shannon Labs, Room C233, 
 180 Park Avenue, Florham Park, NJ 07932-0971
 Email: njas at research.att.com
 Office: 973 360 8415; fax: 973 360 8178
 Home page: http://www.research.att.com/~njas/