OFFSET
1,1
COMMENTS
Positive root of cosh(y) * sqrt(1 - (y/sinh(y))^2) - arccos(y/sinh(y)) = 4*Pi.
In general, all roots of sin(z) = z are given by z = 0 and z = +-(x_k)+-(y_k)*i, where x_k is the root of cos(x) * sqrt((x/sin(x))^2 - 1) - arccosh(x/sin(x)) = 0 in the range (2*k*Pi, (2*k+1/2)*Pi), y_k is the positive root of cosh(y) * sqrt(1 - (y/sinh(y))^2) - arccos(y/sinh(y)) = 2*k*Pi.
EXAMPLE
z = 13.8999597139... + 3.3522098848...*i is the unique root of sin(z) = z in the region {z: 4*Pi <= Re(z) < 6*Pi, Im(z) >= 0}.
MATHEMATICA
RealDigits[y /. FindRoot[{Re[Sin[x + I*y]] == x, Im[Sin[x + I*y]] == y}, {{x, 14}, {y, 3}}, WorkingPrecision -> 100], 10, 90][[1]] (* Amiram Eldar, Oct 10 2021 *)
PROG
(PARI) solve(y=3.3, 3.4, cosh(y) * sqrt(1 - (y/sinh(y))^2) - acos(y/sinh(y)) - 4*Pi)
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Jianing Song, Oct 10 2021
STATUS
approved