OFFSET
0,4
COMMENTS
LINKS
Ralf Stephan, Table of n, a(n) for n = 0..999
FORMULA
0 = b(n)*(+b(n+1) - 4*b(n+2) + 4*b(n+3)) + b(n+1)*(-2*b(n+1) + 9*b(n+2) - 10*b(n+3)) + b(n+2)*(-2*b(n+2) + 4*b(n+3)) if n>=0. - Michael Somos, May 26 2014
b(n) = -Re(Phi(2, 1, n + 1)). - Eric W. Weisstein, Dec 11 2017
G.f. for b(n): -log(1-x)/(2*(1-x/2)). - Vladimir Kruchinin, Nov 14 2022
EXAMPLE
0, 1/2, 1/2, 5/12, 1/3, 4/15, 13/60, 151/840, 16/105, 83/630, 73/630, ...
MATHEMATICA
Table[-Re[LerchPhi[2, 1, n + 1]], {n, 0, 20}] // Numerator (* Eric W. Weisstein, Dec 11 2017 *)
-Re[LerchPhi[2, 1, Range[20]]] // Numerator (* Eric W. Weisstein, Dec 11 2017 *)
RecurrenceTable[{b[n] == b[n - 1]/2 + 1/(2 n), b[0] == 0}, b[n], {n, 20}] // Numerator (* Eric W. Weisstein, Dec 11 2017 *)
PROG
(Sage)
def a():
b = n = 0
while True:
yield numerator(b)
n = n + 1
b = (b/2 + 1/(2*n)) # Ralf Stephan, May 18 2014
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Paul Curtz, May 12 2014
EXTENSIONS
a(14)-a(25) from Jean-François Alcover, May 12 2014
Corrected a(22) and a(24), more terms from Ralf Stephan, May 18 2014
STATUS
approved