関連情報
ホーム > 製品 > NAG数値計算ライブラリ > C#向けNAGライブラリ > サンプルソースコード集 > ガンマ分布から疑似乱数ベクトルを生成

C#による ガンマ分布から疑似乱数ベクトルを生成

C#によるサンプルソースコード
使用関数名:g05sj

Keyword: ガンマ分布, 疑似乱数ベクトル

概要

本サンプルはガンマ分布から疑似乱数ベクトルの生成を行うC#によるサンプルプログラムです。 本サンプルは以下に示される確率密度関数をもつガンマ分布(a=5.0、b=1.0)から5個の疑似乱数を生成し出力します。

ガンマ分布のデータ 

※本サンプルはNAG Library for .NETに含まれる関数 g05sj() のExampleコードです。本サンプル及び関数の詳細情報は g05sj のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで

出力結果

(本関数の詳細はg05sj のマニュアルページを参照)

この出力例をダウンロード
g05sj Example Program Results

   5.0702
   6.1337
   3.1018
   3.9863
   4.9648


  • 3〜7行目に生成された疑似乱数が出力されています。

ソースコード

(本関数の詳細はg05sj のマニュアルページを参照)

※本サンプルソースコードは .NET環境用の科学技術・統計計算ライブラリである「NAG Library for .NET」の関数を呼び出します。
サンプルのコンパイル及び実行方法


このソースコードをダウンロード
//      g05sj Example Program Text
//      C# version, NAG Copyright 2008
using System;
using NagLibrary;
using System.IO;
namespace NagDotNetExamples
{
  public class G05SJE
  {
    static void Main(String[] args)
    {
      StartExample();
    }
    public static void StartExample()
    {
      try
      {
        const int mseed=1;
        const int n=5;
        double a,   b; int genid,  i,  subid;
        double[] x = new double[n];
        int[] seed = new int[mseed];
        int ifail;
        Console.WriteLine("g05sj Example Program Results");
        Console.WriteLine("");
        //      Set the distribution parameters A and B
        a = 5.00e0;
        b = 1.00e0;
        //      Initialise the seed to a repeatable sequence
        seed[0] = 1762543;
        //      genid and subid identify the base generator.
        genid = 1;
        subid = 1;
        //      Initialise the generator to a repeatable sequence
        G05.G05State g05State = new G05.G05State(genid, subid, seed, out ifail);
        if (ifail != 0)
        {
          Console.WriteLine("** Generator initialisation failed with ifail = {0,5}", ifail);
          goto L20;
        }
        //      Generate the variates
        G05.g05sj(n, a, b, g05State, x, out ifail);
        if (ifail != 0)
        {
          Console.WriteLine("** g05sj failed with ifail = {0,5}", ifail);
          goto L20;
        }
        //      Display the variates
        for (i = 1 ; i <= n ; i++)
        {
          Console.WriteLine(" {0,8:f4}", x[i - 1]);
        }
        Console.WriteLine("");
        // 
        L20: ;
        // 
      }
      catch (Exception e)
      {
        Console.WriteLine(e.Message);
        Console.WriteLine( "Exception Raised");
      }
    }
  }
}


Results matter. Trust NAG.

Privacy Policy | Trademarks