FC2ブログ

# 逆双曲線余割(近似) Inverse Hyperbolic Cosecant Approximation
# 引数 値 ($X)
# 戻り値 逆双曲線余割 ($InverseHyperbolicCosecantApproximation)
sub INVERSEHYPERBOLICCOSECANT{
my ($X) = @_;
my $InverseHyperbolicCosecantApproximation = 0;
my $Sum = 0;
my $Num = 0;
my $Den1 = 0;
my $Den2 = 1;
my $Den3 = 0;
my $Sign = 0;
my $NumFactorial = 1;
my $DenFactorial = 1;
my $Expansion = 20;

# 値の確認
if(abs($X) <= 1){
return 0;
}

for(my $i = 0; $i <= $Expansion; $i++){
# 符号
$Sign = (($i % 2) == 0 ? 1: -1);
# 分子
$Num = $X ** (-((2 * $i) + 1));
# 分母
$Den1 = 2 ** (2 * $i);
$Den3 = (2 * $i) + 1;

if($i != 0){
# 分母 階乗
for(my $j = (2 * $i); $j > (2 * ($i - 1)); $j--){
$NumFactorial = $NumFactorial * $j;
}

# 分母
$DenFactorial = $DenFactorial * $i;

$Den2 = $DenFactorial * $DenFactorial;
}

$Sum += ($Sign * ($NumFactorial / ($Den1 * $Den2)) * ($Num / $Den3));
}

# 逆双曲線余割(近似) Inverse Hyperbolic Cosecant
$InverseHyperbolicCosecantApproximation = $Sum;

return $InverseHyperbolicCosecantApproximation;
}


参考URL
Inverse hyperbolic function - Wikipedia, the free encyclopedia
Inverse Hyperbolic Cosecant -- from Wolfram MathWorld
オンライン コンパイラ/インタプリタ
テクニカル分析
プロフィール

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

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

検索フォーム


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