polygamma(a_,x_):=psi(a,x); polygamma(x_):= psi(x);