長い R プログラムも見たが,簡単に描くことができる。
phi <- (1+sqrt(5))/2
radian <- 90:180*pi/180
xy <- rbind(c(cos(radian), -1, 0, 0, -1), c(sin(radian), 1, 1, 0, 0))
par(mar=c(0, 0, 0, 0))
plot(c(-1, phi-1), c(0, 1), type="n", axes=FALSE, xlab="", ylab="", asp=1)
for (m in 1:11) {
lines(t(xy))
xy <- (matrix(c(0, -1, 1, 0), 2)%*%xy+c(0, 1/phi))/phi
}
実行結果は以下の通り
phi <- (1+sqrt(5))/2
radian <- 90:180*pi/180
xy <- rbind(c(cos(radian), -1, 0, 0, -1), c(sin(radian), 1, 1, 0, 0))
par(mar=c(0, 0, 0, 0))
plot(c(-1, phi-1), c(0, 1), type="n", axes=FALSE, xlab="", ylab="", asp=1)
for (m in 1:11) {
lines(t(xy))
xy <- (matrix(c(0, -1, 1, 0), 2)%*%xy+c(0, 1/phi))/phi
}
実行結果は以下の通り
![](https://blogimg.goo.ne.jp/user_image/38/82/b177c967ac6479a7eaee48db9073cf40.png)
※コメント投稿者のブログIDはブログ作成者のみに通知されます