Incanter

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に指が慣れてきたせいで、括弧の入力ミスが極端に減って快適。