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