d13r4.cpp

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

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

void main()
{
    //program d13r4
    //driver for routine ttest
    int npts = 1024, nshft = 10;
    int i,j,mpts = 512;
    double eps = 0.02;
    double data1[1025], data2[513];
    //generate gaussian distributed data
    long idum = -5;
    for (i = 1; i<=npts; i++)
	{
        data1[i] = gasdev(idum);
    }
    idum = -11;
    for (i = 1; i<=mpts; i++)
	{
        data2[i] = (nshft / 2.0) * eps + gasdev(idum);
    }
    cout<<endl;
    cout<<"      shift       t     probability";
    cout<<endl;
	double t, prob,shift;
	cout.setf(ios::fixed|ios::right);
	cout.precision(2);
    for (i = 1; i<=nshft + 1; i++)
	{
        ttest(data1, npts, data2, mpts, t, prob);
        shift = (i - 1) * eps;
		cout.width(10);
        cout<<shift;
		cout.width(10);
		cout<<t;
		cout.width(10);
        cout<<prob<<endl;
        for (j = 1; j<=npts; j++)
		{
            data1[j] = data1[j] + eps;
        }
    }
}