void dcdchi(double,double,double*,double*,long*);
void dgami(double,double,double*,double*,long*);