## Ryan Matthew Ferris 10:26 PM Thursday, November 19, 2015, 2015 Bellingham, WA
## Updated 7/10/2018
## *Notes on 187.5/16 or 11.71875*
##
# In this code, I introduce 'xi' (11.71875) a constant with strong relationships to both pi and # exp(1).
# 'Xi' (11.71875 or (187.5/16)) is a 'translator constant' between powers of 2 and 16.
# Search google for "11.71875mhz" and you can see this. 11.71875mhz is a widely used frequency in 'crystal resonators'.
# More than simply being an efficient means to compute pi or exp(1) (e.g. 'Eulers number'), 'xi' may
# 'interpolate' between pi and e, giving science the means to measure logarithmic growth patterns of
# increasing radii. This may make it useful for gravitational and rotational calculations. It is
# possible that the 187.5 cm-3 pc represents some fixed attenuation limit (see
# http://arxiv.org/abs/1503.05245) for radio waves signal that require either boosting or
# chirping. The relationships between pi, e, and xi(187.5./16 = 11.71875) may represent fixed
# constants descriptive of all EM signaling. In addition, the constants that describe the
# circumscription of an octahedron (sqrt(1/2) and sqrt(1/6) can also be used to derive xi with a
# factor close to 9/16.
## A number of constants and function names are used. Some are reused:
# [1] "almost2" "aR" "c2" "cR" "eR" "F_pe" "F_pi" "findpi" "i" "j" "k" "k1"
# [13] "k2" "l" "l1" "mm1" "p1" "p2" "piR" "pR" "pX" "s4" "seqe" "seqpi"
# [25] "seqpi_e" "t1" "t2" "t3" "Vx" "x" "xi" "xi_sqrt" "xif" "y"
## Set your graphics device to record history if desired. 86 graphs are produced.
## You may want par(mfrow=c(1,2))
## NB : In R 3.5 there is no bigint by default,however there are 22 digits of accuracy.
## See my usage of BIGINT in Powershell
## For this work I used Microsoft Open R 3.5 with Intel's MKL (Math Kernel Library)
options(digits=22)
# main constants
pi
#[1] 3.1415926535897931
exp(1) # Eulers number
#[1] 2.7182818284590451
xi <- 187.5/16 # transdimensional constant
# [1] 11.71875
# xi with pi and exp(1)
pi + exp(1)
# [1] 5.8598744820488378
2 * (pi + exp(1)) # very close to xi
# [1] 11.719748964097676
(pi + exp(1)) / 11.71875 # very close .5
# [1] 0.5000426224681674414185
xi / (pi + exp(1)) # very close to 2
# [1] 1.9998295246594895
((xi /(pi + exp(1))) * (pi + exp(1)))
# [1] 11.71875
# xi without pi or exp(1)
# from regular octahedron constants for finding circumscribed sphere
sqrt(1/6)
# [1] 0.40824829046386302
sqrt(1/2)
# [1] 0.70710678118654757
904319/1600000 # something a little different than 9/16 ???
# 9/16 - 904319/1600000 = -0.0026993749999999483
((10 * (sqrt(1/2) + sqrt(1/6))) + 904319/1600000) # close to xi
# [1] 11.718750091504107
xi_sqrt <- ((10 * (sqrt(1/2) + sqrt(1/6))) + 904319/1600000)
xi_sqrt
# [1] 11.718750091504107
((xi_sqrt /(pi + exp(1))) * (pi + exp(1)))
# [1] 11.718750091504107
# one way of deriving almost pi with exp(1) and xi
xi /(pi + exp(1))
# [1] 1.9998295246594895
almost2 <- xi /(pi + exp(1)) # almost 2
((xi/2) - exp(1)) # two digits right of decimal of pi
# [1] 3.1410931715409549
((xi/almost2) - exp(1)) # 14 digits right of decimal point of pi
# [1] 3.1415926535897927
pi
# [1] 3.1415926535897931
# one way of deriving almost exp(1) with pi and xi
((xi/2) - pi) # two digits right of decimal of exp(1)
#[1] 2.7177823464102069
exp(1)
# [1] 2.7182818284590451
((xi/almost2) - pi) # 14 digits right of decimal point of exp(1)
# [1] 2.7182818284590446
exp(1)
# [1] 2.7182818284590451
# 187.5 is the first fractional member of a series that returns 11.71875 as
# the quotient of what appears to be integer multiples of (3,5,6,10) by powers of 2 (2^(1:16)).
# 11.71875 pops up in google queries for decimal translation of bytes. All of the results below
# equal 11.71875:
cat('
12000/1024
6000/512
3000/256
1500/128
750/64
375/32
187.5/16
46.875/4
23.4375/2
')
# 11.71875 can be seen as a type of hexadecimal 'translator' ratio between binary and decimal:
12000/1024
#[1] 11.71875
6000/512
#[1] 11.71875
3000/256
#[1] 11.71875
1500/128
#[1] 11.71875
750/(187.5/16)
#[1] 64
12000/(187.5/16)
#[1] 1024
6000/(187.5/16)
#[1] 512
3000/(187.5/16)
#[1] 256
1500/(187.5/16)
#[1] 128
750/(187.5/16)
#[1] 64
# However, I can also come up with something more transcendental:
options(digits=22)
(18.75/2)/10
#[1] 0.9375
(.9375/2)/10
#[1] 0.046875
(.046875/2)/10
#[1] 0.0023437499999999999
(.002343749999999999/2)/10
# 0.00011718749999999996
#or
(23.43749999999999/2)
# [1] 11.718749999999995
# ...
# Multiples of 187.5 can be produced by selecting out the products of vector (1:100) multiplication of the
# quotient of 187.5./16. By sixteen,eight,four:
data.frame((187.5/16 * 1:100)[c(seq(0,100,16))])
cat('
1 187.5
2 375.0
3 562.5
4 750.0
5 937.5
6 1125.0
')
# By eight:
cat('
data.frame((187.5/16 * 1:100)[c(seq(0,100,8))])
1 93.75
2 187.50
3 281.25
4 375.00
5 468.75
6 562.50
7 656.25
8 750.00
9 843.75
10 937.50
11 1031.25
12 1125.00
')
# By four:
data.frame((187.5/16 * 1:100)[c(seq(0,100,4))])
cat('
1 46.875
2 93.750
3 140.625
4 187.500
5 234.375
6 281.250
7 328.125
')
# This series or ladder is significant, but I am not quite sure why it accurately generates the range (1:1/10) to five digits.
l <- {};for(i in 1:10) {l <- rbind(l,(exp(pi^-11.71875)/i))};l;plot(l)
cat(' [,1]
[1,] 1.00000149288121576 #1
[2,] 0.50000074644060788 #1/2
[3,] 0.33333383096040525 #1/3
[4,] 0.25000037322030394 #1/4
[5,] 0.20000029857624316 #1/5
[6,] 0.16666691548020263 #1/6
[7,] 0.14285735612588796 #1/7
[8,] 0.12500018661015197 #1/8
[9,] 0.11111127698680175 #1/9
[10,] 0.10000014928812158 #1/10
')
pi # 22 digits R 3.2.2 constant for pi
#[1] 3.1415926535897931
exp(1) # Eulers number in R 3.2.2
#[1] 2.7182818284590451
# 'xi' a 'transdimensional' constant
187.5/16
#[1] 11.71875
xi <- 187.5/16 # 11.71875 - greater than pi or e
k1 <- (spline(sort(c(xi,exp(1),pi,-pi,-exp(1),-xi))))
k2 <- (sort(c(xi,exp(1),pi,-pi,-exp(1),-xi)))
splinefun(diff(k2))
plot(splinefun(diff(k2)))
plot(spline(sort(c(xi,exp(1),pi,-pi,-exp(1),-xi))))
plot(splinefun(sort(c(xi,exp(1),pi,-pi,-exp(1),-xi))))
xi <- (187.5/16)/10 # 1.171875 - less than pi or e
k1 <- (spline(sort(c(xi,exp(1),pi,-pi,-exp(1),-xi))))
k2 <- (sort(c(xi,exp(1),pi,-pi,-exp(1),-xi)))
splinefun(diff(k2))
plot(splinefun(diff(k2)))
plot(spline(sort(c(xi,exp(1),pi,-pi,-exp(1),-xi))))
plot(splinefun(sort(c(xi,exp(1),pi,-pi,-exp(1),-xi))))
pi #pi
#[1] 3.1415926535897931
exp(1) #e
#[1] 2.7182818284590451
187.5/16 #xi
#[1] 11.71875
((11.71875/2) - exp(1))
#[1] 3.1410931715409549 #close to pie
((11.71875/2) - pi)
#[1] 2.7177823464102069 #close to e
((11.71875/2) - exp(1)) + ((11.71875/2) - pi)
#[1] 5.8588755179511622 # close to .5(xi)
(((11.71875/2) - exp(1)) + ((11.71875/2) - pi)) * 2
#[1] 11.717751035902324 # close to xi
11.71875 - 11.717751035902324
exp(11.71875 - 11.717751035902324)
#[1] 1.0009994632285004 # close to 1
(23.43749999999999/2)
exp(11.71875 - 11.718749999999995)
#[1] 1.0000000000000053 # closer to 1
# Since (pi + exp(1)) * 2 ~= (187.5/16) it is possible to use 11.71875 and some
# adjustment to arrive at either pi or e. (See notes and functions seqpi,seqe,seqpi_e (far)below.)
# The products of pi and e doubled roughly equal (187.5/16) or 11.71875:
(((pi + exp(1)) * 2) /(187.5/16))
#[1] 1.0000852449363349 # nearly one again. irrational?
(((pi + exp(1))) /(187.5/16))
#[1] 0.50004262246816744 # nearly five. irrational?
# 11.71875 can be successfully used to calculate pi or e within 3 digits right of the decimal point without any compensation:
((187.5/16)/2) - pi # almost e...
#[1] 2.7177823464102069 # should be 2.7182818284590451
((187.5/16)/2) - exp(1) # almost pi
#[1] 3.1410931715409549 # should be 3.1415926535897931
# Not quite pi or e:
pi - (((187.5/16)/2) - exp(1))
#[1] 0.00049948204883820679
exp(1) - (((187.5/16)/2) - pi)
#[1] 0.00049948204883820679
# The adjustment ( tanh(11.71875) / 2000 ) gives us 5 digits of accuracy for pi or e:
tanh(11.71875) / 2000
[1] 0.00049999999993374431
eR <- ((11.71875/2) - pi + tanh(11.71875) / 2000 )
format(eR - exp(1),scientific=FALSE)
# [1] "0.0000005179510953468025"
eR
#[1] 2.7182823464101404
exp(1)
#[1] 2.7182818284590451
pR <- ((11.71875/2) - exp(1) + tanh(11.71875) / 2000 )
format(pR - pi,scientific=FALSE)
# [1] "0.0000005179510953468025"
pR
# [1] 3.1415931715408885
pi
# [1] 3.1415926535897931
# using -digamma(1)
xi - (gr/exp(1) + -digamma(1))
[1] -0.00058210447924356
# DOES xif makes xi a 'transdimenesional transcendant'?
(18.75/2)/10
#[1] 0.9375
(.9375/2)/10
#[1] 0.046875
(.046875/2)/10
#[1] 0.0023437499999999999
(.002343749999999999/2)/10
# 0.00011718749999999996
(18.75/2)/10
(.9375/2)/10
(.046875/2)/10
(.002343749999999999/2)/10
xif <- (.002343749999999999/2)/10
xif <- (xif * 100000)
# Powershell [BIGINT] will take me this far:
cat('
PS [bigint](.002343749999999999/2 * [MATH]::POW(10,32))
117187499999999960548787617792
PS [bigint](.002343749999999999/2 * [MATH]::POW(10,64))
11718749999999995682656043607972093960765263440537027795222528
PS[bigint](.002343749999999999/2 * [MATH]::POW(10,128))
117187499999999952444488198957480104184273770534709180501147236632873797389362215023609307534290876647973955371666657662468096
PS [bigint](.002343749999999999/2 * [MATH]::POW(10,256))
11718749999999995794845523228309099783974204074624886274076630993851149729639928280300079885632128025211313970818427611098306439684539273451270347793350886454282383635
265648724318526179750550535267831033255464007490486648205861676480655926487782004883456
')
tanh(11.71875) # helps get us close to 1
#[1] 0.99999999986748866
# so does
xif/xi
#[1] 9.9999999999999947
xi/xif
#[1] 0.10000000000000005
eR <- ((11.71875/2) - pi + tanh(11.71875) / 2000 )
format(eR - exp(1),scientific=FALSE)
# [1] "0.0000005179510953468025"
eR
# [1] 2.7182823464101404
exp(1)
#[1] 2.7182818284590451
pR <- ((11.71875/2) - exp(1) + tanh(11.71875) / 2000 )
format(pR - pi,scientific=FALSE)
# [1] "0.0000005179510953468025"
pR
# [1] 3.1415931715408885
pi
# [1] 3.1415926535897931
(18.75/2)/10
(.9375/2)/10
(.046875/2)/10
(.002343749999999999/2)/10
xif <- (.002343749999999999/2)/10
xif <- (xif * 100000)
eR <- ((xif/2) - pi + tanh(xif) / 2000 )
format(eR - exp(1),scientific=FALSE)
eR
exp(1)
pR <- ((xif/2) - exp(1) + tanh(xif) / 2000 )
format(pR - pi,scientific=FALSE)
pR
pi
## By selecting from an array we can get much more accurate pi and e from xi.
((11.71875/2) - pi + tanh((11.71875/4096) * 1:4096) / 2000 )[c(1322)]
#[1] 2.718281828165519
exp(1)
#[1] 2.7182818284590451
((11.71875/2) - exp(1) + tanh((11.71875/4096) * 1:4096) / 2000 )[c(1322)]
#[1] 3.141592653296267
pi
#[1] 3.1415926535897931
## It is then possible for pi or e to converge from series generated from xi (11.71875)
## This series will converge to pi :
c2 <- as.matrix(c(
(((11.71875/2) - exp(1)) + (tanh((11.71875/8) * 1:8)) / 2000 )[3],
(((11.71875/2) - exp(1)) + (tanh((11.71875/16) * 1:16)) / 2000 )[5],
(((11.71875/2) - exp(1)) + (tanh((11.71875/32) * 1:32)) / 2000 )[11],
(((11.71875/2) - exp(1)) + (tanh((11.71875/64) * 1:64)) / 2000 )[21],
(((11.71875/2) - exp(1)) + (tanh((11.71875/128) * 1:128)) / 2000 )[41],
(((11.71875/2) - exp(1)) + (tanh((11.71875/256) * 1:256)) / 2000 )[83],
(((11.71875/2) - exp(1)) + (tanh((11.71875/512) * 1:512)) / 2000 )[165],
(((11.71875/2) - exp(1)) + (tanh((11.71875/1024) * 1:1024)) / 2000 )[331],
(((11.71875/2) - exp(1)) + (tanh((11.71875/2048) * 1:2048)) / 2000 )[661],
(((11.71875/2) - exp(1)) + (tanh((11.71875/4096) * 1:4096)) / 2000 )[1322],
(((11.71875/2) - exp(1)) + (tanh((11.71875/8192) * 1:8192)) / 2000 )[2644],
(((11.71875/2) - exp(1)) + (tanh((11.71875/16384) * 1:16384)) / 2000 )[5288],
(((11.71875/2) - exp(1)) + (tanh((11.71875/32768) * 1:32768)) / 2000 )[10576],
(((11.71875/2) - exp(1)) + (tanh((11.71875/65536) * 1:65536)) / 2000 )[21154]
))
c2
cat(' [,1]
[1,] 3.1415930191734072
[2,] 3.1415925126008446
[3,] 3.1415928546514991
[4,] 3.1415927145681901
[5,] 3.1415926227939854
[6,] 3.1415926707780644
[7,] 3.1415926473344875
[8,] 3.1415926591902781
[9,] 3.1415926532962670
[10,] 3.1415926532962670
[11,] 3.1415926532962670
[12,] 3.1415926532962670
[13,] 3.1415926532962670
[14,] 3.1415926536666201
')
pi
# [1] 3.1415926535897931
# So we can generate series that converge around pi and e :
options(digits=22)
# These print very long sequences/series
# for(x in (2^seq(0,16,1))) {print(format(((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000))}
# for(x in (2^seq(0,16,1))) {print(format(((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000,digits=22,scientific=TRUE))}
# for(x in (2^seq(0,16,1))) {print(((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000)}
for(x in (2^seq(0,16,1))) {(((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000 )}
for(x in (2^seq(0,16,1))) {plot(((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / x ,pty=.25,col="blue",type="p")};abline(h=pi,col="red")
for(x in (2^seq(0,16,1))) {plot(((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000 ,pty=.25,col="blue",type="p")};abline(h=pi,col="red")
for(x in (2^seq(0,16,1))) {plot(((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000,ylim=c(pi - (pi *.000001),pi + (pi *.000001)),lwd=1,col="blue",type="l")};abline(h=pi,col="red")
for(x in (2^seq(8,12,1))) {plot(grep(3.141592,((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000,value=TRUE),lwd=1,col="blue",type="l")};abline(h=pi,col="red")
for(x in (2^seq(8,14,1))) {plot(grep(3.141592,((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000,value=TRUE),lwd=1,col="blue",type="l")};abline(h=pi,col="red")
findpi <- function(x,y) {for(i in (2^seq(x,y,1))) {p1 <- rbind(p1,(grep(pi,(11.71875/2) - exp(1) + (tanh((11.71875/i) * 1:i)) / 2000,value=TRUE)))}}
findpi <- function(x,y) {for(i in (2^seq(x,y,1))) {p1 <- rbind(p1,(grep(pi,(11.71875/2) - exp(1) + (tanh((11.71875/i) * 1:i)) / 2000,value=TRUE)))};plot(pi)}
p1 <- NULL; p1 <- rbind(p1,cbind(grep(3.1415,(11.71875/2) - exp(1) + (tanh((11.71875/8) * 1:8)) / 2000,value=TRUE)))
# template function code A:
x <- 1
y <- 16
options(digits=(x - 1))
p1 <- NULL
for(i in (2^seq(x,y,1)))
{p1 <- rbind(p1,cbind(grep(strtrim(pi,x),(11.71875/2) - exp(1) + (tanh((11.71875/i) * 1:i)) / 2000,value=TRUE)))}
# options(digits = (x + 1))
plot(p1,lwd=1,col="blue",type="l")
abline(h=pi,col="red")
# template function code B:
x <- 8
y <- 16
options(digits=(22))
p1 <- NULL
plot.default(x,y,type="n", xlim=c(0,15000),ylim=c(3.1415,3.1416))
for(i in (2^seq(x,y,1)))
{
lines(
rbind(p1,cbind(grep(3.1415,(11.71875/2) - exp(1) + (tanh((11.71875/i) * 1:i)) / 2000,value=TRUE))),
lwd=1,col=rgb(runif(1),runif(1),runif(1))
)
}
abline(h=pi,lwd=2,col=rgb(1,0,0,.5))
# These functions below (seqpi,seqe,seqpi_e) need three arguments like:
cat('
seqpi(1,1,8)
seqe(1,1,8)
seqpi_e(1,1,8)
seqpi(0,16,16)
seqe(0,16,16)
seqpi_e(0,16,16)
')
seqpi <- function(x,y,z) {
options(digits=(22))
xi <- 11.71875
#trimpi <- 7
p1 <- NULL
plot.default(x,y,type="n", xlim=c(0,z),ylim=c(3.1415,3.1416))
for(i in (2^seq(x,y,1))) {
lines(
rbind(p1,cbind(grep(3.1415,(11.71875/2) - exp(1) + (tanh((11.71875/i) * 1:i)) / 2000,value=TRUE))),
lwd=2,col=rgb(runif(1),runif(1),runif(1)) )}
abline(h=pi,lwd=6,col=rgb(.5,0,0,.25))
}
seqe <- function(x,y,z) {
options(digits=(22))
xi <- 11.71875
p1 <- NULL
e1 <- NULL
plot.default(x,y,type="n", xlim=c(0,z),ylim=c(2.7182,2.7183))
for(i in (2^seq(x,y,1))) {
lines(rbind(e1,cbind(grep(2.7182,(11.71875/2) - pi + (tanh((11.71875/i) * 1:i)) / 2000,value=TRUE))),
lwd=2,col=rgb(runif(1),runif(1),runif(1)))
}
abline(h=exp(1),lwd=6,col=rgb(0,0,.5,.25))
}
seqpi_e <- function(x,y,z) {
options(digits=(22))
xi <- 11.71875
p1 <- NULL
e1 <- NULL
plot.default(x,y,type="n", xlim=c(0,z),ylim=c(2.7182,3.1416))
for(i in (2^seq(x,y,1))) {
lines(rbind(p1,cbind(grep(3.1415,(11.71875/2) - exp(1) + (tanh((11.71875/i) * 1:i)) / 2000,value=TRUE))),
lwd=2,col=rgb(runif(1),runif(1),runif(1)));
lines(rbind(e1,cbind(grep(2.7182,(11.71875/2) - pi + (tanh((11.71875/i) * 1:i)) / 2000,value=TRUE))),
lwd=2,col=rgb(runif(1),runif(1),runif(1)))
}
abline(h=pi,lwd=6,col=rgb(.5,0,0,.25))
abline(h=exp(1),lwd=6,col=rgb(0,0,.5,.25))
}
# plotting these graphs with numbers bigger than 16 requires more CPU and memory than I got...
seqpi(1,1,8)
seqe(1,1,8)
seqpi_e(1,1,8)
seqpi(0,16,16)
seqe(0,16,16)
seqpi_e(0,16,16)
# for(x in (2^seq(0,16,1))) {pi + (min((((11.71875/2) - exp(1)) + (tanh((11.71875/8) * 1:8)) / 2000 ) - pi))}
for(x in (2^seq(0,4,1))) {match(3.141592,((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000 )}
for(x in (2^seq(0,8,1))) {print(grep(3.1415926535,((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000,value=TRUE))}
for(x in (2^seq(0,8,1))) {plot(grep(3.14159,((11.71875/2) - exp(1)) + (tanh((11.71875/x) * 1:x)) / 2000,value=TRUE))}
plot(c2,ylim=c(pi - (pi *.000001),pi + (pi *.000001)),lwd=1,col="blue",type="l");abline(h=pi,col="red")
# A magic series??
l1 <- c(3,5,11,21,41,83,165,331,661,1322,2644,5288,10576,21154)
as.matrix(l1 / (xi * .1))
cat(' [,1]
[1,] 2.5600000000000001
[2,] 4.2666666666666666
[3,] 9.3866666666666667
[4,] 17.9200000000000017
[5,] 34.9866666666666646
[6,] 70.8266666666666680
[7,] 140.8000000000000114
[8,] 282.4533333333333189
[9,] 564.0533333333332848
[10,] 1128.1066666666665697
[11,] 2256.2133333333331393
[12,] 4512.4266666666662786
[13,] 9024.8533333333325572
[14,] 18051.4133333333338669
')
plot(as.matrix(l1 / (xi * .1)))
plot(as.matrix(l1 / (xi * .1)),type="b")
l1 <- c(3,5,11,21,41,83,165,331,661,1322,2644,5288,10576,21154)
mm1 <- data.frame(cbind("Pxi" = l1))
mm1 <- data.frame(mm1,"P2" = 2^(3:16))
for(i in 1:10) {print(((11.71875/2) - pi) + (tanh((11.71875/(mm1$P2[i])) * 1:mm1$P2[i]) / 2000 )[mm1$Pxi[i]])}
F_pi <- function(x) {((11.71875/2) - exp(1)) + (tanh((11.71875/(mm1$P2[x])) * 1:mm1$P2[x]) / 2000 )[mm1$Pxi[x]]}
F_pe <- function(x) {((11.71875/2) - pi) + (tanh((11.71875/(mm1$P2[x])) * 1:mm1$P2[x]) / 2000 )[mm1$Pxi[x]]}
as.matrix(sapply(1:14, F_pi))
cat(' [,1]
[1,] 3.1415930191734072
[2,] 3.1415925126008446
[3,] 3.1415928546514991
[4,] 3.1415927145681901
[5,] 3.1415926227939854
[6,] 3.1415926707780644
[7,] 3.1415926473344875
[8,] 3.1415926591902781
[9,] 3.1415926532962670
[10,] 3.1415926532962670
[11,] 3.1415926532962670
[12,] 3.1415926532962670
[13,] 3.1415926532962670
[14,] 3.1415926536666201
')
as.matrix(mm1$Pxi /mm1$P2) / .1171875
cat(' [,1]
[1,] 3.2000000000000002
[2,] 2.6666666666666665
[3,] 2.9333333333333331
[4,] 2.7999999999999998
[5,] 2.7333333333333334
[6,] 2.7666666666666666
[7,] 2.7500000000000000
[8,] 2.7583333333333333
[9,] 2.7541666666666669
[10,] 2.7541666666666669
[11,] 2.7541666666666669
[12,] 2.7541666666666669
[13,] 2.7541666666666669
[14,] 2.7544270833333333
')
as.matrix(mm1$Pxi /mm1$P2) / 11.71875
cat(' [,1]
[1,] 0.032000000000000001
[2,] 0.026666666666666668
[3,] 0.029333333333333333
[4,] 0.028000000000000001
[5,] 0.027333333333333334
[6,] 0.027666666666666666
[7,] 0.027500000000000000
[8,] 0.027583333333333335
[9,] 0.027541666666666666
[10,] 0.027541666666666666
[11,] 0.027541666666666666
[12,] 0.027541666666666666
[13,] 0.027541666666666666
[14,] 0.027544270833333332
')
# These are some earlier notes in process:
# A reasonable approximation of pi can be made with with aR below
# cR is a correction
((11.71875/2) - exp(1) + exp(-11.71875/(2^14)) / 2000)
format(pi - ((11.71875/2) - exp(1) + exp(-11.71875/(2^14)) / 2000),scientific=FALSE)
#[1] "-0.00000016045116035812157"
aR <- 11.71875/2 - exp(1)
aR
cR <- exp(-11.71875/(2^14)) / 2000
cR
piR <- aR + cR
pi
piR
piR - pi
format(piR - pi, scientific=FALSE)
aR <- 11.71875/2 - exp(1)
aR
#[1] 3.1410931715409549
cR <- exp(-11.71875/(2^14)) / 2000
cR
#[1] 0.00049964249999855234
piR <- aR + cR
pi
#[1] 3.1415926535897931
piR
#[1] 3.1415928140409535
piR - pi
#[1] 1.6045116035812157e-07
format(piR - pi, scientific=FALSE)
#[1] "0.00000016045116035812157"
eR <- ((11.71875/2) - pi + tanh(11.71875) / 2000 )
format(eR - exp(1),scientific=FALSE)
Vx <- as.data.frame(cbind(as.matrix(1 / 1:10),as.matrix(exp(pi^-11.71875) / 1:10)))
Vx$V3 <- with(Vx,cbind(V2 - V1))
Vx$V4 <- with(Vx,cbind(V1 / V2))
Vx$V5 <- with(Vx,cbind(V2 / V1))
Vx
cat(' V1 V2 V3 V4 V5
1 1.00000000000000000 1.00000149288121576 1.4928812157588567e-06 0.99999850712101290 1.0000014928812158
2 0.50000000000000000 0.50000074644060788 7.4644060787942834e-07 0.99999850712101290 1.0000014928812158
3 0.33333333333333331 0.33333383096040525 4.9762707193812261e-07 0.99999850712101290 1.0000014928812158
4 0.25000000000000000 0.25000037322030394 3.7322030393971417e-07 0.99999850712101290 1.0000014928812158
5 0.20000000000000001 0.20000029857624316 2.9857624314622022e-07 0.99999850712101301 1.0000014928812158
6 0.16666666666666666 0.16666691548020263 2.4881353596906131e-07 0.99999850712101290 1.0000014928812158
7 0.14285714285714285 0.14285735612588796 2.1326874510840810e-07 0.99999850712101290 1.0000014928812158
8 0.12500000000000000 0.12500018661015197 1.8661015196985709e-07 0.99999850712101290 1.0000014928812158
9 0.11111111111111110 0.11111127698680175 1.6587569064141494e-07 0.99999850712101290 1.0000014928812158
10 0.10000000000000001 0.10000014928812158 1.4928812157311011e-07 0.99999850712101301 1.0000014928812158
')
# Other Sequences/Series/Prodcuts of interest
23.4375/2
#[1] 11.71875
(23.4375 * 5^1)/1
#[1] 117.1875
(23.4375 * 5^2)/.5
#[1] 1171.875
(23.4375 * 5^3)/.25
#[1] 11718.75
(23.4375 * 5^4)/.125
#[1] 117187.5
(23.4375 * 5^5)/.0625
#[1] 1171875
(23.4375 * 5^6)/.03125
#[1] 11718750
11.71875/2
#[1] 5.859375
(11.71875 * 5^1)/1
#[1] 58.59375
(11.71875 * 5^2)/.5
#[1] 585.9375
(11.71875 * 5^3)/.25
#[1] 5859.375
(11.71875 * 5^4)/.125
#[1] 58593.75
(11.71875 * 5^5)/.0625
#[1] 585937.5
(11.71875 * 5^6)/.03125
#[1] 5859375
exp(pi^-11.71875)
#[1] 1.0000014928812158
exp(pi^-11.71875)/2
#[1] 0.50000074644060788
exp(pi^-11.71875)/3
#[1] 0.33333383096040525
exp(pi^-11.71875)/4
#[1] 0.25000037322030394
exp(pi^-11.71875)/5
#[1] 0.20000029857624316
exp(pi^-11.71875)/6
#[1] 0.16666691548020263
exp(pi^-11.71875)/7
#[1] 0.14285735612588796
exp(pi^-11.71875)/8
#[1] 0.12500018661015197
exp(pi^-11.71875)/9
#[1] 0.11111127698680175
1/9
#[1] 0.1111111111111111
exp(pi^-11.71875)/10
#[1] 0.10000014928812158
p2 <- 2^(1:16)
pX <- data.frame(cbind(p2,(187.5/16) * p2)) # becomes V2
pX$V3 <- cbind((((pi + exp(1)) * 2) * p2))
pX$V4 <- with(pX,cbind(V3 - V2))
pX$V5 <- with(pX,cbind(V4 * 1000))
pX$V6 <- with(pX,cbind(p2/V5))
plot(pX$p2,type="p",col="red");par(new=T);plot(pX$V5,type="l",col="blue");par(new=F);
pX
cat(' p2 V2 V3 V4 V5 V6
1 2 23.4375 23.439497928195351 0.0019979281953510508 1.9979281953510508 1.0010369765308733
2 4 46.8750 46.878995856390702 0.0039958563907021016 3.9958563907021016 1.0010369765308733
3 8 93.7500 93.757991712781404 0.0079917127814042033 7.9917127814042033 1.0010369765308733
4 16 187.5000 187.515983425562808 0.0159834255628084065 15.9834255628084065 1.0010369765308733
5 32 375.0000 375.031966851125617 0.0319668511256168131 31.9668511256168131 1.0010369765308733
6 64 750.0000 750.063933702251234 0.0639337022512336262 63.9337022512336262 1.0010369765308733
7 128 1500.0000 1500.127867404502467 0.1278674045024672523 127.8674045024672523 1.0010369765308733
8 256 3000.0000 3000.255734809004935 0.2557348090049345046 255.7348090049345046 1.0010369765308733
9 512 6000.0000 6000.511469618009869 0.5114696180098690093 511.4696180098690093 1.0010369765308733
10 1024 12000.0000 12001.022939236019738 1.0229392360197380185 1022.9392360197380185 1.0010369765308733
11 2048 24000.0000 24002.045878472039476 2.0458784720394760370 2045.8784720394760370 1.0010369765308733
12 4096 48000.0000 48004.091756944078952 4.0917569440789520741 4091.7569440789520741 1.0010369765308733
13 8192 96000.0000 96008.183513888157904 8.1835138881579041481 8183.5138881579041481 1.0010369765308733
14 16384 192000.0000 192016.367027776315808 16.3670277763158082962 16367.0277763158082962 1.0010369765308733
15 32768 384000.0000 384032.734055552631617 32.7340555526316165924 32734.0555526316165924 1.0010369765308733
16 65536 768000.0000 768065.468111105263233 65.4681111052632331848 65468.1111052632331848 1.0010369765308733
')
t1 <- ((pi + exp(1)) /(187.5/16))^-(1:16)
t2 <- 2^(1:16)
t3 <- t1/t2
s4 <- cbind(t1,t2,t3)
plot(s4[,3],type="l")
s4
cat(' t1 t2 t3
[1,] 1.9998295246594897 2 0.99991476232974486
[2,] 3.9993181276998007 4 0.99982953192495017
[3,] 7.9979544702799732 8 0.99974430878499665
[4,] 15.9945454865482386 16 0.99965909290926491
[5,] 31.9863642975083522 32 0.99957388429713601
[6,] 63.9672757086713943 64 0.99948868294799054
[7,] 127.9236465742348372 128 0.99940348886120967
[8,] 255.8254853212606292 256 0.99931830203617433
[9,] 511.6073587057998679 512 0.99923312247226537
[10,] 1023.1275009729168914 1024 0.99914795016886415
[11,] 2046.0805839367201315 2048 0.99906278512535163
[12,] 4091.8123615891813643 4096 0.99897762734110873
[13,] 8182.9271700727167627 8192 0.99889247681551718
[14,] 16364.4593528497443913 16384 0.99880733354795803
[15,] 32726.1289689190489298 32768 0.99872219753781277
[16,] 65446.6789398585315212 65536 0.99863706878446246
')
t1 <- ((pi + exp(1)) /(187.5/16))^-(1:1023)
t2 <- 2^(1:1023)
t3 <- t1/t2
s4 <- cbind(t1,t2,t3)
plot(s4[,3],type="l")
#print(s4) #
print("Table that shows a relationship of 2^1:1023 produced by the dividend of (pi + e) / xi^-(1:1023):")
cat(' t1 t2 t3
[1,] 1.9998295246594897e+00 2.0000000000000000e+00 0.99991476232974486
[2,] 3.9993181276998007e+00 4.0000000000000000e+00 0.99982953192495017
[3,] 7.9979544702799732e+00 8.0000000000000000e+00 0.99974430878499665
[4,] 1.5994545486548239e+01 1.6000000000000000e+01 0.99965909290926491
[5,] 3.1986364297508352e+01 3.2000000000000000e+01 0.99957388429713601
[6,] 6.3967275708671394e+01 6.4000000000000000e+01 0.99948868294799054
[7,] 1.2792364657423484e+02 1.2800000000000000e+02 0.99940348886120967
[8,] 2.5582548532126063e+02 2.5600000000000000e+02 0.99931830203617433
')
t1 <- ((pi + exp(1)) /(187.5/16))^-(1:1023)
t2 <- 2^(1:1023)
t3 <- t1/t2
s4 <- cbind(t1,t2,t3)
plot(s4[,3])
l <- {};
j <- t3;
for(i in t3) {k <- (j - i);j <- i;l <-rbind(l,k)};
# long time for plot render
#plot(l[-1])
# l
cat(' [,1]
k 0.0000000000000000e+00
k 8.5230404794689463e-05
k 8.5223139953516380e-05
k 8.5215875731736723e-05
k 8.5208612128906402e-05
....
k 7.8152951918264435e-05
k 7.8146290342417934e-05
k 7.8139629334228466e-05
k 7.8132968894695232e-05
k 7.8126309022596985e-05
')
for (i in 1:16) {print(((pi + exp(1)) /(187.5/16))^-i)}
#[1] 1.9998295246594897
#[1] 3.9993181276998007
#[1] 7.9979544702799732
#[1] 15.994545486548239
#[1] 31.986364297508352
#[1] 63.967275708671394
#[1] 127.92364657423484
#[1] 255.82548532126063
#[1] 511.60735870579987
#[1] 1023.1275009729169
#[1] 2046.0805839367201
#[1] 4091.8123615891814
#[1] 8182.9271700727168
#[1] 16364.459352849744
#[1] 32726.128968919049
#[1] 65446.678939858532
as.matrix(2^(1:16))
cat(' [,1]
[1,] 2
[2,] 4
[3,] 8
[4,] 16
[5,] 32
[6,] 64
[7,] 128
[8,] 256
[9,] 512
[10,] 1024
[11,] 2048
[12,] 4096
[13,] 8192
[14,] 16384
[15,] 32768
[16,] 65536
')
t1 <- ((pi + exp(1)) /(187.5/16))^-(1:16)
t2 <- 2^(1:16)
cbind(t1,t2)
cat(' t1 t2
[1,] 1.9998295246594897 2
[2,] 3.9993181276998007 4
[3,] 7.9979544702799732 8
[4,] 15.9945454865482386 16
[5,] 31.9863642975083522 32
[6,] 63.9672757086713943 64
[7,] 127.9236465742348372 128
[8,] 255.8254853212606292 256
[9,] 511.6073587057998679 512
[10,] 1023.1275009729168914 1024
[11,] 2046.0805839367201315 2048
[12,] 4091.8123615891813643 4096
[13,] 8182.9271700727167627 8192
[14,] 16364.4593528497443913 16384
[15,] 32726.1289689190489298 32768
[16,] 65446.6789398585315212 65536
')
cbind(t1,t2,t1/t2,"t2-t1"=t2-t1)
cat(' t1 t2
[1,] 1.9998295246594897 2 0.99991476232974486 1.7047534051028990e-04
[2,] 3.9993181276998007 4 0.99982953192495017 6.8187230019933764e-04
[3,] 7.9979544702799732 8 0.99974430878499665 2.0455297200268063e-03
[4,] 15.9945454865482386 16 0.99965909290926491 5.4545134517614002e-03
[5,] 31.9863642975083522 32 0.99957388429713601 1.3635702491647805e-02
[6,] 63.9672757086713943 64 0.99948868294799054 3.2724291328605659e-02
[7,] 127.9236465742348372 128 0.99940348886120967 7.6353425765162797e-02
[8,] 255.8254853212606292 256 0.99931830203617433 1.7451467873937077e-01
[9,] 511.6073587057998679 512 0.99923312247226537 3.9264129420013205e-01
[10,] 1023.1275009729168914 1024 0.99914795016886415 8.7249902708310856e-01
[11,] 2046.0805839367201315 2048 0.99906278512535163 1.9194160632798685e+00
[12,] 4091.8123615891813643 4096 0.99897762734110873 4.1876384108186357e+00
[13,] 8182.9271700727167627 8192 0.99889247681551718 9.0728299272832373e+00
[14,] 16364.4593528497443913 16384 0.99880733354795803 1.9540647150255609e+01
[15,] 32726.1289689190489298 32768 0.99872219753781277 4.1871031080951070e+01
[16,] 65446.6789398585315212 65536 0.99863706878446246 8.9321060141468479e+01
')
# Some interesting sequences and explorations
(pi + exp(1)) /(187.5/16)
#[1] 0.50004262246816744
# correction factors gives 11.71875 :
exp(1.004057)^pi /2
#[1] 11.718759272520591
exp(pi^1.0035369) /2
#[1] 11.718759367405029
# interesting because give us nearly 1
exp(pi^-11.718759367405029)
#[1] 1.0000014928652075
exp(pi^-11.71875)
#[1] 1.0000014928812158
# power of 2 sequence resulting in almost 1, 1/2, 1/4, 1/8, 1/16
187.5/8 / exp(pi)
#[1] 1.0128262093071623
187.5/16 / exp(pi)
#[1] 0.50641310465358114
187.5/32 / exp(pi)
#[1] 0.25320655232679057
187.5/64 / exp(pi)
#[1] 0.12660327616339528
187.5/128 / exp(pi)
#[1] 0.063301638081697642
# ...
exp(11.718759272520591) /2048
#[1] 59.987744700119116 # 60
exp(11.718759272520591) /1024
#[1] 119.97548940023823 # 120
exp(11.718759272520591) /512
#[1] 239.95097880047646 # 240
exp(11.718759272520591) /256
#[1] 479.90195760095293 #480
exp(pi) / (187.5/2)
#[1] 0.24683405474964551 #.25
exp(pi) / (187.5/4)
#[1] 0.49366810949929102 #.5
exp(pi) / (187.5/8)
#[1] 0.98733621899858204 #1
exp(pi) / (187.5/16)
#[1] 1.9746724379971641 #2
exp(pi) / (187.5/32)
#[1] 3.9493448759943282 #4
exp(pi) / (187.5/64)
#[1] 7.8986897519886563 #8
exp(pi) / (187.5/128)
#[1] 15.797379503977313 #16
exp(pi) / (187.5/256)
#[1] 31.594759007954625 #32
exp(pi) / (187.5/512)
#[1] 63.18951801590925 #64
exp(pi) / (187.5/1024)
#[1] 126.3790360318185 #128
exp(pi) / (187.5/2048)
#[1] 252.758072063637 #256
(187.5/16)
#[1] 11.71875
((pi + exp(1)) /(187.5/16)) * 2
#[1] 1.0000852449363349
exp(((pi + exp(1)) /(187.5/16)) * 2)
#[1] 2.7185135580972011
exp(1)
#[1] 2.7182818284590451
exp(((pi + exp(1)) /(187.5/16)) * 1.999)
#[1] 2.7171545252638842
exp(1)
#[1] 2.7182818284590451
187.5/8 / exp(pi)
#[1] 1.0128262093071623
187.5/16 / exp(pi)
#[1] 0.50641310465358114
exp(pi)
#[1] 23.140692632779267
exp(pi) / 2
#[1] 11.570346316389633
exp(pi) / exp(11.570346316389633)
#[1] 0.0002184936375100267
exp(pi) / exp(187.5/16)
#[1] 0.00018835965834893231
tanh(187.5/16) + (((pi + exp(1)) * 2 ) - 187.5/16)
#[1] 1.0009989639651642
tanh(187.5/16)
#[1] 0.99999999986748866
(((pi + exp(1)) * 2 ) - 187.5/16)
#[1] 0.00099896409767552541
(((pi + exp(1)) * 2 ) - 187.5/16) / tanh(187.5/16)
#[1] 0.00099896409780789956
(((pi + exp(1)) * 2 ) - 187.5/16) - tanh(187.5/16)
#[1] -0.99900103576981314
acosh(187.5/16)
#[1] 3.1525118680526645
asinh(187.5/16)
#[1] 3.1561527971624654
acosh(187.5/16) - pi
#[1] 0.010919214462871363
asinh(187.5/16) - pi
#[1] 0.014560143572672324
asinh(187.5/16) / pi
#[1] 1.004634637643437
(pi + exp(1))
#[1] 5.8598744820488378
(pi + exp(1)) /(187.5/16)
#[1] 0.50004262246816744
(187.5/16) /2
#[1] 5.859375
pi + exp(1)
#[1] 5.8598744820488378
t1 <- (187.5/16) /2
t2 <- pi + exp(1)
t1 - exp(1)
#[1] 3.1410931715409549
# golden ratio stuff
digits(options=22)
gr <- (sqrt(5) + 1) / 2
xi <- (187.5 /16) / 10
XI <- 187.5 /16
pi
#[1] 3.1415926535897931
exp(1)
#[1] 2.7182818284590451
gr
#[1] 1.6180339887498949
xi
#[1] 1.171875
((pi*xi - gr) - (pi/xi - gr))
[1] 1.0007281598622488
gr1 <- ((pi*xi - gr) - (pi/xi - gr))
((pi*xi - gr) / (pi/xi - gr))
#[1] 1.9416032511543122
gr
#[1] 1.6180339887498949
pi
#[1] 3.1415926535897931
xi
#[1] 1.171875
exp(1)
#2.7182818284590451
grxi <- gr + xi
#[1] 2.7899089887498949
exp(1) - (gr + xi)
#[1] -0.071627160290849812
(gr + xi) - exp(1)
#[1] 0.071627160290849812
((xi - exp(1)) + gr)
#[1] 0.071627160290849812
grxi <- gr + xi
format(grxi - exp(1), scientific=FALSE)
# [1] "0.071627160290849812"
exp(1) - ((gr + xi) - (xi/gr))
# [1] 0.652631420275433
exp(1) - ((gr + xi) - ((xi/gr) / 10))
#[1] 0.00079869776577856655
exp(1)
#[1] 2.7182818284590451
((gr + xi) - ((xi/gr) / 10))
#[1] 2.7174831306932665
#[1] 2.7063293868263711
exp(1)
#[1] 2.7182818284590451
sqrt(exp(1))
[1] 1.6487212707001282
gr
#[1] 1.6180339887498949
xi/sqrt(3/16)/sqrt(exp(1))
#[1] 1.6414717483914856
gr/exp(1)
#[1] 0.59524143957771114
exp(1)/gr
#[1] 1.679990560988901
gr
#[1] 1.6180339887498949
sqrt(exp(1)/(xi/sqrt(3/16)))
#[1] 1.0022058057366372
(xi/sqrt(3/16))/sqrt(exp(1))
#[1] 1.6414717483914856
-digamma(1) @ R 'built-in' Euler–Mascheroni constant approximation; Good for 15 digits
#[1] 0.57721566490153231
exp(1)/gr - ((gr/exp(1))/10)
#[1] 1.6204664170311298
gr
#[1] 1.6180339887498949
gr/exp(1) + -digamma(1)
#[1] 1.1724571044792436
xi
#[1] 1.171875
options(digits=15)
xi
#[1] 1.171875
gr/exp(1) + -digamma(1)
#[1] 1.1724571044792
xi - (gr/exp(1) + -digamma(1))
#[1] -0.00058210447924356
sxi <- (gr/exp(1) + -digamma(1)) * 10
pR <- ((sxi/2) - exp(1) + (-digamma(1))/1000 )
format(pR - pi,scientific=FALSE)
#[1] "0.002988256012280921"
pR <- ((Xi/2) - exp(1) + tanh(Xi) / 2000 )
format(pR - pi,scientific=FALSE)
# [1] "0.0000005179510953468025"
No comments:
Post a Comment