srandomdev();                       //always return zero
randinit(random());                 //to change the freefem seed
gslrng r;
gslrngset(r, random());             // to change the r seed of gsl

real ran = gslranflat(r,1,0);  //unifor distributed