login
Iteration x := x - (sin(x) - cos(x)*x - Pi/2)/(sin(x)*x).
3

%I #37 Aug 18 2023 23:51:36

%S 1,9,0,5,6,9,5,7,2,9,3,0,9,8,8,3,8,9,4,8,8,2,6,6,6,4,3,7,1,6,0,9,6,6,

%T 7,0,3,4,9,5,0,4,3,1,2,1,6,1,2,8,0,3,2,1,2,1,9,3,5,6,4,5,5,9,9,9,4,5,

%U 4,4,2,4,0,9,9,5,7,9,5,0,2,2,7,5,7,1,6,1,6,6,3,4,6,4,6,3,0,3,9,7,1,5,3,9,8

%N Iteration x := x - (sin(x) - cos(x)*x - Pi/2)/(sin(x)*x).

%C Decimal expansion of psi, the unique solution on (0,Pi) of sin(psi) - psi*cos(psi) = Pi/2, an auxiliary constant used in the Hall-Tenenbaum inequality applied to real multiplicative functions. - _Jean-François Alcover_, Sep 05 2014

%H Steven R. Finch, <a href="http://arxiv.org/abs/2001.00578">Errata and Addenda to Mathematical Constants</a>, p. 28.

%H M. Fraser, <a href="http://www.jstor.org/stable/2690163">A tale of two goats</a>, Math. Mag., 55 (1982), 221-227.

%H Gerd Lamprecht, <a href="http://www.gerdlamprecht.de/Roemisch_JAVA.htm#ZZZZZ0004">Iterationsrechner mit Algorithmus</a>

%H Gerd Lamprecht, <a href="http://www.gerdlamprecht.de/langeZahlen/19056957293.pdf">10000 digits</a>

%H Gerd Lamprecht, <a href="http://www.gerdlamprecht.de/Zahlenfolgen.html">Zahlenfolgen (sequence)</a>

%F x := x - (sin(x) - cos(x)*x - Pi/2)/(sin(x)*x).

%F Equals 2*arccos(A133731/2).

%e A133731 = cos(1.9056957293.../2)*2.

%t psi = x /. FindRoot[Sin[x] - x*Cos[x] == Pi/2, {x, 2}, WorkingPrecision -> 102]; RealDigits[psi] // First (* _Jean-François Alcover_, Sep 05 2014 *)

%o (Other) Gerd Lamprecht Iterationsrechner: #(sin(x)-cos(x)*x-PI/2)/(sin(x)*x)@Na=2;@Nb=a;a=a-Fx(a);@N@Aa-b)%3C%204e-16@N0@N1@Nc=cos(a/2)*2;

%o (PARI) solve(x=1,2,sin(x)-x*cos(x)-Pi/2) \\ _Charles R Greathouse IV_, Mar 03 2021

%Y Cf. A072112, A133731.

%K cons,nonn

%O 1,2

%A Gerd Lamprecht (gerdlamprecht(AT)googlemail.com), Feb 12 2010