Keyword: プサイ関数, 対数関数
概要
本サンプルはプサイ関数と対数関数との差分（Ψ(x)−ln x）を求めるC#によるサンプルプログラムです。 本サンプルは引数xを読み込み、xの各値についてΨ(x)−ln xを求めて出力します。この場合Ψは以下で表わされます。
※本サンプルはnAG Library for .NETに含まれる関数 s14ac() のExampleコードです。本サンプル及び関数の詳細情報は s14ac のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
（本関数の詳細はs14ac のマニュアルページを参照）- １行目はタイトル行で読み飛ばされます。
- ２〜５行目に引数xの値を指定しています。
出力結果
（本関数の詳細はs14ac のマニュアルページを参照）| この出力例をダウンロード |
s14ac Example Program Results
x psi(x)-log(x) ifail
0.1000 -8.1212 0
0.5000 -1.2704 0
3.6000 -0.1453 0
8.0000 -0.0638 0
- ５〜８行目に引数xの値、Ψ(x)−ln xの値とエラーコードが出力されています。エラーコード"0"はエラーがなかったことを意味します。
ソースコード
（本関数の詳細はs14ac のマニュアルページを参照）
※本サンプルソースコードは .NET環境用の科学技術・統計計算ライブラリである「nAG Library for .NET」の関数を呼び出します。
サンプルのコンパイル及び実行方法
| このソースコードをダウンロード |
// s14ac Example Program Text
// C# version, nAG Copyright 2008
using System;
using NagLibrary;
namespace NagDotNetExamples
{
public class S14ACE
{
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/s14ace.d");
}
else
{
sr = new DataReader(datafile);
}
double f, x; int ifail;
Console.WriteLine("s14ac Example Program Results");
// Skip heading in data file
sr.Reset();
Console.WriteLine("");
Console.WriteLine(" {0}", " x psi(x)-log(x) ifail");
Console.WriteLine("");
try
{
while (true)
{
sr.Reset();
x = double.Parse(sr.Next());
//
f = S.s14ac(x, out ifail);
//
if (ifail >= 0)
{
Console.WriteLine(" {0,12:f4}{1,15:f4}{2,7}", x, f, ifail);
}
else
{
Console.WriteLine("** s14ac failed with ifail = {0,5}", ifail);
}
}
}
catch
{
}
//
}
catch (Exception e)
{
Console.WriteLine(e.Message);
Console.WriteLine("Exception Raised");
}
}
static void discardmessage(String message)
{
}
}
}