OFFSET
0,3
COMMENTS
The asymptotic expansion of the higher order exponential integral E(x,m=2,n=8) ~ exp(-x)/x^2*(1 - 17/x + 242/x^2 - 3382/x^3 + 48504/x^4 - 725592/x^5 + 11393808/x^6 - ...) leads to the sequence given above. See A163931 and A028421 for more information. - Johannes W. Meijer, Oct 20 2009
REFERENCES
Mitrinovic, D. S. and Mitrinovic, R. S. see reference given for triangle A051379.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..440
FORMULA
a(n) = A051379(n, 2)*(-1)^(n-1).
E.g.f.: -log(1-x)/(1-x)^8.
a(n) = n!*Sum_{k=0..n-1} ((-1)^k*binomial(-8,k)/(n-k)), for n>=1. - Milan Janjic, Dec 14 2008
a(n) = n!*[7]h(n), where [k]h(n) denotes the k-th successive summation of the harmonic numbers from 0 to n. - Gary Detlefs, Jan 04 2011
Conjecture: a(n) +(-2*n-13)*a(n-1) +(n+6)^2*a(n-2)=0. - R. J. Mathar, Aug 04 2013
MATHEMATICA
f[k_] := k + 7; 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
easy,nonn
AUTHOR
STATUS
approved