A316249 - OEIS

%I #7 Aug 22 2018 18:03:32

%S 2,7,4,5,9,9,3,4,7,5,7,9,4,0,2,7,5,2,2,5,1,6,4,4,1,7,0,0,1,6,3,4,0,3,

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

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

%N Decimal expansion of the least x such that 1/x + 1/(x+1) + 1/(x+3) = 3.

%C Equivalently, the least root of 3*x^3 + 9*x^2 + x - 3;

%C Middle root: A316250;

%C Greatest root: A316251.

%C See A305328 for a guide to related sequences.

%F greatest root: -1 + (2/3)*sqrt(2)*cos(Pi/3 - (1/3)*arctan(sqrt(431)/9))) + (2/3)*sqrt(2)*cos(-Pi/3 + (1/3)*arctan(sqrt(431)/9)))

%F ****

%F middle: -1 - (1/3)*sqrt(2)*cos(Pi/3 - (1/3)*arctan(sqrt(431)/9))) - (1/3)*sqrt(2)*cos(-Pi/3 + (1/3)*arctan(sqrt(431)/9))) + sqrt(2/3)*sin(Pi/3 - (1/3)*arctan(sqrt(431)/9))) - sqrt(2/3)*sin(-Pi/3 + (1/3)*arctan(sqrt(431)/9)))

%F ****

%F least: -1 - (1/3)*sqrt(2)*cos(Pi/3 - (1/3)*arctan(sqrt(431)/9))) - (1/3)*sqrt(2)*cos(-Pi/3 + (1/3)*arctan(sqrt(431)/9))) - sqrt(2/3)*sin(Pi/3 - (1/3)*arctan(sqrt(431)/9))) + sqrt(2/3)*sin(-Pi/3 + (1/3)*arctan(sqrt(431)/9)))

%e greatest root: 0.4896787901169438959...

%e middle root: -0.7436853143229163734...

%e least root: -2.745993475794027522...

%t a = 1; b = 1; c = 1; u = 0; v = 1; w = 3; d = 3;

%t r[x_] := a/(x + u) + b/(x + v) + c/(x + w);

%t t = x /. ComplexExpand[Solve[r[x] == d, x]]

%t N[t, 20]

%t y = Re[N[t, 200]];

%t RealDigits[y[[1]]] (* A316251, greatest *)

%t RealDigits[y[[2]]] (* A316249, least *)

%t RealDigits[y[[3]]] (* A316250, middle *)

%Y Cf. A305328, A316250, A316251.

%K nonn,cons

%O 1,1

%A _Clark Kimberling_, Aug 21 2018