// Author:  Douglas Wilhelm Harder
// Copyright (c) 2012 by Douglas Wilhelm Harder.  All rights reserved.

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include "Box_muller.h"

#define N 10000

int main() {
	double s = 0.0;
	double ss = 0.0;

	srand48( time( NULL ) );

	int i;

	for ( i = 0; i < N; ++i ) {
		double x = randn_bm();

		s += x;
		ss += x*x;
	}

	double mean = s/N;

	printf( "Mean:      %f\n", mean );
	printf( "Std dev:   %f\n", sqrt( ss/N - mean*mean ) );

	return 0;
}