# 黄金比 psi <- (1+sqrt(5))/2 # 黄金角 psi.angle <- (2*pi)/(1+psi) # シリアルナンバー t <- 0:5 # シリアル角 angles <- psi.angle * t # 角は y = ax 的に増える plot(t,angles) # 角に対応する単位円周上の座標 x <- cos(angles) y <- sin(angles) # …

で、に黄金角に近い値を用いて、n=1,2,...で打点する。 その点分布に対してボロノイ図を描く n <- 200 # 点の数 x <- (1+sqrt(5))/2 # 黄金比 delta <- 0.1 # 黄金比からちょっとずらす phi <- 2*pi * (x + delta) # 黄金角 x <- sqrt(1:n) * cos((1:n) * ph…