OFFSET
4,2
COMMENTS
Column k=10 of triangle A236830. - Philippe Deléham, Feb 02 2014
LINKS
G. C. Greubel, Table of n, a(n) for n = 4..1000
FORMULA
a(n) = A026848(n). - Philippe Deléham, Feb 02 2014
G.f.: (x^4*C(x)^10)/(1-x*C(x)^3) where C(x) is the g.f. of A000108. - Philippe Deléham, Feb 02 2014
D-finite with recurrence (n+4)*(3421*n+2687)*a(n) +(3421*n^2-245139*n-819238)*a(n-1) +3*(-126201*n^2+820641*n+1451992)*a(n-2) +(1944367*n^2-12105285*n+5094446)*a(n-3) +6*(-438489*n^2+3204217*n-6453730)*a(n-4) -12*(2*n-7)*(30601*n-111490)*a(n-5)=0. - R. J. Mathar, Jul 22 2025
MATHEMATICA
Drop[CoefficientList[Series[(1-Sqrt[1-4*x])^10/(128*x^4*(8*x^2 -(1 - Sqrt[1-4*x])^3 )), {x, 0, 40}], x], 4] (* G. C. Greubel, Jul 17 2019 *)
PROG
(PARI) my(x='x+O('x^40)); Vec((1-sqrt(1-4*x))^10/(128*x^4*(8*x^2 -(1 - sqrt(1-4*x))^3 ))) \\ G. C. Greubel, Jul 17 2019
(Magma) R<x>:=PowerSeriesRing(Rationals(), 40); Coefficients(R!( (1-Sqrt(1-4*x))^10/(128*x^4*(8*x^2 -(1-Sqrt(1-4*x))^3 )) )); // G. C. Greubel, Jul 17 2019
(SageMath) a=((1-sqrt(1-4*x))^10/(128*x^4*(8*x^2 -(1-sqrt(1-4*x))^3 ))).series(x, 45).coefficients(x, sparse=False); a[4:40] # G. C. Greubel, Jul 17 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved
