%I #20 Sep 08 2022 08:44:44
%S 1,0,0,0,0,3,17,104,740,6005,54717,553253,6147263,74450193,976124760,
%T 13774204958,208143541601,3353423725801,57380805974819,
%U 1039230152655057,19860842917407765,399423618672311729
%N Integer part of Gamma(n+1/9)/Gamma(1/9).
%H G. C. Greubel, <a href="/A020069/b020069.txt">Table of n, a(n) for n = 0..445</a>
%p Digits := 64:f := proc(n,x) trunc(GAMMA(n+x)/GAMMA(x)); end;
%p seq(trunc(pochhammer(1/9,n)), n = 0..25); # _G. C. Greubel_, Nov 13 2019
%t IntegerPart[Pochhammer[1/9, Range[0, 25]]] (* _G. C. Greubel_, Nov 13 2019 *)
%o (PARI) P(n,x) = gamma(x+n)/gamma(x);
%o vector(26, n, truncate(P(n-1, 1/9)) ) \\ _G. C. Greubel_, Nov 13 2019
%o (Magma) [Truncate(Gamma(n+1/9)/Gamma(1/9)): n in [0..25]]; // _G. C. Greubel_, Nov 13 2019
%o (Sage) [int(rising_factorial(1/9, n)) for n in (0..25)] # _G. C. Greubel_, Nov 13 2019
%Y Cf. A020064, A020065, A020066, A020067, A020068.
%Y Cf. A020024, A020114.
%K nonn
%O 0,6
%A _Simon Plouffe_
|