a(n+1) is the moment of order n for the probability density function rho(x) = Pi^(-3/2)*sqrt(x/2)*exp(x/2)/(1-erf^2(I*sqrt(x/2))) on the interval 0..infinity, with erf the error function and I=sqrt(-1). [Groux Roland, Nov 10 2009]