A241885
Write the coefficient of x^n/n! in the expansion of (x/(exp(x)-1))^(1/2) as f(n)/g(n); sequence gives f(n).
Original entry on oeis.org
1, -1, 1, 1, -3, -19, 79, 275, -2339, -11813, 14217, 95265, -4634445, -193814931, 131301607, 1315505395, -3890947599, -136146236611, 46949081169401, 124889801445461, -10635113572583999, -158812278992229461, 56918172351554857, 8484151253958927197
Offset: 0
For n=1, B_1=-1/2 and B_1^(1/2)=-1/4 so a(1)=-1.
For n=6, B_6=1/6 and B_6^(1/2)=79/86016 so a(6)=79.
1/1, -1/4, 1/48, 1/64, -3/1280, -19/3072, 79/86016, 275/49152, -2339/2949120, -11813/1310720, 14217/11534336 = A241885 / A242225.
-
g := proc(f, n) option remember; local g0, m; g0 := sqrt(f(0));
if n=0 then g0 else if n=1 then 0 else add(binomial(n, m)*g(f,m)*g(f,n-m), m=1..n-1) fi; (f(n)-%)/(2*g0) fi end:
a := n -> numer(g(bernoulli, n));
seq(a(n), n = 0..23); # Peter Luschny, May 07 2014
-
a := 1
g[0] := Sqrt[f[0]]
f[k_] := BernoulliB[k]
g[1] := f[1]/(2 g[0]^1);
g[k_] := (f[k] - Sum[Binomial[k, m] g[m] g[k - m], {m, 1, k - 1}])/(2 g[0])
Table[Factor[g[k]], {k, 0, 15}] // TableForm
(* Alternative: *)
Table[Numerator@NorlundB[n, 1/2, 0], {n, 0, 23}] (* Peter Luschny, Feb 18 2024 *)
-
a(n)=numerator(sum(k=0,n,binomial(-1/2,k)*binomial(n+1/2,n-k)*stirling(n+k,k,2)/binomial(n+k,k))) \\ Tani Akinari, Oct 08 2024
A222411
Numerators in Taylor series expansion of (x/(exp(x) - 1))^(3/2)*exp(x/2).
Original entry on oeis.org
1, -1, -1, 5, 7, -19, -869, 715, 2339, -200821, -12863, 2117, 7106149, -64604977, -131301607, 7629931291, 174053933, -19449462373, -46949081169401, 355455588729389, 10635113572583999, -6511303438681407901, -349640201588122693, 9112944418860287
Offset: 0
The first few fractions are 1, -1/4, -1/32, 5/384, 7/10240, -19/40960, -869/61931520, 715/49545216, ... = A222411/A222412. - _Petros Hadjicostas_, May 14 2020
- Alois P. Heinz, Table of n, a(n) for n = 0..300
- David Broadhurst, Relations between A241885/A242225, A222411/A222412, and A350194/A350154.
- F. J. Dyson, N. E. Frankel and M. L. Glasser, Lehmer's Interesting Series, arXiv:1009.4274 [math-ph], 2010-2011.
- F. J. Dyson, N. E. Frankel and M. L. Glasser, Lehmer's interesting series, Amer. Math. Monthly, 120 (2013), 116-130.
- D. H. Lehmer, Interesting series involving the central binomial coefficient, Amer. Math. Monthly, 92(7) (1985), 449-457.
-
gf:= (x/(exp(x)-1))^(3/2)*exp(x/2):
a:= n-> numer(coeff(series(gf, x, n+3), x, n)):
seq(a(n), n=0..25); # Alois P. Heinz, Mar 02 2013
-
Series[(x/(Exp[x]-1))^(3/2)*Exp[x/2], {x, 0, 25}] // CoefficientList[#, x]& // Numerator (* Jean-François Alcover, Mar 18 2014 *)
A222412
Denominators in Taylor series expansion of (x/(exp(x) - 1))^(3/2)*exp(x/2).
Original entry on oeis.org
1, 4, 32, 384, 10240, 40960, 61931520, 49545216, 7927234560, 475634073600, 1993133260800, 177167400960, 48753634065776640, 195014536263106560, 39002907252621312000, 842462796656620339200, 2204424056667635712000, 79359266040034885632000
Offset: 0
The first few fractions are 1, -1/4, -1/32, 5/384, 7/10240, -19/40960, -869/61931520, 715/49545216, ... = A222411/A222412. - _Petros Hadjicostas_, May 14 2020
- Alois P. Heinz, Table of n, a(n) for n = 0..300
- David Broadhurst, Relations between A241885/A242225, A222411/A222412, and A350194/A350154.
- F. J. Dyson, N. E. Frankel and M. L. Glasser, Lehmer's Interesting Series, arXiv:1009.4274 [math-ph], 2010-2011.
- F. J. Dyson, N. E. Frankel and M. L. Glasser, Lehmer's interesting series, Amer. Math. Monthly, 120 (2013), 116-130.
- D. H. Lehmer, Interesting series involving the central binomial coefficient, Amer. Math. Monthly, 92(7) (1985), 449-457.
-
gf:= (x/(exp(x)-1))^(3/2)*exp(x/2):
a:= n-> denom(coeff(series(gf, x, n+3), x, n)):
seq(a(n), n=0..25); # Alois P. Heinz, Mar 02 2013
-
Series[(x/(Exp[x]-1))^(3/2)*Exp[x/2], {x, 0, 25}] // CoefficientList[#, x]& // Denominator (* Jean-François Alcover, Mar 18 2014 *)
A242225
Write the coefficient of x^n/n! in the expansion of (x/(exp(x)-1))^(1/2) as f(n)/g(n); sequence gives g(n).
Original entry on oeis.org
1, 4, 48, 64, 1280, 3072, 86016, 49152, 2949120, 1310720, 11534336, 4194304, 1526726656, 2348810240, 12079595520, 3221225472, 73014444032, 51539607552, 137095356088320, 5772436045824, 3809807790243840, 725677674332160, 2023101395107840, 3166593487994880
Offset: 0
For n=1, B_1=-1/2 and B_1^(1/2)=-1/4 so a(1)=4.
For n=6, B_6=1/6 and B_6^(1/2)=79/86016 so a(6)=86016.
-
g := proc(f, n) option remember; local g0, m; g0 := sqrt(f(0));
if n=0 then g0 else if n=1 then 0 else add(binomial(n, m)*g(f, m)*g(f, n-m), m=1..n-1) fi; (f(n)-%)/(2*g0) fi end:
a := n -> denom(g(bernoulli, n));
seq(a(n), n=0..23);
-
a := 1
g[0] := Sqrt[f[0]]
f[k_] := BernoulliB[k]
g[1] := f[1]/(2 g[0]^1);
g[k_] := (f[k] - Sum[Binomial[k, m] g[m] g[k - m], {m, 1, k - 1}])/(2 g[0])
Table[Denominator[Factor[g[k]]], {k, 0, 15}] // TableForm
(* Alternative: *)
Table[Denominator@NorlundB[n, 1/2, 0], {n, 0, 23}] (* Peter Luschny, Feb 18 2024 *)
A350194
Numerators of a power series characterizing how powers of the cosine function converge to the Gaussian function.
Original entry on oeis.org
1, -1, 1, 1, -1, -19, 79, 11, -2339, -11813, 677, 2117, -308963, -64604977, 131301607, 263101079, -5614643, -1768132943, 46949081169401, 9606907803497, -10635113572583999, -158812278992229461, 8131167478793551, 9112944418860287, -40395223967437706149
Offset: 0
Showing 1-5 of 5 results.
Comments