A316136 Decimal expansion of the greatest x such that 1/x + 1/(x+2) + 1/(x+3) = 1.

1, 8, 6, 6, 1, 9, 8, 2, 6, 2, 5, 0, 9, 0, 2, 2, 5, 0, 5, 4, 7, 3, 9, 9, 4, 2, 3, 0, 9, 8, 6, 5, 9, 1, 4, 0, 6, 2, 3, 5, 9, 0, 0, 7, 5, 5, 7, 7, 2, 8, 8, 0, 4, 6, 1, 9, 7, 8, 4, 1, 6, 7, 0, 4, 8, 6, 1, 1, 1, 7, 0, 7, 6, 1, 0, 8, 7, 1, 6, 1, 9, 7, 8, 2, 6, 3, 9, 5, 0, 3, 0, 6, 6, 1, 2, 8, 2, 2, 5, 0, 2, 8, 0, 0, 0

Offset

1,2

Comments

Equivalently, the greatest root of x^3 + 2*x^2 - 4*x - 6;

Middle root: A316135;

Greatest root: A316136.

See A305328 for a guide to related sequences.

Links

Table of n, a(n) for n=1..105.

Formula

greatest root: -(2/3) + 8/3 cos[1/3 arctan[(3 sqrt[303])/37]]

middle: -(2/3) - 4/3 cos[1/3 arctan[(3 sqrt[303])/37]] + (4 sin[1/3 arctan[(3 sqrt[303])/37]])/sqrt[3]

least: -(2/3) - 4/3 cos[1/3 arctan[(3 sqrt[303])/37]] - (4 sin[1/3 arctan[(3 sqrt[303])/37]])/sqrt[3]

Example

greatest root: 1.8661982625090225055...
middle root: -1.2107558809591917224...
least root: -2.6554423815498307831...

Mathematica

a = 1; b = 1; c = 1; u = 0; v = 2; w = 3; d = 1;
r[x_] := a/(x + u) + b/(x + v) + c/(x + w);
t = x /. ComplexExpand[Solve[r[x] == d, x]]
N[t, 20]
u = N[t, 200];
RealDigits[u[[1]]]  (* A316134, least *)
RealDigits[u[[2]]]  (* A316135, middle *)
RealDigits[u[[3]]]  (* A316136, greatest *)
RealDigits[8 Cos[ArcTan[3 Sqrt[303]/37]/3]/3 - 2/3, 10, 111][[1]] (* Robert G. Wilson v, Jul 21 2018 *)

Crossrefs

Cf. A305328, A316134, A316135.

Sequence in context: A010527 A270137 A269846 * A102887 A067970 A003675

Adjacent sequences: A316133 A316134 A316135 * A316137 A316138 A316139

Keyword

nonn,cons

Author

Clark Kimberling, Jul 21 2018

Status

approved