Keyword: 多ガンマ関数, ポリガンマ関数, 実数
概要
本サンプルは実数xの多ガンマ関数（プサイ関数の導関数）を求めるC#によるサンプルプログラムです。 本サンプルは引数xを読み込み、xの各値について以下に示される多ガンマ関数を求めて出力します。
※本サンプルはnAG Library for .NETに含まれる関数 s14ae() のExampleコードです。本サンプル及び関数の詳細情報は s14ae のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
（本関数の詳細はs14ae のマニュアルページを参照）- １行目はタイトル行で読み飛ばされます。
- ２行目に引数xの値と導関数の階数kを指定しています。
出力結果
（本関数の詳細はs14ae のマニュアルページを参照）- ５行目に引数xの値、導関数の階数、多ガンマ関数の値とエラーコードが出力されています。エラーコード"0"はエラーがなかったことを意味します。
ソースコード
（本関数の詳細はs14ae のマニュアルページを参照）
※本サンプルソースコードは .NET環境用の科学技術・統計計算ライブラリである「nAG Library for .NET」の関数を呼び出します。
サンプルのコンパイル及び実行方法
| このソースコードをダウンロード |
// s14ae Example Program Text
// C# version, nAG Copyright 2008
using System;
using NagLibrary;
namespace NagDotNetExamples
{
public class S14AEE
{
static bool defaultdata = true;
static string datafile = "";
static void Main(String[] args)
{
if (args.Length == 1)
{
defaultdata = false;
datafile = args[0];
}
StartExample();
}
public static void StartExample()
{
try
{
PrintManager.Warning = new PrintManager.MessageLogger(discardmessage);
DataReader sr = null;
if (defaultdata)
{
sr = new DataReader("exampledata/s14aee.d");
}
else
{
sr = new DataReader(datafile);
}
double x, y; int k;
int ifail;
Console.WriteLine("s14ae Example Program Results");
// Skip heading in data file
sr.Reset();
Console.WriteLine("");
Console.WriteLine(" {0}", " x k (d^k/dx^k)psi(x) ifail");
Console.WriteLine("");
try
{
while (true)
{
sr.Reset();
x = double.Parse(sr.Next());
k = int.Parse(sr.Next());
//
y = S.s14ae(x, k, out ifail);
//
if (ifail >= 0)
{
Console.WriteLine(" {0,5:f1}{1,5} {2,12:e4} {3,7}", x, k, y, ifail);
}
else
{
Console.WriteLine("** s14ae failed with ifail = {0,5}", ifail);
}
}
}
catch
{
}
//
}
catch (Exception e)
{
Console.WriteLine(e.Message);
Console.WriteLine("Exception Raised");
}
}
static void discardmessage(String message)
{
}
}
}