trigamma(x_):= psi(1,x);