The Watson integral for the d-dimensional hypercubic lattice

and the associated logarithmic integral

are investigated. In particular, a new method is developed which enables one to
calculate the numerical values of {Ld : d = 1,2,...} and
{Wd : d = 3,4,...} with extremely high precision. The asymptotic behaviour of
Ld and Wd as d→∞ is also determined. Finally, some
generalizations of the results are briefly discussed.