OFFSET
1,2
LINKS
FORMULA
7*a(n) = (9*n-2)(!^9) := Product_{j=1..n} (9*j-2).
E.g.f.: (-1+(1-9*x)^(-7/9))/7.
D-finite with recurrence: a(1) = 1, a(n) = (9*n - 2)*a(n-1) for n > 1. - Georg Fischer, Feb 15 2020
a(n) = (1/7) * 9^n * Pochhammer(n, 7/9). - G. C. Greubel, Oct 19 2022
From Amiram Eldar, Dec 21 2022: (Start)
a(n) = A147631(n+1)/7.
Sum_{n>=1} 1/a(n) = 7*(e/9^2)^(1/9)*(Gamma(7/9) - Gamma(7/9, 1/9)). (End)
MAPLE
f := gfun:-rectoproc({(9*n - 2)*a(n - 1) - a(n) = 0, a(1) = 1}, a(n), remember);
map(f, [$ (1 .. 20)]); # Georg Fischer, Feb 15 2020
MATHEMATICA
Table[9^n*Pochhammer[7/9, n]/7, {n, 40}] (* G. C. Greubel, Oct 19 2022 *)
PROG
(Magma) [n le 1 select 1 else (9*n-2)*Self(n-1): n in [1..40]]; // G. C. Greubel, Oct 19 2022
(SageMath) [9^n*rising_factorial(7/9, n)/7 for n in range(1, 40)] # G. C. Greubel, Oct 19 2022
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
EXTENSIONS
Terms a(15) onward added by G. C. Greubel, Oct 19 2022
STATUS
approved