Гамма-распределение на Перле

Маленькая perl-программа, возвращающая случайное число по Гамма-распределению.

#!/usr/bin/perl
# Генератор случайных чисел по закону гамма-распределения

$k = 20;
$a = 10;

$tmp = 1;

for ($i = 0; $i < $k; $i++){
$tmp *= rand[- $tmp * $a,$tmp * $a];
}

$gamma = - log($tmp) / $a;

print $gamma;

Поблагодарить автора

Оставить комментарий

Your email address will not be published. Required fields are marked *

*

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>