ボックス＝ミューラー法

よしコクのある乱数の実際のサンプルをアップした。こうやって並べてみると、味の違いは一目瞭然。 https://t.co/esbYaYWYLh

— 深津 貴之 (@fladdict) November 4, 2016

ボックスミューラー法の養殖された乱数とはｗｗ？

10000.times {
  x = rand()
  y = rand()

  z1 = Math.sqrt(-2 * Math.log(x)) * Math.cos(2 * Math::PI * y)
  z2 = Math.sqrt(-2 * Math.log(x)) * Math.sin(2 * Math::PI * y)

  puts "#{z1}\t#{z2}"
}

ruby bm.rb | pbcopy 

googleスプレッドシートに貼り付けてヒストグラムのグラフを作る（重い…）

rand - Google スプレッドシート

1万回程度ではあんまり綺麗なグラフにならないすね

---

おまけ

• Ruby gem → GitHub - s-osa/random_bell
• Javaクラスライブラリ → Random (Java Platform SE 8 )
• npm → gaussian