|
%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
|
