//***********************************
double CPSD060415Dlg::Data(double a,double b,long int *seed)  
{
	using namespace std;
//double a,b;
//long int *seed;
{
	double t;
	*seed=2045*(*seed)+1;
	*seed=*seed-(*seed/1048576)*1048576;
	t=(*seed)/1048576.0;
	t=a+(b-a)*t;
	return(t);  }

}
//**********************************
//**********************************
double CPSD060415Dlg::gauss(double mean,double sigma,long int *s)
{
	using namespace std;
//long int *s;

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