OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: 1/sqrt(1-8*x^2)+(1-sqrt(1-8*x^2))/(4*x*sqrt(1-8*x^2)) = (1+4*x-sqrt(1-8*x^2))/(4*x*sqrt(1-8*x^2)).
a(n) = binomial(n, floor(n/2))2^floor(n/2).
a(n+1) = (1/Pi)*int(x^n*(x+4)/sqrt(8-x^2),x,-2*sqrt(2),2*sqrt(2)) if n is odd [corrected by Vaclav Kotesovec, Nov 13 2017].
Conjecture: (n+1)*a(n) +(n-1)*a(n-1) -n*a(n-2) +(2-n)*a(n-3) = 0. - R. J. Mathar, Nov 15 2011
MATHEMATICA
nmax = 30; CoefficientList[Series[BesselI[0, 2*Sqrt[2]*x] + BesselI[1, 2*Sqrt[2]*x]/Sqrt[2], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Nov 13 2017 *)
PROG
(PARI) x='x+O('x^30); Vec((1+4*x-sqrt(1-8*x^2))/(4*x*sqrt(1-8*x^2))) \\ G. C. Greubel, Aug 17 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); Coefficients(R!((1+4*x-Sqrt(1-8*x^2))/(4*x*Sqrt(1-8*x^2)))); // G. C. Greubel, Aug 17 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paul Barry, Sep 20 2004
STATUS
approved