%I #12 Sep 13 2018 08:04:37
%S 6,7,2,8,3,3,2,4,6,5,5,3,1,6,6,6,0,7,9,8,8,4,6,6,0,9,4,2,3,6,6,1,6,0,
%T 6,3,8,2,1,9,8,0,1,8,5,6,5,6,7,6,2,1,1,6,0,1,8,2,6,7,1,5,5,8,1,7,4,9,
%U 1,6,1,9,7,0,7,0,9,6,1,8,3,7,1,5,4,4
%N Decimal expansion of the middle x such that 1/x + 1/(x+1) + 1/(x+3) = 2.
%C Equivalently, the least root of 2*x^3 + 5*x^2 - 2*x - 3;
%C Middle root: A316165;
%C Greatest root: A316166.
%C See A305328 for a guide to related sequences.
%F greatest root: -(5/6) + (1/6) sqrt(37) cos((1/3)(Pi - arctan((6 sqrt(1329))/53))) + (1/6) sqrt(37) cos((1/3)(-Pi + arctan((6 sqrt(1329))/53)))
%F ****
%F middle: -(5/6) - (1/12) sqrt(37) cos((1/3)(Pi - arctan((6 sqrt(1329))/53))) -
%F (1/12) sqrt(37) cos((1/3)(-Pi + arctan((6 sqrt(1329))/53))) +
%F (1/4) sqrt(37/3) sin((1/3)(Pi - arctan((6 sqrt(1329))/53))) -
%F (1/4) sqrt(37/3) sin((1/3)(-Pi + arctan((6 sqrt(1329))/53)))
%F ****
%F least: -(5/6) - (1/12) sqrt(37) cos((1/3)(Pi - arctan((6 sqrt(1329))/53))) -
%F (1/12) sqrt(37) cos(1/3(-Pi + arctan((6 sqrt(1329))/53))) -
%F (1/4) sqrt(37/3) sin((1/3)(Pi - arctan((6 sqrt(1329))/53))) +
%F (1/4) sqrt(37/3) sin(1/3)(-Pi + arctan((6 sqrt(1329))/53)))
%e greatest root: 0.83684889130097120054...
%e middle root: -0.67283324655316660799...
%e least root: -2.6640156447478045925...
%t a = 1; b = 1; c = 1; u = 0; v = 1; w = 3; d = 2;
%t r[x_] := a/(x + u) + b/(x + v) + c/(x + w);
%t t = Re[x /. ComplexExpand[Solve[r[x] == d, x]]]
%t N[t, 20]
%t u = N[t, 200];
%t u1 = RealDigits[u[[1]]] (* A316166, greatest *)
%t u2 = RealDigits[u[[2]]] (* A316164, least *)
%t u3 = RealDigits[u[[3]]] (* A316165, middle *)
%o (PARI) solve(x=-1, 0, 2*x^3+5*x^2-2*x-3) \\ _Altug Alkan_, Aug 27 2018
%Y Cf. A305328, A316164, A316166.
%K nonn,cons
%O 0,1
%A _Clark Kimberling_, Aug 08 2018