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