void jyndd ( int n, double x, double &bjn, double &djn );
double *jn_zeros ( int n, int nt );