Computation of the derivative of the Hurwitz zeta-function and the higher Kinkelin constants

I use the numerical values of the generalised Glaisher-Kinkelin-Bendersky (GKB) constants to give numerical values for the derivatives of the Hurwitz zeta function at negative integers, rather than the other way round. I point out that both Glaisher's numerical approach and Bendersky's recursion for the generalised gamma function were anticipated by Jeffery in 1862 who gave the value of the second constant as an example. I therefore propose that GKB become GKBJ.