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

void main()
{
    //program d13r3
    //driver for routine mdian2
    int i,npts = 50;
	double xmed;
    double data[51];
    long idum = -5;
    for (i = 1; i<=npts; i++)
	{
        data[i] = gasdev(idum);
    }
    mdian2(data, npts, xmed);
	cout.setf(ios::fixed|ios::left);
	cout.precision(6);
    cout<<endl;
    cout<<"Gaussian distrib., zero mean, unit variance"<<endl;
    cout<<"Median according to mdian2 is   ";
	cout.width(12);
    cout<<xmed<<endl;
    mdian1(data, npts, xmed);
    cout<<"Median according to mdian1 is   ";
	cout.width(12);
    cout<<xmed<<endl;
}