FC2ブログ

# ラゲール多項式 Laguerre Polynomial
# 引数 次数 値 ($N, $X)
# 戻り値 ラゲール多項式 ($LaguerrePolynomial)
sub LAGUERREPOLYNOMIAL{
my ($N, $X) = @_;
my $LaguerrePolynomial = 0;
my $NextLaguerre = 0;
my $Laguerre = 1 - $X;
my $PrevLaguerre = 1;

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

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

# 計算
for(my $i = 1; $i < $N; $i++){
$NextLaguerre = ((((2 * $i) + 1 - $X) * $Laguerre) - ($i * $PrevLaguerre)) / ($i + 1);

# 漸化式
$PrevLaguerre = $Laguerre;
$Laguerre = $NextLaguerre;
}

# ラゲール多項式 Laguerre Polynomial
$LaguerrePolynomial = $Laguerre;

return $LaguerrePolynomial;
}


参考URL
ラゲール多項式(グラフ) - 高精度計算サイト
オンライン コンパイラ/インタプリタ
テクニカル分析
プロフィール

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

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

検索フォーム


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