  The use of Bateman method for solving the two-variable version of the
Lippmann-Schwinger equation without recourse to partial-wave decomposition is
investigated. Bateman method is based on a special kind of interpolation of the
momentum representation of the potential on a multi-variate grid. A suitable
scheme for the generation of a multi-variate Cartesian grid is described. The
method is tested in the nucleon-nucleon scattering employing a model
two-nucleon potential. Our results show that the Bateman method is capable of
producing quite accurate solutions with relatively small number of grid points.
