d6r8.cpp

Visual C++ 常用数值算法集 源代码

#include <iostream.h>
#include <math.h>
#include <iomanip.h>
#include <stdlib.h>

void main()
{
    //program d6r8
    //driver for routine poidev
    int idum,n = 20, npts = 10000, iscal = 200, llen = 50;
	int xm,klim,i,j,k,m;
    double dist[21];
	char text[50];
    idum = -13;
    do
	{
        for (j = 1; j<=21; j++)
            dist[j] = 0.0;
        cout<<"Order of Poisson distribution [n=1..20] -1 to end."<<endl;
        xm = 5;
        cout<<xm<<endl;
        if (xm <= 0) exit(1);
	}
    while (xm > 20);
    for (i = 1; i<=npts; i++)
	{
        j = int(poidev(xm, idum)) + 1;
        if ((j >= 1) && (j <= 21)) dist[j] = dist[j] + 1;
    }
    cout<< "Poisson distribution deviate, order "<<xm<<"of "<< npts<< " points"<<endl;
    cout<< "    x     p[x]     graph:"<<endl;
    for (j = 1; j<=20; j++)
	{
        dist[j] = dist[j] / npts;
        for (k = 0; k<=50; k++)
            text[k] =' ';
        klim = int(iscal * dist[j]);
        if (klim > llen) klim = llen;
        for (k = 1; k<=klim; k++)
            text[k] = '*';
        cout<< setw(5)<<float(j);
        cout<< setw(10)<<dist[j]<<"    ";
        
        for (m = 1; m<=50; m++)
            cout<< text[m];
		cout<<endl;
     }
}