A study is made of the numerical evaluation to high precision

of the three branches of the Lambert W-function and of related

definite integrals. The latter are evaluated by nonstandard

Gaussian quadrature based on known integral

representations or by explicit formulae involving the

trigamma function.