12 float* Gauss(
const float width) {
16 for (
int irand = 0; irand < 100; irand++){
21 frand = frand/RAND_MAX;
22 x1 = 2.0*(frand) - 1.0;
25 frand = frand/RAND_MAX;
29 w = sqrt( (-2.0 * log( w ) ) / w );
30 y[0] = width*x1 * w + 1.0;
31 y[1] = width*x2 * w + 1.0;
38 float* EqualDistribution(
const float width){
42 while ( y[0] < (1.-width) || y[0] > (1.+width) )
44 y[0] = 0.5+rand()/(RAND_MAX*1.0);