n=10; mu=2; mu_plus=3; p=0.2;

  for(i=1; i<=n; i++)
    {
    x1 = gasdev(&idum);
    x2 = gasdev(&idum);
    z  = sqrt(1-p*p)*x1 + p*x2 + mu;
    z_p= x2 + mu_plus;

    if (z > 2.17 || z<-2.17)
     {
        Z[i]=z;
        Z_plus[count]=z_p;
     }
  printf("%d %lf %lfn", i, Z[i], Z_plus[i]);
}

while(count < n)
    {
      x1 = gasdev(&idum);
      x2 = gasdev(&idum);
      z  = sqrt(1-p*p)*x1 + p*x2 + mu;
      z_p= x2 + mu_plus;

   if (z > 2.17 || z<-2.17)
   {
       count++;
       Z[count]=z;
        Z_plus[count]=z_p;
    }
    printf("%d %lfn",count, Z[count]);

    if (count >n) break;
  }