perl

perlで正規分布

ke-sato 2017年2月2日 / 2020年12月15日

乱数の出力はボックスミューラー法。

#!/usr/bin/perl

$Pi = 3.14159265359;

sub RandBoxMuller {
    my ($u, $sigma) = @_;
    my ($r1, $r2) = (rand(), rand());

    while ($r1 == 0) { $r1 = rand(); }

    return ($sigma * sqrt(-2 * log($r1)) * sin(2 * $Pi * $r2)) + $u;
}

#!/usr/bin/perl

$Pi = 3.14159265359;

sub RandBoxMuller {

my ($u, $sigma) = @_;

my ($r1, $r2) = (rand(), rand());

while ($r1 == 0) { $r1 = rand(); }

return ($sigma * sqrt(-2 * log($r1)) * sin(2 * $Pi * $r2)) + $u;

}

C言語で簡易メモリ管理

errno一覧

Share LaTeXで数学の図形・グラフを作成

gnuplotで正規分布を表示

ソースコードと完全に同じ文字列を出力する