cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A161771 Decimal expansion of (70*exp(Pi*sqrt(163)))^2.

Original entry on oeis.org

3, 3, 7, 7, 3, 6, 8, 7, 5, 8, 7, 6, 9, 3, 5, 4, 7, 1, 4, 6, 6, 3, 1, 9, 6, 3, 2, 5, 0, 6, 0, 2, 4, 4, 6, 3, 2, 0, 0, 0, 0, 0, 0, 0, 0, 8, 0, 2, 3, 1, 9, 3, 5, 6, 6, 2, 5, 2, 4, 9, 5, 7, 7, 1, 0, 4, 4, 1, 2, 4, 0, 6, 5, 9, 7, 4, 0, 9, 9, 7, 1, 0, 0, 6, 8, 5, 9, 8, 5, 1, 9, 3, 7, 0, 6, 5, 2, 2, 3, 2, 2, 8, 1, 6, 9
Offset: 39

Views

Author

Mark A. Thomas, Jun 18 2009

Keywords

Comments

Where exp^(Pi*sqrt163) is the Ramanujan constant and 70^2 is related to the norm vector 0 of the Leech lattice where 1^2 + 2^2 + 3^2 + ... + 22^2 + 23^2 + 24^2 = 70^2. A curiosity is: exp^2(Pi*sqrt163)*70^2 ~ hc/piGm^2 where all physics values are CODATA 2006 and m = neutron mass and exp^2(Pi*sqrt163)*70^2 = 3.377368...x 10^38 and hc/piGm^2 = 3.37700 x 10^38 (+- 0.00050) where 0.00050 = u_c which is the combined standard uncertainty.
This can also be expressed in a symmetric form in terms of the square of the neutron mass in units of Planck mass: where hc/2PiGm^2 = (Mp/m)^2 (Mp = Planck mass and m = neutron mass) and (exp^2(Pi*sqrt163)70^2)/2 ~ (Mp/m)^2. Note the divisor 2 in this case, which yields (exp^2(Pi*sqrt163)*70^2)/2 = 168868437938467735733159816253012231600.00000040115967. - Mark A. Thomas, Jul 02 2009

Examples

			337736875876935471466319632506024463200.00000080231935662524957710...
		

Crossrefs

Near relation to A160514 and A160515.

Programs

  • Magma
    SetDefaultRealField(RealField(100)); R:= RealField(); Exp(2*Pi(R)*Sqrt(163))*70^2; // G. C. Greubel, Oct 24 2018
  • Maple
    evalf((70*exp(Pi*sqrt(163)))^2,120); # Muniru A Asiru, Oct 25 2018
  • Mathematica
    First@ RealDigits[Exp[Pi Sqrt[163]]^2 70^2, 10, 105] (* Mark A. Thomas, Jun 18 2009, edited by Michael De Vlieger, Feb 19 2018 *)
  • PARI
    default(realprecision, 100); exp(2*Pi*sqrt(163))*70^2 \\ G. C. Greubel, Oct 24 2018
    

Formula

Equals exp(2*Pi*sqrt(163))*70^2.

A002937 An exotic continued fraction (for the real root of x^3 - 8*x - 10).

Original entry on oeis.org

3, 3, 7, 4, 2, 30, 1, 8, 3, 1, 1, 1, 9, 2, 2, 1, 3, 22986, 2, 1, 32, 8, 2, 1, 8, 55, 1, 5, 2, 28, 1, 5, 1, 1501790, 1, 2, 1, 7, 6, 1, 1, 5, 2, 1, 6, 2, 2, 1, 2, 1, 1, 3, 1, 3, 1, 2, 4, 3, 1, 35657
Offset: 0

Views

Author

Keywords

Comments

From Peter Bala, Oct 02 2013: (Start)
A theorem of Kuzmin in the measure theory of continued fractions says that for a random real number alpha, the probability that some given partial quotient of alpha is equal to a positive integer k is given by (1/log(2))*( log(1 + 1/k) - log(1 + 1/(k+1)) ). For example, almost all real numbers have 41.5% of their partial quotients equal to 1, 17% equal to 2, 9.3% equal to 3 and so on. Thus large partial quotients are the exception in continued fraction expansions.
Let now alpha denote the real root of Brillhart's cubic equation x^3 - 8*x - 10 = 0. Then alpha = (5 + (1/9)*sqrt(3*163))^(1/3) + (5 - (1/9)*sqrt(3*163))^(1/3) = 3.31862 82177 50185 65910 .... The continued fraction expansion of alpha is unusual in that there are 8 surprisingly large partial quotients in the first 200 terms, namely [22986, 1501790, 35657, 49405, 53460, 16467250, 48120, 325927]. The explanation for this behavior was given by Stark.
Roberts, p. 227, gives a remarkable product representation for beta := 2 + alpha, namely, beta = q^(1/24) * Product_{n >= 1} (1 + (1/q)^(2*n+1)), where q = exp(Pi*sqrt(163)), and notes that the first factor exp((Pi/24)*sqrt(163)) = 5.31862 82177 50185 63885 ... gives 16 decimal places of beta. Cf. A160514.
Some powers beta^k of beta share with beta the property of having exceptionally large partial quotients early on in their continued fraction expansion. Particularly noteworthy is beta^2, which, like beta, has 8 large partial quotients [126425, 8259853, 1620, 271730, 294038, 90569882, 264667, 1792603] in its first 200 terms. Strangely, apart from the third term 1620, these numbers are almost exactly 5.5 times larger than the corresponding large partial quotients [22986, 1501790, 35657, 49405, 53460, 16467250, 48120, 325927] occurring in the continued fraction expansion of beta.
Also noteworthy is beta^8 = 640320.00000000062437..., very nearly an integer, whose continued fraction expansion begins [640320, 1601600400, 320160, 2135467200, 261949, 10, 1, 20533337, 1, 1, 4, 3, 1, 3, 1369, 2, 2, 14, 3, 9535605, 1, 3, 2, 1, 2, 1, ...].
The cubic irrationals r*beta, where r is rational, also appear to have large partial quotients early on in their continued fraction expansions, which appear to be related to the partial quotients of beta. For example, (2/3)*beta has 8 exceptionally large partial quotients [137919, 9010749, 213951, 296433, 320769, 98803508, 288729, 1955567] in the first 200 terms of its continued fraction expansion and each of these partial quotients is almost exactly 6 times larger than the corresponding large partial quotient in the expansion of beta. (End)

Examples

			3.318628217750185659109680153... = 3 + 1/(3 + 1/(7 + 1/(4 + 1/(2 + ...)))).
		

References

  • R. P. Brent, A. J. van der Poorten, and H. J. Te Riele (1996, May), A comparative study of algorithms for computing continued fractions of algebraic numbers. In ANTS-II: International Algorithmic Number Theory Symposium (pp. 35-47), LNCS Vol. 1122, Springer, Berlin, Heidelberg. See Tables 3 and 4.
  • A. Ya. Khinchin, Continued Fractions, Dover Publications, 1997.
  • J. Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 227.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • H. M. Stark, An explanation of some exotic continued fractions found by Brillhart, pp. 21-35 of A. O. L. Atkin and B. J. Birch, editors, Computers in Number Theory. Academic Press, NY, 1971.

Crossrefs

Cf. A160332 (decimal expansion).

Programs

  • PARI
    { allocatemem(932245000); default(realprecision, 21000); x=NULL; p=x^3 - 8*x - 10; rs=polroots(p); r=real(rs[1]); c=contfrac(r); for (n=1, 20001, write("b002937.txt", n-1, " ", c[n])); } \\ Harry J. Smith, May 11 2009
    
  • PARI
    contfrac(polrootsreal(x^3 - 8*x - 10)[1]) \\ Charles R Greathouse IV, Apr 14 2014

A160515 Decimal expansion of the integer (640320^3 + 744)^2 * 70^2 = 337736875876935471466319632507953926400.

Original entry on oeis.org

3, 3, 7, 7, 3, 6, 8, 7, 5, 8, 7, 6, 9, 3, 5, 4, 7, 1, 4, 6, 6, 3, 1, 9, 6, 3, 2, 5, 0, 7, 9, 5, 3, 9, 2, 6, 4, 0, 0
Offset: 39

Views

Author

Mark A. Thomas, May 16 2009

Keywords

Comments

640320^3 + 744 is an integer that is very close to Ramanujan's constant e^(Pi*sqrt(163)) = 640320^3 + 743.99999999999925 and 70^2 is related to the norm vector 0 of the hyperbolic version of the Leech lattice since 1^2 + 2^2 + 3^3 + ... + 22^2 + 23^2 + 24^2 = 70^2.
337736875876935471466319632507953926400 = 2^8*3^2*5^2*7^2*10939058860032031^2 and 10939058860032031 is a prime number that can be decomposed to 2^15*3^2*5^3*23^3*29^3 + 31.

Programs

  • Mathematica
    IntegerDigits[(640320^3 + 744)^2*70^2] (* Michael De Vlieger, Dec 19 2015 *)

Formula

Equals Product_{i=1..5} A160514(i).

Extensions

Partially edited by R. J. Mathar, May 30 2009
Edited by N. J. A. Sloane, Dec 24 2010

A212131 Decimal expansion of k such that e^(k*sqrt(163)) = round(e^(Pi*sqrt(163))).

Original entry on oeis.org

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

Views

Author

Omar E. Pol, Jun 25 2012

Keywords

Comments

Decimal expansion of log(262537412640768744)/sqrt(163).
First differs from A000796 at a(32).
Note that 262537412640768744 = 24*10939058860032031 = 2^3 * 3 * 10939058860032031, is the nearest integer to the value of Ramanujan's constant e^(Pi*sqrt(163)) = 262537412640768743.999999999999250... = A060295.

Examples

			3.14159265358979323846264338327972661934754988... (very close to Pi).
		

Crossrefs

Programs

  • Mathematica
    RealDigits[Log[Round[E^(Pi Sqrt[163])]]/Sqrt[163], 10, 105][[1]] (* Bruno Berselli, Jun 26 2012 *)

Formula

k = log(round(e^(Pi*sqrt(163))))/sqrt(163).

Extensions

More terms from Alois P. Heinz, Jun 25 2012
Showing 1-4 of 4 results.