|
%I #5 Mar 30 2012 18:57:57
%S 1,4,0,7,1,0,6,6,3,9,7,6,3,0,3,9,9,0,5,9,9,6,1,6,6,8,8,9,0,3,5,9,5,7,
%T 1,2,8,5,0,8,7,2,6,3,1,0,5,4,0,1,0,4,5,4,1,3,6,1,7,6,9,4,2,3,1,7,2,2,
%U 7,3,2,3,1,9,6,9,6,4,2,7,3,6,7,4,7,2,9,1,3,6,5,2,1,3,2,7,1,6,5
%N Decimal expansion of x<0 satisfying 2*x^2+3*sin(x)=1.
%C See A198866 for a guide to related sequences. The Mathematica program includes a graph.
%e negative: -1.40710663976303990599616688903595712850...
%e positive: 0.28352555144412274861845117327862135032...
%t a = 2; b = 3; 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, -1.5, -1.4}, WorkingPrecision -> 110]
%t RealDigits[r] (* A199071 *)
%t r = x /. FindRoot[f[x] == g[x], {x, .28, .29}, WorkingPrecision -> 110]
%t RealDigits[r] (* A199072 *)
%Y Cf. A198866.
%K nonn,cons
%O 1,2
%A _Clark Kimberling_, Nov 02 2011
|
