OFFSET
0,1
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..444
FORMULA
E.g.f.: exp(x)*(37-30*x+9*x^2-x^3)/(1-x)^7.
a(n) = (n*(n^6 + 21*n^5 + 172*n^4 + 705*n^3 + 1522*n^2 + 1623*n + 653)*a(n-1) - (n^3 + 12*n^2 + 41*n + 37))/(n^6 + 15*n^5 + 82*n^4 + 207*n^3 + 244*n^2 + 105*n - 1), with a(0) = 37. - G. C. Greubel, Oct 08 2023
MATHEMATICA
a[n_]:= If[n<1, 37, (n*(n^6+21*n^5+172*n^4+705*n^3+1522*n^2+1623*n +653)*a[n-1] -(n^3+12*n^2+41*n+37))/(n^6+15*n^5+82*n^4+207*n^3 +244*n^2+105*n-1)];
Table[a[n], {n, 0, 40}] (* G. C. Greubel, Oct 08 2023 *)
PROG
(PARI) my(x='x+O('x^25)); Vec(serlaplace(exp(x)*(37-30*x+9*x^2-x^3)/(1-x)^7)) \\ Michel Marcus, Apr 06 2019
(Magma) R<x>:=PowerSeriesRing(Rationals(), 40); Coefficients(R!(Laplace( (37-30*x+9*x^2-x^3)*Exp(x)/(1-x)^7 ))); // G. C. Greubel, Oct 08 2023
(SageMath)
def A144497_list(prec):
P.<x> = PowerSeriesRing(QQ, prec)
return P( (37-30*x+9*x^2-x^3)*exp(x)/(1-x)^7 ).egf_to_ogf().list()
A144497_list(40) # G. C. Greubel, Oct 08 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
David Applegate and N. J. A. Sloane, Dec 13 2008
STATUS
approved