Re Bernoulli numbers

N. J. A. Sloane njas at research.att.com
Tue Feb 3 21:39:10 CET 2004


Thanks, Edwin!

%I A090495
%S A090495 574,1185,1240,1269,1376,1906,1910
%N A090495 Numbers n such that numerator(Bernoulli(2*n)/(2*n)) is different from numerator(Bernoulli(2*n)/(2*n*(2*n-1))).
%C A090495 n such that A001067 is different from A046968.
%O A090495 1,1
%K A090495 nonn,more
%Y A090495 Cf. A090496, A001067, A046968.
%p A090495 a:=n->numer(bernoulli(2*n)/(2*n)): b:=n->numer(bernoulli(2*n)/(2*n*(2*n-1))): for n from 1 to 2000 do if a(n)<>b(n) then print(n,a(n)/b(n));  fi; od:
%t A090495 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]] ] ]
%A A090495 njas, Feb 03 2004
%E A090495 a(1) discovered by Michael Somos, Feb 01 2004. a(2)-a(7) from Edwin Clark, Feb 03, 2004.

%I A090496
%S A090496 37,103,37,59,131,37,67
%N A090496 Ratio of numerator(Bernoulli(2*n)/(2*n)) to numerator(Bernoulli(2*n)/(2*n*(2*n-1))) for n's for which they are different.
%C A090496 A001067(n) / A046968(n) when they are different.
%O A090496 1,1
%K A090496 nonn,more
%Y A090496 Cf. A090496, A001067, A046968.
%A A090496 njas, Feb 03 2004
%E A090496 a(1) discovered by Michael Somos, Feb 01 2004. a(2)-a(7) from Edwin Clark, Feb 03, 2004.

NJAS





More information about the SeqFan mailing list