The Lagrange-mesh numerical method has the simplicity of a mesh calculation
and the accuracy of a variational calculation. A three-dimensional Lagrang
e-mesh method based on zeros of Laguerre polynomials is applied tome study
of the ground states of the helium atom, the hydrogen and positronium negat
ive ions and the hydrogen molecular ion. The calculations lead to big spars
e matrices. Highly accurate energies are obtained, which essentially agree
with recently published results or improve them. The accuracy is estimated
by varying the mesh size and two nonlinear scale parameters. Different mean
values of observables are obtained very simply with the Gauss quadrature a
pproximation associated with the mesh.