dist_uniformtest.cpp

随机数 基于c++的Monte Carlo源代码

/*

Send a stream number as the first command line argument when running
this code.  Suppose the binary of this program is named "main".  Then
type "./main 2 >> rand_2.txt".  Then run Octave and type "load
rand_2.txt".  Finally, type "hist(M)" to see a relatively flat
histogram of the values produced by this program.  Slightly different
histograms are produced by different streams, but this is to be
expected.

*/

#include <iostream>
#include <cstdlib>

#include "../include/mcsLibre.h"

int main(int argc, char *argv[])
{
  mcsLibre mytest;

  int samples = 100000;
  int i;
  double V;

  cout << "# Created by Octave 2.1.33, Wed Feb 06 19:40:34 2002 PST <[e-mail]>" << endl;
  cout << "# name: M" << endl;
  cout << "# type: matrix" << endl;
  cout << "# rows: " << samples << endl;
  cout << "# columns: 1" << endl;

  for (i=0; i<samples; i++)
    {
      V = mytest.uniform(3, 8, atoi(argv[1]));
      cout << V << endl;
    }

  return 0;
}