/*powell - spectrum minimization routine - logic closely based on Press et al.
 * but totally recoded to be redistributable */
void Powell(float p[], 
  float *xi, 
  long int *n, 
  long int np, 
  double ftol, 
  long int *iter, 
  float *fret);