MATLABよりも簡単にできるんじゃないか、と思って、
Interface Oct,2007
の記事のコードを変換。
(use '(incanter core stats charts)) ; 平均μと分散σ^2を指定して、正規分布から標本を抽出する。 (defn randncq [mu variance n-samples] (let [sigma (sqrt variance)] (for [s (sample-normal n-samples)] (+ (* sigma s) mu)))) ; ヒストグラム表示 (view histogram (randncq 1 4 1000))
Emacsのpareditに指が慣れてきたせいで、括弧の入力ミスが極端に減って快適。