OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..220
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem., arXiv:quant-phys/0402027, 2004.
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, Phys. Lett. A 309 (2003) 198-205.
M. Schork, On the combinatorics of normal ordering bosonic operators and deforming it, J. Phys. A 36 (2003) 4651-4665.
FORMULA
a(n) = (2*n)!*hypergeom([2*n+1], [3], 1)/(2*exp(1)), n=1, 2, ... Special values of the confluent hypergeometric function 1F1.
a(n) = sum(A090438(n, k), k=2..2*n)= sum((1/k!)*product(fallfac(k+(j-1)*(4-2), 2), j=1..n), k=1..infinity)/exp(1), n>=1. From eq.(9) of the Blasiak et al. reference with r=4, s=2. fallfac(n, m) := A008279(n, m) (falling factorials triangle). a(0) := 1 may be added.
8*n*(2*n-1)*(2*n+1)*(n+1)^2*(n+3)*(n+2)*a(n)+(2*(n+1))*(8*n^3+32*n^2+42*n+13)*a(n+1)*(n+3)*(n+2)-(8*n^2+38*n+51)*(n+3)*(n+2)*a(n+2)+(n+3)*(n+2)*a(n+3) = 0. - Robert Israel, May 23 2016
a(n) = A052852(2*n-1). - Mark van Hoeij, Sep 05 2022
MAPLE
f:= n -> simplify((2*n)!*hypergeom([2*n+1], [3], 1)/(2*exp(1))):
map(f, [$1..30]); # Robert Israel, May 23 2016
MATHEMATICA
a[n_] := n*(2n-1)!*Hypergeometric1F1[2-2n, 3, -1]; Array[a, 30] (* Jean-François Alcover, Sep 01 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Karol A. Penson, Jul 01 2002
EXTENSIONS
Edited by Wolfdieter Lang, Dec 23 2003
STATUS
approved