login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A173986 a(n) = numerator((Psi(1, 2/3) - Psi(1, n+2/3))/9), where Psi(1, z) is the Trigamma function. 6
0, 1, 29, 489, 60769, 3026081, 884023809, 890877733, 474015890357, 80471258049933, 67921427083803253, 1089963588226225073, 1092655876391630769, 395273284628034202009, 665644988593672027490729 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

a(n+1)/A173987(n+1) gives, for n >= 0, the partial sum Sum_{k=0..n} 1/(3*k+2)^2. The limit n -> infinity is given in A294967 as the Hurwitz Zeta function or the Trigamma function (1/9)*Zeta(2, 2/3) = (1/9)*Psi(1, 2/3) = 0.3404306010 ... - Wolfdieter Lang, Nov 12 2017

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..300

FORMULA

a(n) = numerator(r(n)) with r(n) = (1/9)*(4*(Pi^2)/3 - Zeta(2, 1/3) - Zeta(2, (3*n+2)/3)) = (1/9)*(Zeta(2, 2/3) - Zeta(2, (3*n+2)/3)) with the Hurwitz Zeta function Zeta(2, q). This becomes the formula given in the name. - Wolfdieter Lang, Nov 13 2017

a(n) = numerator of (1/9)*(2(Pi^2)/3 - J - Zeta(2, (3n+2)/3)) where J is the constant A173973 [which becomes the preceding formula].

a(n) = numerator of Sum_{k=0..(n-2)} 2/(3*k+2)^2. - G. C. Greubel, Aug 23 2018

EXAMPLE

The rationals a(n)/A173987(n) begin 0/1, 1/4, 29/100, 489/1600, 60769/193600, 3026081/9486400, 884023809/2741569600, 890877733/2741569600, ... - Wolfdieter Lang, Nov 12 2017

MAPLE

r := n -> (Psi(1, 2/3) - Psi(1, n+2/3))/9:

seq(numer(simplify(r(n))), n=0..14); # Peter Luschny, Nov 13 2017

MATHEMATICA

Table[Numerator[FunctionExpand[(4*Pi^2/3 - Zeta[2, 1/3] - Zeta[2, (3*n + 2)/3])/9]], {n, 0, 20}] (* Vaclav Kotesovec, Nov 14 2017 *)

Numerator[Table[Sum[2/(3*k + 2)^2, {k, 0, n - 2}], {n, 1, 20}]] (* G. C. Greubel, Aug 23 2018 *)

PROG

(PARI) for(n=1, 20, print1(numerator(sum(k=0, n-2, 2/(3*k+2)^2)), ", ")) \\ G. C. Greubel, Aug 23 2018

(MAGMA) [0] cat [Numerator((&+[2/(3*k+2)^2: k in [0..n-2]])): n in [2..20]]; // G. C. Greubel, Aug 23 2018

CROSSREFS

For denominators see A173987.

For 9*A173986 see A173985.

Cf. A006752, A120268, A173945, A173947, A173948, A173949, A173953, A173955, A173973, A173982, A173983, A173984, A173985, A294967.

Sequence in context: A022657 A182014 A261540 * A258462 A320550 A211833

Adjacent sequences:  A173983 A173984 A173985 * A173987 A173988 A173989

KEYWORD

frac,nonn,easy

AUTHOR

Artur Jasinski, Mar 04 2010

EXTENSIONS

Name simplified by Peter Luschny, Nov 13 2017

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 18 02:57 EST 2020. Contains 332006 sequences. (Running on oeis4.)