正規乱数を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);
}

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