double derf(double);
double derfc(double);
void dgami(double,double,double*,double*,long*);
void dgamik(double,double,double,long);