FC2ブログ

# ルジャンドル多項式 Legendre Polynomial
# 引数 次数 値 ($N, $X)
# 戻り値 ルジャンドル多項式 ($LegendrePolynomial)
sub LEGENDREPOLYNOMIAL{
my ($N, $X) = @_;
my $LegendrePolynomial = 0;
my $NextLegendre = 0;
my $Legendre = $X;
my $PrevLegendre = 1;

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

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

# 計算
for(my $i = 1; $i < $N; $i++){
$NextLegendre = ((((2 * $i) + 1) * $X * $Legendre) - ($i * $PrevLegendre)) / ($i + 1);

# 漸化式
$PrevLegendre = $Legendre;
$Legendre = $NextLegendre;
}

# ルジャンドル多項式 Legendre Polynomial
$LegendrePolynomial = $Legendre;

return $LegendrePolynomial;
}


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

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

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

検索フォーム


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