digamma(x_):= psi(x);