login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Integer part of Gamma(n+4/9)/Gamma(4/9).
7

%I #23 Sep 08 2022 08:44:44

%S 1,0,0,1,5,24,130,842,6274,52988,500442,5226841,59818300,744405515,

%T 10008118594,144561713037,2232675345794,36715105686390,

%U 640474621418148,11813198572823622,229701083360459334,4696111037591613052

%N Integer part of Gamma(n+4/9)/Gamma(4/9).

%H G. C. Greubel, <a href="/A020067/b020067.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: seq(f(n,4/9), n=0..25);

%p # alternative:

%p seq(trunc(pochhammer(4/9,n)), n = 0..25); # _G. C. Greubel_, Nov 13 2019

%t IntegerPart[Pochhammer[4/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, 4/9)) ) \\ _G. C. Greubel_, Nov 13 2019

%o (Magma) [Truncate(Gamma(n+4/9)/Gamma(4/9)): n in [0..25]]; // _G. C. Greubel_, Nov 13 2019

%o (Sage) [int(rising_factorial(4/9, n)) for n in (0..25)] # _G. C. Greubel_, Nov 13 2019

%Y Cf. A020064, A020065, A020066, A020068, A020069.

%Y Cf. A020022, A020112.

%K nonn

%O 0,5

%A _Simon Plouffe_