int main() { double d = psigamma(2.0, 4); return d != 0 ? 0 : 1; }
/* n'th order derivative of log gamma function */ double D_lgamma(double x, double n){ if(n<.5)return lgammafn(x); else return psigamma(x,n-1.0); }