Abstract
The resolution of a system for modeling the competition between opponents in a chemostat when one of these can produce a toxin has been studied. We propose a novel method to overcome the analytical difficulties of standard mathematical methods. The method is based on the variational iteration method and combined with the Gauss-Seidel technique for increasing the convergence rate. Numerical examples are considered to demonstrate the practicality and improve the convergence of the proposed method.