BEGIN { Machines = 10
	Mean     = 10
	Sd       = Sd ? Sd : 0.25
 	N=10000; while(N--) {
 	    Show = Sep = ""
 	    for(M=1;M<=Machines;M++)  {
		S    = int(0.5+ normal(Mean,Sd))
   	        Show = Show Sep S
		Sep = ","
            }
	    print Show   
	}
}
function any(n) {return int(0.5 + n*rand())  }

function gaussianPdf(mean,standardDev,x,  pi) {
    pi= 1068966896 / 340262731; #: good to 17 decimal places
    e = 2.71828182845904523;
    return 1/(standardDev*sqrt(2*pi)) * \
               e^ (-1*(x-mean)^2/(2*standardDev*standardDev))
}
function normal(m,s) {    return m+box_muller()*s; }
function box_muller(m,s,    n,x1,x2,w) {
    w=1;
    while (w >= 1) {
	x1= 2.0 * rand() - 1;
	x2= 2.0 * rand() - 1;
	w = x1*x1 + x2*x2};
    w = sqrt((-2.0 * log(w))/w);
    return x1 * w;
}