- 126 名前：nobodyさん mailto:sage [2009/11/26(木) 10:05:40 ID:???]
- >>125
ニートのオイラがPerlでその関数を作ってみた。
引数はNORMDISTと似せているけど、最後の引数は精度。
積分は以下サイトの半無限積分を参考
www.cis.twcu.ac.jp/~osada/rikei/rikei2008-8.pdf
my $probability = normdist( x => 160, mean => 170, sigma => 10, precision => 7 );
print $probability;
sub normdist {
my %param = @_;
use Math::Trig qw( pi );
my $h = 1/1_000;
my $precision = 10**-$param{precision};
my $z = ( $param{x} - $param{mean} ) / $param{sigma};
my $const = 1 / sqrt( 2*pi() );
return 1 - calc_integral( $z, $h, $precision, sub { $const * exp( -0.5 * $_[0]**2 ) } );
}
|
 |
