最新記事
- perl で特殊関数 40 対数積分 li(x) (Logarithmic Integral li: LOGARITHMICINTEGRALLI)
- perl で特殊関数 39 指数積分 Ei(x) (Exponential Integral Ei: EXPONENTIALINTEGRALEI)
- perl で特殊関数 38 指数積分 En(x) (Exponential Integral En: EXPONENTIALINTEGRALEN)
- perlで統計 34 ロジスティック分布 (逆関数) (Logistic Distribution (Inverse Function): LOGISTICDISTRIBUTIONINVERSE)
- perlで統計 33 コーシー分布 (逆関数) (Cauchy Distribution (Inverse Function): CAUCHYDISTRIBUTIONINVERSE)
- 全ての記事タイトル一覧
- ランキングに参加しています。この記事が参考になりましたら、応援お願い致します。
# 変形シュトルーベ関数 Modified Struve Function
# 引数 値 値 ($N, $X)
# 戻り値 変形シュトルーベ関数 ($ModifiedStruveFunction)
sub MODIFIEDSTRUVEFUNCTION{
my ($N, $X) = @_;
my $ModifiedStruveFunction = 0;
my $ModifiedStruve = 0;
my $PrevModifiedStruve = 0;
my $Num = 0;
my $Den1 = 0;
my $Den2 = 0;
my $Limit = 100;
my $Epsilon = 1.0e-10;
# 計算
for(my $i = 0; $i < $Limit; $i++){
# 分子
$Num = (($X / 2) ** ((2 * $i) + $N + 1));
# 分母
$Den1 = &GAMMAFUNCTION($i + (3 / 2));
$Den2 = &GAMMAFUNCTION($i + $N + (3 / 2));
# 一つ前
$PrevModifiedStruve = $ModifiedStruve;
# 変形シュトルーベ関数 Modified Struve Function
$ModifiedStruve += ($Num / ($Den1 * $Den2));
# 収束判定
last if(abs($ModifiedStruve - $PrevModifiedStruve) < $Epsilon);
}
# 変形シュトルーベ関数 Modified Struve Function]
$ModifiedStruveFunction = $ModifiedStruve;
return $ModifiedStruveFunction;
}
# ガンマ関数(近似) Gamma Function Approximation
# 引数 値 ($X)
# 戻り値 ガンマ関数(近似) ($GammaFunction)
sub GAMMAFUNCTION{
my ($X) = @_;
my $GammaFunction = 0;
my $Factorial = 1;
my $Diff = abs($X - int($X));
my $Pi = atan2(1, 1) * 4;
my $Temp = 0;
# 値の確認
if(($X <= 0) && ($Diff == 0)){
return 0;
}
if($Diff == 0){
for(my $i = $X - 1; $i >= 2; $i--){
$Factorial *= $i;
}
# X-1の階乗 ガンマ関数 Gamma Function
$GammaFunction = $Factorial;
}elsif($Diff == 0.5){
my $Start = ($X > 0 ? ((2 * int($X)) - 1) : ((2 * int(abs($X) + 1)) - 1));
for(my $i = $Start; $i > 1; $i -= 2){
$Factorial *= $i;
}
if($X > 0){
$Temp = ($Factorial / (2 ** int($X))) * sqrt($Pi);
}else {
$Temp = ((((-1) ** int(abs($X) + 1)) * (2 ** int(abs($X) + 1))) / $Factorial) * sqrt($Pi);
}
# 半整数 ガンマ関数 Gamma Function
$GammaFunction = $Temp
}else {
if($X > 0){
$Temp = exp(&LOGGAMMAFUNCTION($X));
}else {
$Temp = $Pi / (sin($Pi * $X) * exp(&LOGGAMMAFUNCTION(1 - $X)));
}
# ガンマ関数(近似) Gamma Function Approximation
$GammaFunction = $Temp
}
return $GammaFunction;
}
# ログガンマ関数(近似) Log Gamma Function Approximation
# 引数 値 ($x)
# 戻り値 ログガンマ関数(近似) ($LogGammaFunction)
sub LOGGAMMAFUNCTION{
my ($x) = @_;
my $LogGammaFunction = 0;
my @Bernoulli = ((1 / 12), (1 / 360), (1 / 1260), (1 / 1680), (1 / 1188), (691 / 360360), (1 / 156), (3617 / 122400), (43867 / 244188), (174611 / 125400), (77683 / 5796), (236364091 / 1506960));
my $Pi = atan2(1, 1) * 4;
my $X = $x;
my $Sum = 0;
my $Sign = 0;
my $Power1 = 1;
my $Power2 = 1;
my $Count = @Bernoulli - 1;
for($X = $x; $X <= $Count; $X++){
$Power1 *= $X;
}
$Power2 = 1 / ($X * $X);
for(my $i = $Count; $i >= 0; $i--){
# 符号
$Sign = (($i % 2) == 0 ? 1: -1);
$Sum = $Sum + ($Sign * $Bernoulli[$i]);
$Sum = $Sum * $Power2 if($i != 0);
}
$Sum = $Sum / $X;
# ログガンマ関数(近似) Log Gamma Function Approximation
$LogGammaFunction = ((1 / 2) * log(2 * $Pi)) - log($Power1) - $X + (($X - (1 / 2)) * log($X)) + $Sum;
return $LogGammaFunction;
}
参考URL
Modified Struve Functions
変形シュトルーベ関数 - 高精度計算サイト
スポンサードリンク :
Tags :
Home :
シリーズ
月別アーカイブ カテゴリ
オンライン コンパイラ/インタプリタ
Perl リファレンス
株式
テクニカル分析
計算式ライブラリ
プロフィール
Author:雨宮
Firefoxを使用しているので気づかなかったけど、IE6でソースコードを上手くコピーできない
5/3
携帯用ならIE6でもソースコードをコピーできる
携帯用