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

void main()
{
    //program d13r6a
    //driver for routine tutest
    int i,j,npts = 5000, nshft = 10;
    int mpts = 1000;
    double t,prob,shift,eps = 0.02;
    double var1 = 1.0, var2 = 4.0;
    double data1[5001], data2[1001];
    //generate two gaussian distributions of different variance
    long idum = -51773;
    double fctr1 = sqrt(var1);
    for (i = 1; i<=npts; i++)
	{
        data1[i] = fctr1 * gasdev(idum);
    }
    double fctr2 = sqrt(var2);
    for (i = 1; i<=mpts; i++)
	{
        data2[i] = (nshft / 2.0) * eps + fctr2 * gasdev(idum);
    }
    cout<<endl;
    cout<<"Distribution #1: variance = ";
	cout<<var1<<endl;
    cout<<"Distribution #2: variance = ";
	cout<<var2;
    cout<<endl;
    cout<<"      Shift      T     Probability"<<endl;
	cout.setf(ios::fixed|ios::right);
	cout.precision(2);
    for (i = 1; i<=nshft + 1; i++)
	{
        tutest(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;
		}
    }
}