OFFSET
0,1
COMMENTS
Not equal to exp(-4/9), which agrees with the first 16 decimal places. Related to Jacobi theta constant theta_2 and Dedekind's eta(x^2)^2/eta(x): Sum_{n>=0} x^(n*(n+1)/2) = 1.9873697... (A106334). This constant divided by constant A106334 equals constant A106335, the radius of convergence of the g.f. of A106336.
FORMULA
Sum_{n>=0} (1 - n*(n+1)/2)*x^(n*(n+1)/2) = 0.
EXAMPLE
0 = 1 - 2*x^3 - 5*x^6 - 9*x^10 - 14*x^15 - 20*x^21 - 27*x^28 - ...
x=0.641180388429954579645644888628301106553419618910071190877560305051317278
MATHEMATICA
digits = 105; g[x_?NumericQ] := NSum[(1 - n*(n+1)/2)*x^(n*(n+1)/2), {n, 0, Infinity}, WorkingPrecision -> digits+5, NSumTerms -> 100]; x /. FindRoot[g[x], {x, 1/2}, WorkingPrecision -> digits+5] // RealDigits[#, 10, digits]& // First (* Jean-François Alcover, Feb 12 2013 *)
PROG
(PARI) solve(x=.6, .7, sum(n=0, 100, (1-n*(n+1)/2)*x^(n*(n+1)/2)))
CROSSREFS
KEYWORD
cons,nonn
AUTHOR
Paul D. Hanna, Apr 29 2005
STATUS
approved