# 	psi(x) is digamma(x) function, second polygamma(2,x) function;
#	psi(n,x) is the nth derivative of psi(x), generalized polygamma function, balanced polygamma function;
#	psi(n,x)=d^n/dx^n psi(x);
#psi(-1,x_):=loggamma(x);