  A method to compute the bound state eigenvalues and eigenfunctions of a
Schr\"{o}dinger equation or a spinless Salpeter equation with central
interaction is presented. This method is the generalization to the
three-dimensional case of the Fourier grid Hamiltonian method for
one-dimensional Schr\"{o}dinger equation. It requires only the evaluation of
the potential at equally spaced grid points and yields the radial part of the
eigenfunctions at the same grid points. It can be easily extended to the case
of coupled channel equations and to the case of non-local interactions.
