  The method for analytic evaluation of four-particle integrals, proposed by
Fromm and Hill, is generalized to include complex exponential parameters. An
original procedure of numerical branch tracking for multiple valued functions
is developed. It allows high precision variational solution of the Coulomb
four-body problem in a basis of exponential-trigonometric functions of
interparticle separations. Numerical results demonstrate high efficiency and
versatility of the new method.
