/*
* compute the median of an array of signed integers, using the idea
* in <<Numerical Recipes>>.
*/
static long
median(a0, eps_ptr, x, xlim, gnuf)
double a0; /* initial guess for the median */
float *eps_ptr; /* spacing near the median */
long *x, *xlim; /* the data */
unsigned int gnuf; /* good enough estimate */
{