OFFSET
0,2
COMMENTS
Hankel transform is A006720(n+3).
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: (sqrt((1-2x)*(1-2x-4x^2+4x^3))-2x^2+4x-1)/(2x*(1-4x+3x^2)).
Conjecture: (n+1)*a(n) +2*(-4*n-1)*a(n-1) +(19*n-5)*a(n-2) -36*a(n-3) +8*(-7*n+26)*a(n-4) +2*(34*n-143)*a(n-5) +24*(-n+5)*a(n-6)=0. - R. J. Mathar, Oct 10 2014
MATHEMATICA
CoefficientList[Series[(Sqrt[(1-2*x)*(1-2*x-4*x^2+4*x^3)]-2*x^2+4*x-1)/( 2 x*(1 - 4 x + 3 x^2)), {x, 0, 50}], x] (* G. C. Greubel, Sep 22 2018 *)
PROG
(PARI) my(x='x+O('x^50)); Vec((sqrt((1-2*x)*(1-2*x-4*x^2+4*x^3))-2*x^2+4*x-1)/(2*x*(1-4*x+3*x^2))) \\ G. C. Greubel, Sep 22 2018
(Magma) m:=50; R<x>:=PowerSeriesRing(Rationals(), m); Coefficients(R!((Sqrt((1-2*x)*(1-2*x-4*x^2+4*x^3))-2*x^2+4*x-1)/(2*x*(1-4*x+3*x^2)))); // G. C. Greubel, Sep 22 2018
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Mar 04 2010
STATUS
approved