rubyで安定分布乱数の発生法 - arupaka-

rubyでの安定分布乱数の発生法

これだけででるのか？？
でるとしたら、ものすごく楽。
http://rb-gsl.rubyforge.org/randist.html

最低限

require 'rubygems'
require 'narray'
require 'gsl'

#乱数の初期化
#メルセンヌツイスタ、種は0　
 rng=GSL::Rng.alloc

#安定分布を乱数を1000個発生
data2=rng.levy_skew(c=1,alpha=1,beta=1,10000)

グラフで確認

require 'rubygems'
require 'narray'
 require 'gsl'

#乱数の初期化
#メルセンヌツイスタ、種は0　
 rng=GSL::Rng.alloc

#安定分布を発生
data2=rng.levy_skew(c=1,alpha=1,beta=1,10000)


#ヒストグラムを作成
#箱の数50, -30から30の範囲
 h=GSL::Histogram.alloc(50,[-30,30])
 d= h.increment(data2)

#グラフを出力　xterm
 h.graph()