FC2ブログ

# ベル多項式 Bell Polynomial
# 引数 整数 値 ($n, $X)
# 戻り値 ベル多項式 ($BellPolynomial)
sub BELLPOLYNOMIALl{
my ($n, $X) = @_;
my $BellPolynomial = 0;
my @Bell = (1);
my $Combination = 0;
my $BinomialCoefficient = 0;
my $Sum = 0;
my $N = int($n);

# 整数の確認
if($N < 0){
return 0;
}

# 計算
for(my $i = 1; $i <= $N; $i++){
$Combination = ($i - 1);
$BinomialCoefficient = 1;
$Sum = ($BinomialCoefficient * $Bell[0] * $X);

for(my $j = 1; $j <= ($i - 1); $j++){
# 二項係数 Binomial Coefficient
$Combination = $Combination / $j;
$BinomialCoefficient = $Combination;
$Combination = $Combination * (($i - 1) - $j);

$Sum += $BinomialCoefficient * $Bell[$j] * $X;
}

$Bell[$i] = $Sum;
}

# ベルヌーイ多項式 Bernoulli Polynomial
$BellPolynomial = $Bell[$N];

return $BellPolynomial;
}


参考URL
Bell Polynomial -- from Wolfram MathWorld
オンライン コンパイラ/インタプリタ
テクニカル分析
プロフィール

Author:雨宮
Firefoxを使用しているので気づかなかったけど、IE6でソースコードを上手くコピーできない

5/3
携帯用ならIE6でもソースコードをコピーできる
携帯用

検索フォーム


あわせて読みたいブログパーツ
一寸先は闇 RSS