OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
FORMULA
a(n) ~ -(2/3)*Gamma(1/3)^-1*n^(-5/3)*3^(2*n)*(1 + (5/9)*n^-1 + ...).
G.f.: (1-9*x)^(2/3).
D-finite with recurrence: n*a(n) +3*(-3*n+5)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
Sum_{n>=0} 1/a(n) = 99/128 - 5*sqrt(3)*Pi/512 - 15*log(3)/512. - Amiram Eldar, Dec 02 2022
From Peter Bala, Oct 14 2024: (Start)
a(n) = -1/(sqrt(3)*Pi) * 9^n * Gamma(2/3)*Gamma(n-2/3)/Gamma(n+1).
For n >= 1, a(n) = - Integral_{x = 0..9} x^n * w(x) dx, where w(x) = sqrt(3)/(2*Pi) * x^(1/3)*(9 - x)^(2/3)/x^2. (End)
MAPLE
a:= n-> (3^n/n!)*product(3*k-2, k=0..n-1); seq(a(n), n=0..25); # G. C. Greubel, Aug 22 2019
MATHEMATICA
Table[9^n*Pochhammer[-2/3, n]/n!, {n, 0, 25}] (* G. C. Greubel, Aug 22 2019 *)
PROG
(PARI) a(n)=if(n<0, 0, prod(k=0, n-1, 3*k-2)*3^n/n!)
(Magma) [1] cat [3^n*(&*[3*k-2: k in [0..n-1]])/Factorial(n): n in [1..25]]; // G. C. Greubel, Aug 22 2019
(Sage) [9^n*rising_factorial(-2/3, n)/factorial(n) for n in (0..25)] # G. C. Greubel, Aug 22 2019
(GAP) List([0..25], n-> 3^n*Product([0..n-1], k-> 3*k-2)/Factorial(n) ); # G. C. Greubel, Aug 22 2019
CROSSREFS
KEYWORD
sign,easy,changed
AUTHOR
Joe Keane (jgk(AT)jgk.org)
STATUS
approved