login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A306495 Expansion of e.g.f. (2-exp(-x))*exp(x)/(x-1)^2. 2

%I #9 Feb 19 2019 18:50:53

%S 1,4,16,74,402,2542,18446,151482,1390738,14126582,157365222,

%T 1908110866,25022451482,352918443438,5327630246542,85716034274282,

%U 1464281837606946,26470821156031462,504879319309407158,10132393298394712002,213441590598213760042

%N Expansion of e.g.f. (2-exp(-x))*exp(x)/(x-1)^2.

%H Alois P. Heinz, <a href="/A306495/b306495.txt">Table of n, a(n) for n = 0..448</a>

%F a(n) = Sum_{k=-n..n} A324224(n+1,k).

%F a(n) = (2*n+1)*a(n-1) - (n+2)*(n-1)*a(n-2) + (n-1)*(n-2)*a(n-3) for n > 2, a(n) = 4^n for n < 3.

%p egf:= (2-exp(-x))*exp(x)/(x-1)^2:

%p a:= n-> n! * coeff(series(egf, x, n+1), x, n):

%p seq(a(n), n=0..23);

%p # second Maple program:

%p a:= proc(n) option remember; `if`(n<3, 4^n,

%p (2*n+1)*a(n-1)-(n+2)*(n-1)*a(n-2)+(n-1)*(n-2)*a(n-3))

%p end:

%p seq(a(n), n=0..23);

%Y Row sums of A324224(n+1).

%K nonn

%O 0,2

%A _Alois P. Heinz_, Feb 19 2019

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 00:26 EDT 2024. Contains 371264 sequences. (Running on oeis4.)