OFFSET
0,1
COMMENTS
The higher order exponential integrals, see A163931, are defined by E(x,m,n) = x^(n-1)*int(E(t,m-1,n)/t^n, t=x..infinity) for m>=1 and n>=1, with E(x,m=0,n) = exp(-x).
The series expansions of the higher order exponential integrals are dominated by the gamma(k,n) and the alpha(k,n) constants, see A163927.
The values of the gamma(k,n) = G(k,n) coefficients can be determined with the Maple program.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..5000
J. W. Meijer and N. H. G. Baken, The Exponential Integral Distribution, Statistics and Probability Letters, Volume 5, No.3, April 1987. pp 209-211.
FORMULA
G(2,1) = gamma(2,1) = gamma^2/2+Pi^2/12.
G(k,n) = (1/k)*(gamma*G(k-1,n)) - (1/k)*Sum_{p=1..n-1}(p^(-1))* G(k-1,n) + (1/k)* Sum_{i=0..k-2}(Zeta(k-i) * G(i,n)) - (1/k)*Sum_{i=0..k-2}(Sum_{p=1..n-1}(p^(i-k)) * G(i,n)) with G(0,n) = 1 for k>=0 and n>=1.
G(k,n+1) = G(k,n) -G(k-1,n)/n.
GF(z,n) = GAMMA(n-z)/GAMMA(n).
EXAMPLE
G(2,1) = 0.9890559953279725553953956515...
MAPLE
ncol:=1; nmax:=5; kmax:=nmax; for n from 1 to nmax do G(0, n):=1 od: for n from 1 to nmax do for k from 1 to kmax do G(k, n):= expand((1/k)*((gamma-sum(p^(-1), p=1..n-1))* G(k-1, n)+sum((Zeta(k-i)-sum(p^(-(k-i)), p=1..n-1))*G(i, n), i=0..k-2))) od; od: for k from 0 to kmax do G(k, ncol):=G(k, ncol) od;
MATHEMATICA
RealDigits[ N[ EulerGamma^2/2 + Pi^2/12, 105]][[1]] (* Jean-François Alcover, Nov 07 2012, from 1st formula *)
CROSSREFS
KEYWORD
dead
AUTHOR
Johannes W. Meijer and Nico Baken, Aug 13 2009, Aug 17 2009
STATUS
approved