OFFSET
0,4
COMMENTS
From Johannes W. Meijer, Oct 20 2009: (Start)
The asymptotic expansion of the higher order exponential integral E(x,m=3,n=8) ~ exp(-x)/x^3*(1 - 27/x + 539/x^2 - 9850/x^3 + 176554/x^4 + ...) leads to the sequence given above. See A163931 and A163932 for more information.
(End)
REFERENCES
Mitrinovic, D. S. and Mitrinovic, R. S. see reference given for triangle A051379.
FORMULA
a(n) = A051379(n, 2)*(-1)^n; e.g.f.: ((log(1-x))^2)/(2*(1-x)^8).
If we define f(n,i,a)=sum(binomial(n,k)*stirling1(n-k,i)*product(-a-j,j=0..k-1),k=0..n-i), then a(n) = |f(n,2,8)|, for n>=1. - Milan Janjic, Dec 21 2008
MATHEMATICA
With[{nn=20}, CoefficientList[Series[(Log[1-x])^2/(2(1-x)^8), {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Jul 10 2013 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
STATUS
approved