A monotone approximation is proposed for a system without monotone non
linearity. A new concept of ordered pair of supersolution and subsolution is intro
duced, and then the existence of numerical solutions is studied. A new monotone
iteration is provided for solving the resulting problem. An approximation with
high accuracy is investigated. The corresponding iteration possesses geometric
convergence rate. The numerical results support the theoretical analysis.