d13r2.cpp

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

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

void main()
{
    //program d13r2
    //driver for routine mdian1
    int i,j,npts = 50;
	double xmed;
    double data[51];
    long idum = -5;
    for (i = 1; i<=npts; i++)
	{
        data[i] = gasdev(idum);
    }
    mdian1(data, npts, xmed);
	cout.setf(ios::fixed|ios::right);
	cout.precision(6);
    cout<<endl;
    cout<<"Gaussian distrib., zero mean, unit variance"<<endl;
    cout<<"Median of data set is ";
	cout.width(10);
	cout<<xmed<<endl;
    cout<<"Sorted data"<<endl;
    for (i = 1; i<=10; i++)
	{
        j = (i - 1) * 5;
		cout.width(12);
        cout<<data[j + 1];
		cout.width(12);
        cout<<data[j + 2];
		cout.width(12);
        cout<<data[j + 3];
		cout.width(12);
        cout<<data[j + 4];
		cout.width(12);
        cout<<data[j + 5];
		cout<<endl;
    }
}