double dcsevl(double,double[],long);
double dgam1(double);
double dgamma(double);
void dinits(double[],long,double,long*);
double drcomp(double,double);
double drlog(double);
double dxparg(long);