login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A229555
An exotic continued fraction for the real root of 6y^3 + 4y^2 - 4y - 7.
1
1, 22, 1, 31, 2, 3, 1, 63, 1, 10, 1, 2, 1, 7, 1, 160905, 2, 1, 4, 58, 2, 2, 1, 2, 1, 7, 3, 1, 3, 1, 4, 3, 1, 47, 1, 214540, 1, 2, 9, 1, 45, 1, 3, 1, 48, 1, 21, 1, 9, 1, 8, 1, 2, 249610, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 20, 1, 4, 19, 1, 2, 1, 1, 1, 1, 3, 4, 1, 1, 1
OFFSET
0,2
COMMENTS
Among its 148 initial numbers, this sequence has 8 unexpectedly big terms which are considerably bigger than those of the well-known exotic sequence A002937. It is peculiar that the ratio of the biggest of these denominators for the two examples, as well as the ratio of the smallest ones, is almost exactly 7: 115270760/16467250 = 7.000000607... and 160905/22986 = 7.000130514...
Deleting the first term in the sequence, we get a continued fraction for the algebraic integer given as the real root of the equation x^3 - 22*x^2 - 22*x - 6 = 0. It has the same set of big denominators, only shifted one position towards the beginning. - Sergei Duzhin, Oct 03 2013
For all integer cubic polynomials a*y^3 + b*y^2 + c*y + d with 0 < a <= 7, |b| <= 7, |c| <= 7, |d| <= 7, the given one (together with its three modifications obtained by a change of one sign and a palindromic reversal of coefficients) is the only polynomial in this set that, among its first 200 denominators, contains a number greater than 10^8, thus establishing a record. - Sergei Duzhin, Oct 04 2013
FORMULA
y = ((2906 - 126*sqrt(3*163))^(1/3) + (2906 + 126*sqrt(3*163))^(1/3) - 4) / 18. - Andrey Zabolotskiy, Jan 21 2023
MAPLE
Digits:=500:
with(numtheory):
x:=fsolve(6*y^3 + 4*y^2 - 4*y - 7);
cfrac(x, 200, 'quotients');
MATHEMATICA
r = Roots[6 x^3 + 4 x^2 - 4 x - 7 == 0, x][[1, 2]]; ContinuedFraction[r, 115] (* T. D. Noe, Oct 02 2013 *)
PROG
(PARI)
\p 250
contfrac(real(polroots(Pol([6, 4, -4, -7]))[1])) \\ Charles R Greathouse IV, Oct 01 2013
CROSSREFS
Sequence in context: A040502 A040503 A040504 * A040505 A040478 A040477
KEYWORD
nonn,cofr,easy
AUTHOR
Sergei Duzhin, Oct 01 2013
STATUS
approved