在螺旋上的質數

下圖是把1到10000的質數畫在阿基米德螺旋線上


Ｒ的程式碼參考《Plotting Times of Discrete Events》，依序執行

• library(numbers)
• library(ggplot2)
• 到上面的連結載入函式 spiral.length
• 到上面的連結載入函式 spiral.plot
• spiral.plot(Primes(1,10000), 20)

觀察上面黏在一起的點，就是所謂的孿生質數（twin prime）

2017-Jan-19更新：

找到一篇文章《Prime and Random 'Nixon' Spirals》，動畫很有趣
_EOF_