double dgam1(double);
double dgamma(double);