void dbesjn(double,double,long,double[]);
void dbespq(double,double,double*,double*);
void derv1(STRING,double,byte);
double dgamma(double);
void ermsg(STRING,long,long,STRING,byte);
void ierv1(STRING,long,byte);