#include <stdlib.h>
#include <time.h>
#include <math.h>
#include "randn.h"

double randn(const double Mean, const double SquareMargin)
{
	const double PI = 3.1415926;
	double model = 4294967296.0;
	double multiplicator = 663608941.0;
	double temp_a, temp_b;
	static double seed = (double)time(NULL);
	
	seed *= multiplicator;
	seed -= (int)(seed/model) * model;
	temp_a = seed/model;

	seed *= multiplicator;
	seed -= (int)(seed/model) * model;
	temp_b = seed/model;

	double norm_randn = sqrt(-2*log(temp_a))*cos(2*PI*temp_b);
	double gauss_rand = Mean + norm_randn*SquareMargin;
	return(gauss_rand);
}