N($mean, $deviation^2)なる正規分布に従った乱数生成。一様乱数はメルセンヌツイスタ。 任意の確率分布の乱数をしたい場合は、、、確か逆関数法とか棄却法とかで頑張るんだった気が。 use strict; use warnings; use Math::Random::MT qw(rand srand); use Math::Trig; srand(1); my $mean = 10.0; my $deviation = 1.0; my $numData = 10000; my @data; for(my $i = 0; $i < $numData; $i++){ my $rnd = normalRnd($mean, $deviation); push(@data, $rnd); print "$rnd\n"; } sub normalRnd{ my ($mean, $dev