OFFSET
0,1
EXAMPLE
0.489635226684303373081541660578468619322416625...
MATHEMATICA
val = -s*Log[r*s] / Sqrt[2*Pi*((-2 - 3*Log[r*s] + 2*Log[1 - r*s])* QPolyGamma[0, 1, r*s] + QPolyGamma[0, 1, r*s]^2 - 4*ArcTanh[1 - 2*r*s]*(Log[r*s] - Log[1 - r*s]/2 - r*(s/(1 - r*s))) - 2*(Log[1 - r*s]/(1 - r*s)) - QPolyGamma[1, 1, r*s] + r*s*Log[r* s]*((-r)*s^2*Log[r*s]* Derivative[0, 2][QPochhammer][r*s, r*s] + 2*Derivative[0, 0, 1][QPolyGamma][0, 1, r*s]))] /. FindRoot[{s == 1/QPochhammer[r*s], 1/s + r*s*Derivative[0, 1][QPochhammer][r*s, r*s] == (Log[1 - r*s] + QPolyGamma[0, 1, r*s])/(s* Log[r*s])}, {r, 1/5}, {s, 1}, WorkingPrecision -> 1000]; RealDigits[Chop[val], 10, -Floor[Log[10, Abs[Im[val]]]] - 3][[1]]
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Vaclav Kotesovec, Sep 26 2023
STATUS
approved