In this paper, we concern ourselves with the study of a class of stationary states for reaction-diffusion systems
with densities having disjoint supports. Major contribution of this work is computing the numerical solution
of problem as the rate of interaction between two different species tend to infinity. The main difficulty is the
nonlinearity nature of problem. To do so, an efficient iterative method is proposed by hybrid of the radial
basis function (RBF) collocation and finite difference (FD) methods to approximate the solution. Numerical
results with good accuracies are achieved where the shape parameter is carefully selected. Finally, some
numerical examples are given to illustrate the good performance of the method.