void boxrnd(double m,double sig)
{
int r1,r2;

r1=rand();
r2=rand();
x =sig*sqrt(-2*log(r1))*cos(2*3.14159*r2)+m;
y =sig*sqrt(-2*log(r1))*sin(2*3.14159*r2)+m;

}