FC2ブログ

# ベッセル多項式 Bessel Polynomial
# 引数 次数 値 ($N, $X)
# 戻り値 ベッセル多項式 ($BesselPolynomial)
sub BESSELPOLYNOMIAL{
my ($N, $X) = @_;
my $BesselPolynomial = 0;
my $Bessel = 0;
my $PrevBessel1 = $X + 1;
my $PrevBessel2 = 1;

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

if($N == 0){
return 1;
}elsif ($N == 1){
return ($X + 1);
}

# 計算
for(my $i = 2; $i <= $N; $i++){
$Bessel = (((2 * $i) - 1) * $X * $PrevBessel1) + $PrevBessel2;

# 漸化式
$PrevBessel2 = $PrevBessel1;
$PrevBessel1 = $Bessel;
}

# ベッセル多項式 Bessel Polynomial
$BesselPolynomial = $Bessel;

return $BesselPolynomial;
}


参考URL
Bessel polynomials - Wikipedia, the free encyclopedia
perlで逆ベッセル多項式 Reverse Bessel Polynomial: perl: 一寸先は闇
オンライン コンパイラ/インタプリタ
テクニカル分析
プロフィール

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

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

検索フォーム


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