The Lagrange mesh method is a very powerful procedure to compute eigenvalue
s and eigenfunctions of nonrelativistic Hamiltonians. The trial eigenstates
are developed in a basis of well-chosen functions and the computation of H
amiltonian matrix elements requires only the evaluation of the potential at
grid points. It is shown that this method can be used to solve semirelativ
istic two-body eigenvalue equations. As in the nonrelativistic case, it is
very accurate, fast, and very simple to implement.