Power Law:(17)Graph with Lucas number

Power Law:(11)のデータに、リュカ数の基準線を引いてみる。
Rのabline(a, b)を使用する。
切片(A=0)にリュカ数の何番目かをセットし、傾きをexp(5.008) -> exp(5)とする。
The base point lines of the number of Lucas are pulled to the data of Power Law:(11).
abline(a, b) of R is used.
The number of Lucas is set in cut (A=0), and the inclination is assumed the exp(5.008) -> exp(5).
~~~
[1]リュカ数の場合,in case of Lucas number
data_ab <- read.csv("powerlaw-11-a-b.txt");
plot(log(abs(data_ab$A)), log(data_ab$B), xlab="log(abs(a))", ylab="log(b)", xlim=c(-2, 2), ylim=c(3, 15), col="orange", pch=20, main="Scatter chart:a-b, Power Law(11) with Lucas")

Lucas:n=5
abline(log(11), log(exp(5)), col="gray", lty=2)

Lucas:n=6
abline(log(18), log(exp(5)), col="gray", lty=2)

Lucas:n=7
abline(log(29), log(exp(5)), col="gray", lty=2)

Lucas:n=8
abline(log(47), log(exp(5)), col="gray", lty=2)

Lucas:n=9
abline(log(76), log(exp(5)), col="gray", lty=2)

Lucas:n=10
abline(log(123), log(exp(5)), col="gray", lty=2)

~~~
Lucas:n=12
abline(log(322), log(exp(5)), col="gray", lty=2)

Lucas:n=14
abline(log(843), log(exp(5)), col="gray", lty=2)

Lucas:n=18
abline(log(5778), log(exp(5)), col="gray", lty=2)

Lucas:n=19
abline(log(9349), log(exp(5)), col="gray", lty=2)

Lucas:n=21
abline(log(24476), log(exp(5)), col="gray", lty=2)

Lucas:n=22
abline(log(39603), log(exp(5)), col="gray", lty=2)

Lucas:n=25
abline(log(167761), log(exp(5)), col="gray", lty=2)

Lucas:n=26
abline(log(271443), log(exp(5)), col="gray", lty=2)

abline(v = log(1), col="red")


~~~
[2]フィボナッチ数の場合,in case of Fibonacci number
Power Law:(15)では、グラフから読み取った数字からリュカ数が近いと見ているが、
フィボナッチのグラフではどうなっているかをみる
It sees how to become it in Fibonacci's graph though it is thought that the
number of Lucas is near from the figure read from the graph in Power Law:(15).

data_ab <- read.csv("powerlaw-11-a-b.txt");
plot(log(abs(data_ab$A)), log(data_ab$B), xlab="log(abs(a))", ylab="log(b)", xlim=c(-2, 2), ylim=c(3, 15), col="orange", pch=20, main="Scatter chart:a-b, Power Law(11) with Fibonacci")

Fibonacci:n=7
abline(log(13), log(exp(5)), col="gray", lty=2)

Fibonacci:n=8
abline(log(21), log(exp(5)), col="gray", lty=2)

Fibonacci:n=9
abline(log(34), log(exp(5)), col="gray", lty=2)

Fibonacci:n=10
abline(log(55), log(exp(5)), col="gray", lty=2)

Fibonacci:n=11
abline(log(89), log(exp(5)), col="gray", lty=2)

Fibonacci:n=12 :???
abline(log(144), log(exp(5)), col="gray", lty=2)

Fibonacci:n=14
abline(log(377), log(exp(5)), col="gray", lty=2)

Fibonacci:n=15
abline(log(610), log(exp(5)), col="gray", lty=2)

~~~
Fibonacci:n=20
abline(log(6765), log(exp(5)), col="gray", lty=2)

Fibonacci:n=21
abline(log(10946), log(exp(5)), col="gray", lty=2)

Fibonacci:n=22
abline(log(17711), log(exp(5)), col="gray", lty=2)

Fibonacci:n=23
abline(log(28657), log(exp(5)), col="gray", lty=2)

Fibonacci:n=24
abline(log(46368), log(exp(5)), col="gray", lty=2)

Fibonacci:n=26
abline(log(121393), log(exp(5)), col="gray", lty=2)

Fibonacci:n=27
abline(log(196418), log(exp(5)), col="gray", lty=2)

abline(v = log(1), col="red")


~~~
end