« 山崎駅でオーバーラン/いつまで復旧にかかるのか。 | トップページ | 東寺に行って来た。 »

2010年2月15日 (月)

正規乱数をkeisan.casio.jpにUP!

最近確率の計算を続けているのでおまけとしてこんなのつくってみた。

正規乱数

線形合同法+Box-Muller法を使っているだけ。

プログラムはこれ↓

av = 平均;
sig = 標準偏差;
n = 乱数の個数;
seed=乱数の種;
for (j = 1 ; j <= n; j=j+1) {
    seed=mod(1566083941*seed+1,2^32);
    r=seed / 2^32;
    t = sqrt(-2.*ln(1.-r));
    seed=mod(1566083941*seed+1,2^32);
    r=seed / 2^32;
    u = 2.*pi*r;
    println(j, t*cos(u)*sig + av, t*sin(u)*sig + av);
}

まあ乱数のアルゴリズムがこれなんで、実用性はないかもしれないが。

« 山崎駅でオーバーラン/いつまで復旧にかかるのか。 | トップページ | 東寺に行って来た。 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/47569188

この記事へのトラックバック一覧です: 正規乱数をkeisan.casio.jpにUP!:

« 山崎駅でオーバーラン/いつまで復旧にかかるのか。 | トップページ | 東寺に行って来た。 »

最近のコメント

2019年1月
    1 2 3 4 5
6 7 8 9 10 11 12
13 14 15 16 17 18 19
20 21 22 23 24 25 26
27 28 29 30 31    
フォト
無料ブログはココログ