//gauss_gen.c



double gauss(double mean,double sigma,long int *s)
{	
	double uniform(double,double,long int *);
	
	int i;
	double x,y;
	for(x=0,i=0;i<12;i++)
		x+=uniform(0.0,1.0,s);
	x=x-6.0;
	y=mean+x*sigma;	
	return(y);
}