login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A295519 a(n) = e^3 * Sum_{k=0..n-1} Gamma(k + 1, 3). 0
0, 1, 5, 22, 100, 493, 2701, 16678, 116704, 923473, 8204077, 81069166, 882762292, 10503611245, 135576241957, 1886597854894, 28151936397856, 448397396131969, 7592570340752053, 136187683731334054, 2579494839314653540, 51445637954467827661 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Table of n, a(n) for n=0..21.

FORMULA

a(n) = (3*n-6)*a(n-3)+(4-4*n)*a(n-2)+(3+n)*a(n-1) for n >= 3.

E.g.f.: exp(x+2)*(Ei(1,2-2*x)-Ei(1,2)). - Robert Israel, Dec 17 2017

MAPLE

a := proc(n) option remember; if n=0 then 0 elif n=1 then 1 elif n=2 then 5 else

(3*n-6)*a(n-3)+(4-4*n)*a(n-2)+(3+n)*a(n-1) fi end: seq(a(n), n=0..21);

MATHEMATICA

a[n_] := E^3 Sum[Gamma[k + 1, 3], {k, 0, n - 1}]; Table[a[n], {n, 0, 21}]

PROG

(PARI) a=vector(1000); a[1]=1; a[2]=5; a[3]=22; for(n=4, #a, a[n] = (n+3)*a[n-1]+(4-4*n)*a[n-2]+(3*n-6)*a[n-3]); va=concat(0, vector(1000, n, a[n])) \\ Altug Alkan, Dec 17 2017

CROSSREFS

Cf A053486, A295518.

Sequence in context: A123347 A087439 A033452 * A179602 A262440 A296044

Adjacent sequences:  A295516 A295517 A295518 * A295520 A295521 A295522

KEYWORD

nonn

AUTHOR

Peter Luschny, Dec 17 2017

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 22 12:33 EST 2020. Contains 332136 sequences. (Running on oeis4.)