FC2ブログ

# 逆ベッセル多項式 Reverse Bessel Polynomial
# 引数 次数 値 ($N, $X)
# 戻り値 逆ベッセル多項式 ($ReverseBesselPolynomial)
sub REVERSEPOLYNOMIAL{
my ($N, $X) = @_;
my $ReverseBesselPolynomial = 0;
my $ReverseBessel = 0;
my $PrevReverseBessel1 = $X + 1;
my $PrevReverseBessel2 = 1;

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

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

# 計算
for(my $i = 2; $i <= $N; $i++){
$ReverseBessel = (((2 * $i) - 1) * $PrevReverseBessel1) + (($X * $X) * $PrevReverseBessel2);

# 漸化式
$PrevReverseBessel2 = $PrevReverseBessel1;
$PrevReverseBessel1 = $ReverseBessel;
}

# 逆ベッセル多項式 ReverseBessel Polynomial
$ReverseBesselPolynomial = $ReverseBessel;

return $ReverseBesselPolynomial;
}


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

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

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

検索フォーム


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