This paper deals with the construction of nonstandard finite difference method (NFSD) for nonlinear initial value
problems modeled by a system of nonlinear ordinary differential equations. The proposed scheme preserves the
positivity property as well as the requirement of conservation law and boundedness. In order to illustrate the
accuracy of the new scheme, the numerical results compared with the standard ones.