OFFSET
0,2
COMMENTS
Binomial transform of A081672.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..500
FORMULA
E.g.f. exp(3*x) - exp(x)*(1-BesselI_0(2*x)).
Conjecture: n*(2*n - 7)*a(n) +(-12*n^2 + 50*n - 33)*a(n-1) +(16*n^2 - 76*n + 87)*a(n-2) +3*(4*n^2 - 22*n + 27)*a(n-3) -9*(2*n-5)*(n-3)*a(n-4)=0. - R. J. Mathar, Nov 24 2012
a(n) ~ 3^n * (1 + sqrt(3)/(2*sqrt(Pi*n))). - Vaclav Kotesovec, Jul 02 2015
From Robert Israel, Jun 03 2016: (Start)
E.g.f. A(x) satisfies
-27 A + (-63 x + 9) A' + (-18 x^2 + 42 x + 39) A'' + (12 x^2 + 68 x - 25) A''' + (16 x^2 - 58 x + 4) A'''' + (-12 x^2 + 11 x) A''''' + 2 x^2 A'''''' = 0.
This implies Mathar's conjectured recursion. (End)
MAPLE
Egf:= exp(3*x)-exp(x)*(1-BesselI(0, 2*x)):
S:= series(Egf, x, 101):
seq(coeff(S, x, n)*n!, n=0..100); # Robert Israel, Jun 03 2016
MATHEMATICA
CoefficientList[Series[E^(3*x)-E^x*(1-BesselI[0, 2*x]), {x, 0, 50}], x] * Range[0, 50]! (* Vaclav Kotesovec, Jul 02 2015 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Mar 28 2003
STATUS
approved