%I #17 Sep 08 2022 08:44:45
%S 1,0,0,1,3,14,79,505,3710,30920,288591,2982109,33797240,416832633,
%T 5557768440,79661347641,1221473997164,19950741953681,345812860530487,
%U 6339902443058941,122571447232472865,2492286093726948268
%N Integer part of Gamma(n + 1/3)/Gamma(1/3).
%H G. C. Greubel, <a href="/A020089/b020089.txt">Table of n, a(n) for n = 0..449</a>
%p Digits:= 64: x:=1/3: f:= proc(n,x) trunc(GAMMA(n+x)/GAMMA(x)); end;
%p # alternative
%p f:= n -> trunc(simplify(GAMMA(n+1/3)/GAMMA(1/3))): map(f, [$0..30]); # _Robert Israel_, Nov 20 2018
%t Table[IntegerPart[Pochhammer[1/3, n]], {n, 0, 30}] (* _G. C. Greubel_, Nov 19 2018 *)
%o (PARI) default(realprecision, 100); for(n=0,30, print1(truncate( gamma(n+1/3)/gamma(1/3)), ", ")) \\ _G. C. Greubel_, Nov 19 2018
%o (Magma) [Truncate(Gamma(n + 1/3)/Gamma(1/3)): n in [0..30]]; // _G. C. Greubel_, Nov 19 2018
%o (Sage) [int(gamma(n + 1/3)/gamma(1/3)) for n in range(30)] # _G. C. Greubel_, Nov 19 2018
%Y Cf. A020044, A020134.
%K nonn
%O 0,5
%A _Simon Plouffe_
|