OFFSET
0,2
COMMENTS
The asymptotic expansion of the higher order exponential integral E(x,m=2,n=5) ~ exp(-x)/x^2*(1 - 11/x + 107/x^2 - 1066/x^3 + 11274/x^4 - 127860/x^5 + 1557660/x^6 - ... ) leads to the sequence given above. See A163931 and A028421 for more information. - Johannes W. Meijer, Oct 20 2009
REFERENCES
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
T. D. Noe, Table of n, a(n) for n = 0..100
D. S. Mitrinovic and R. S. Mitrinovic, Tableaux d'une classe de nombres reliƩs aux nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz. No. 77 1962, 77 pp.
FORMULA
a(n) = Sum_{k=0..n} (-1)^(n+k)*binomial(k+1, 1)*5^k*Stirling1(n+1, k+1). - Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004
a(n) = n!*Sum_{k=0..n-1} (-1)^k*binomial(-5,k)/(n-k). - Milan Janjic, Dec 14 2008
a(n) = n!*[4]h(n), where [k]h(n) denotes the k-th successive summation of the harmonic numbers from 0 to n. With offset 1. - Gary Detlefs, Jan 04 2011
E.g.f.: (1 + 5*log(1/(1-x)))/(1 - x)^6. - Ilya Gutkovskiy, Jan 23 2017
MATHEMATICA
f[k_] := k + 4; t[n_] := Table[f[k], {k, 1, n}]; a[n_] := SymmetricPolynomial[n - 1, t[n]]; Table[a[n], {n, 1, 16}] (* Clark Kimberling, Dec 29 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004
STATUS
approved