ABSTRACT:  Over time, several methods were developed to deal with neutral particle transport problems. The interest in these problems is related to their wide range of applications, from neutron transport and heat transfer in nuclear reactors to radiative transfer in atmospheric clouds. Unlike the discrete ordinates or discrete ordinates–like methods, integral methods do not require discretization of angular variables. Instead, angular variables are completely eliminated by an integration procedure over the solid angle, which allows elimination of the ray effect. That said, this paper presents a new approach to estimate the scalar flux in two dimensional fixed-source neutron transport problems in a heterogeneous medium, considering isotropic scattering and vacuum and reflective boundary conditions. Here, the Nyström method is combined with the singularity-subtraction technique to present an integral formulation for the scalar flux in a mesh grid over all regions of the domain. The iterative method of the Neumann series is used as an alternative to direct methods to solve the resulting system of equations generated from the domain discretization. Numerical results are given to verify the offered method.