G. C. Greubel, <a href="/A203428/b203428_1.txt">Table of n, a(n) for n = 1..33</a>
G. C. Greubel, <a href="/A203428/b203428_1.txt">Table of n, a(n) for n = 1..33</a>
reviewed
approved
proposed
reviewed
editing
proposed
G. C. Greubel, <a href="/A203428/b203428_1.txt">Table of n, a(n) for n = 1..33</a>
a(n) = (-3)^binomial(n,2) * (Gamma(n+1))^(n-1) / BarnesG(n+1). - G. C. Greubel, Sep 28 2023
(* First program *)
f[j_] := 1/(3 *j); z = 16;
v[n_] := Product[Product[f[k] - f[j], {j, 1, k - 1}], {k, 2, n}]
1/Table[v[n], {n, 1, z}] (* A203428 *)
Table[v[n]/(3 *v[n + 1]), {n, 1, z - 1}] (* A203429 *)
(* Second program *)
Table[(-3)^Binomial[n, 2]*(Gamma[n+1])^(n-1)/BarnesG[n+1], {n, 20}] (* G. C. Greubel, Sep 28 2023 *)
(Magma)
Barnes:= func< n | (&*[Factorial(j): j in [1..n-1]]) >;
A203428:= func< n | (-3)^Binomial(n, 2)*(Factorial(n))^n/Barnes(n+1) >;
[A203428(n): n in [1..25]]; // G. C. Greubel, Sep 28 2023
(SageMath)
def barnes(n): return product(factorial(j) for j in range(n))
def A203428(n): return (-3)^binomial(n, 2)*(factorial(n))^n/barnes(n+1)
[A203428(n) for n in range(1, 21)] # G. C. Greubel, Sep 28 2023
approved
editing
_Clark Kimberling (ck6(AT)evansville.edu), _, Jan 02 2012
proposed
approved
editing
proposed
allocated for Clark KimberlingReciprocal of Vandermonde determinant of (1/3,1/6,...,1/(3n)).
1, -6, -486, 839808, 42515280000, -80335512599040000, -6890065294166289123840000, 31601087581187838970614157148160000, 8925080517850366815864624583251321642024960000
1,2
Each term divides its successor, as in A203429.
allocated
sign
Clark Kimberling (ck6(AT)evansville.edu), Jan 02 2012
approved
editing
allocated for Clark Kimberling
allocated
approved