double *helmholtz_exact ( double a, int m, int n, double alpha, double beta, 
  double gamma, int nxy, double x[], double y[] );
void jyndd ( int n, double x, double &bjn, double &djn );
double *jn_zeros ( int n, int nt );