d13r5.cpp

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

#include "iostream.h"
#include "math.h"
#include "stdlib.h"

void main()
{
     //program d13r5
     //driver for routine avevar
     int npts = 1000;
     double ave,var,shift,eps = 0.1;
     double data[1001];
     //generate gaussian distributed data
     long idum = -5;
     cout<<endl;
     cout<<"     Shift    Average   Variance"<<endl;
	 cout.setf(ios::fixed|ios::right);
	 cout.precision(2);
     for (int i = 1; i<=11; i++)
	 {
         shift = (i - 1) * eps;
         for (int j = 1; j<=npts; j++)
		 {
             data[j] = shift + i * gasdev(idum);
         }
         avevar(data, npts, ave, var);
		 cout.width(10);
         cout<<shift;
		 cout.width(10);
         cout<<ave;
		 cout.width(10);
         cout<<var<<endl;
	 }
}