%I #8 Apr 10 2021 08:05:23
%S 1,4,1,3,9,2,2,5,4,0,9,0,9,2,9,6,7,4,0,4,2,4,4,5,3,3,3,3,0,3,6,0,3,3,
%T 1,1,3,0,4,0,9,0,1,9,1,5,7,1,0,0,0,8,3,1,5,0,5,5,0,3,1,6,0,0,5,8,0,6,
%U 3,7,8,3,6,7,5,4,0,2,7,3,0,1,2,4,9,0,2,5,7,2,8,1,9,1,2,2,6,1,8,7
%N Decimal expansion of the least x > 0 satisfying 1 = x*sin(x - Pi/5).
%e x=1.41392254090929674042445333303603311304090191571000...
%t Plot[{1/x, Sin[x], Sin[x - Pi/2], Sin[x - Pi/3], Sin[x - Pi/4]}, {x,
%t 0, 2 Pi}]
%t t = x /. FindRoot[1/x == Sin[x], {x, 1, 1.2}, WorkingPrecision -> 100]
%t RealDigits[t] (* A133866 *)
%t t = x /. FindRoot[1/x == Sin[x - Pi/2], {x, 1, 2}, WorkingPrecision -> 100]
%t RealDigits[t] (* A196767 *)
%t t = x /. FindRoot[1/x == Sin[x - Pi/3], {x, 1, 2}, WorkingPrecision -> 100]
%t RealDigits[t] (* A196768 *)
%t t = x /. FindRoot[1/x == Sin[x - Pi/4], {x, 1, 2}, WorkingPrecision -> 100]
%t RealDigits[t] (* A196769 *)
%t t = x /. FindRoot[1/x == Sin[x - Pi/5], {x, 1, 2}, WorkingPrecision -> 100]
%t RealDigits[t] (* A196770 *)
%t t = x /. FindRoot[1/x == Sin[x - Pi/6], {x, 1, 2}, WorkingPrecision -> 100]
%t RealDigits[t] (* A196771 *)
%Y Cf. A196772.
%K nonn,cons
%O 1,2
%A _Clark Kimberling_, Oct 06 2011
|