|
%I #5 Mar 30 2012 18:57:57
%S 3,3,6,4,8,2,7,0,1,9,2,3,3,5,2,8,1,5,7,7,0,3,9,4,9,3,7,6,1,1,0,6,7,7,
%T 8,1,4,4,3,6,5,3,0,1,1,7,8,4,0,0,3,6,7,9,4,6,8,5,6,3,5,3,2,4,2,5,3,4,
%U 9,3,1,1,2,9,0,3,6,8,3,7,2,5,6,4,9,3,2,1,7,3,9,8,2,0,0,1,7,2,7
%N Decimal expansion of x>0 satisfying 3*x^2+2*sin(x)=1.
%C See A198866 for a guide to related sequences. The Mathematica program includes a graph.
%e negative: -0.93194453919657480875799482221903577743...
%e positive: 0.336482701923352815770394937611067781443...
%t a = 3; b = 2; c = 1;
%t f[x_] := a*x^2 + b*Sin[x]; g[x_] := c
%t Plot[{f[x], g[x]}, {x, -2, 2}, {AxesOrigin -> {0, 0}}]
%t r = x /. FindRoot[f[x] == g[x], {x, -.94, -.93}, WorkingPrecision -> 110]
%t RealDigits[r] (* A199152 *)
%t r = x /. FindRoot[f[x] == g[x], {x, .33, .34}, WorkingPrecision -> 110]
%t RealDigits[r] (* A199153 *)
%Y Cf. A198866.
%K nonn,cons
%O 0,1
%A _Clark Kimberling_, Nov 03 2011
|
