A160106 Decimal representation of Bernays's number, 67^257^729.
3, 2, 5, 9, 1, 5, 3, 8, 4, 7, 9, 8, 6, 1, 8, 9, 9, 1, 6, 2, 5, 7, 1, 7, 8, 0, 6, 8, 8, 5, 0, 1, 5, 8, 8, 1, 1, 2, 0, 8, 6, 2, 1, 1, 5, 0, 5, 8, 7, 7, 8, 0, 7, 7, 6, 4, 5, 0, 7, 2, 8, 7, 0, 2, 5, 0, 4, 8, 0, 9, 7, 0, 2, 1, 2, 5, 0, 1, 7, 4, 1, 3, 3, 1, 7, 5, 6, 1, 7, 3, 6, 5, 3, 6, 4, 5, 8, 7, 4, 3, 2, 3, 5, 3, 7
Offset: 1
Links
- Robert G. Wilson v, Table of n, a(n) for n = 1..10000 (first 2223 digits from Lewis Mammel)
- P. Bernays, Sur le platonisme dans les mathématiques, L'Enseignement Mathématique 34 (1935), pp. 52-69. (See p. 61.)
- P. Bernays, Platonism in mathematics, translated from the French by C. D. Parsons.
- Efunda Engineering Fundamentals, Series Expansion of Exponential and Logarithmic Functions
- K. Podnieks, Towards a real Finitism?
Programs
-
Mathematica
nbrdgt = 105; f[base_, exp_] := RealDigits[ 10^FractionalPart[ N[exp*Log10[base], nbrdgt + Floor[ Log10[ exp]] + 2]], 10, nbrdgt][[1]]; f[67, 257^729] (* and the last 100 digits computed by PowerMod[67, 257^729, 10^100] *) (* Robert G. Wilson v, Aug 09 2016 *)
-
bc
/* bc script for cygwin bash shell or other UNIX environment */ /* Explicitly scale to 4000 decimal places */ i=10^4000 /* Calculate natural log of 100/67 */ b=100*i/67 c=i*(b-i)/(b+i) x=c s=0 for( j=1 ; x/j >0 ; j=j+2 ){ s = s + x/j x = x*c/i*c/i j /* progress mark */ } s=2*s /* Now s is the natural log of 100/67 */ /* Calculate natural log of sqrt(sqrt(10)) */ b=sqrt(10*i*i) b=sqrt(b*i) c=i*(b-i)/(b+i) x=c t=0 for( j=1 ; x/j >0 ; j=j+2 ){ t = t + x/j x = x*c/i*c/i j /* progress mark */ } t=2*t /* Now t is the integer part of 10^4000 * ln sqrt(sqrt(10)) */ ln10=4*t /* ln10 is the integer part of 10^4000 * ln 10 */ ln67 = 2*ln10 - s /* ln67 is the integer part of 10^4000 * ln 67 */ lg67 = ln67*i/ln10 /* lg67 is the integer part of 10^4000 * log_10 67 */ a=257^729 lgb = a*lg67 /* lgb is the integer part of 10^4000 * log_10 67^257^729 */ absc = lgb/i /* absc is the abscissa of lgb, and its value is one less than the number of decimal digits in Bernays's number */ mant = lgb - i*absc /* Find number of digits in abscissa */ x=absc for( nab=0 ; x>0 ; nab++ ) x = x/10 /* reduce the scale by nab */ mant = mant/10^nab ln10 = ln10/10^nab i = i/10^nab /* find ln 10^mant */ lnmant = mant*ln10/i /* calculate exp(lnmant) to get leading digits of Bernays's number */ fac=1 x=i n=0 for( j=0 ; x/fac > 0 ; j++ ){ n = n + x/fac x=x*lnmant/i fac = fac*(j+1) j /* progress mark */ } /* display abscissa of log_10 ( Bernays's number ) */ absc /* Display leading digits of Bernays's number. ( Truncation is to avoid displaying round-off error )*/ n / 10^20
Formula
Bernays's number is 67^257^729. The length and values of the sequence of its decimal representation is found by calculating its common logarithm by the formula, 257^729 * log_10(67) using an extended precision of 4000 digits. The number of digits of Bernays's number is given by the abscissa plus one, and the initial sequence is calculated from exponentiation of the mantissa.
Comments