OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..29
FORMULA
a(n) ~ Pi * exp(n^2/2 - n + 1/12) * n^(n^2 + n + 17/12) / (A * 2^(2*n^2 - 25/12)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Oct 15 2019
MAPLE
seq(`if`(n=0, 0, 2*factorial(n^2)*mul(factorial(k)/factorial(n+k), k = 1 .. n-1)), n = 0..10); # G. C. Greubel, Oct 12 2019
MATHEMATICA
Table[If[n==0, 0, 2*(n^2)!*Product[k!/(n+k)!, {k, 1, n-1}]], {n, 0, 10}] (* G. C. Greubel, Oct 12 2019 *)
Flatten[{0, Table[2*BarnesG[n+1]*BarnesG[n+2]*(n^2)!/BarnesG[2*n + 1], {n, 1, 10}]}] (* Vaclav Kotesovec, Oct 15 2019 *)
PROG
(PARI) concat([0], vector(10, n, 2*(n^2)!*prod(k=1, n-1, k!/(n+k)!) )) \\ G. C. Greubel, Oct 12 2019
(Magma) [0, 2] cat [2*Factorial(n^2)*(&*[Factorial(k)/Factorial(n+k): k in [1..n-1]]): n in [2..10]]; // G. C. Greubel, Oct 12 2019
(Sage) [0]+[2*factorial(n^2)*product(factorial(k)/factorial(n+k) for k in (1..n-1)) for n in (1..10)] # G. C. Greubel, Oct 12 2019
(GAP) Concatenation([0], List([1..10], n-> 2*Factorial(n^2)*Product([1..n-1], k-> Factorial(k)/Factorial(n+k)) )); # G. C. Greubel, Oct 12 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Feb 13 2002
EXTENSIONS
Definition corrected by Neven Juric, Nov 08 2012
STATUS
approved