#include <iostream>
#include <cassert>
#include <ctime>

#include "Box_muller.h"
#include "Central_limit_theorem.h"
#include "Marsaglia_polar.h"
#include "Ratio.h"

void build_maze( int, int );

int main() {
	srand48( time(0) );
	std::cout.precision( 16 );

	std::cout << "Box-Muller method" << std::endl;

	for ( int i = 0; i < 10; ++i ) {
		std::cout << Normal_distributions::Box_muller::randn() << std::endl;
	}

	std::cout << std::endl << "Central Limit Theorem" << std::endl;

	for ( int i = 0; i < 10; ++i ) {
		std::cout << Normal_distributions::Central_limit_theorem::randn() << std::endl;
	}

	std::cout << std::endl << "Marsaglia polar method" << std::endl;

	for ( int i = 0; i < 10; ++i ) {
		std::cout << Normal_distributions::Marsaglia_polar::randn() << std::endl;
	}

	std::cout << std::endl << "Ratio method" << std::endl;

	for ( int i = 0; i < 10; ++i ) {
		std::cout << Normal_distributions::Ratio::randn() << std::endl;
	}

	return 0;
}