polygamma(x, delta^2/2+lambda^2/2+1) for x=0