“…are just like the definitions of [1]. The mathematical controlling model for two-phase incompressible flow in porous media is given by There hold some basic physical conditions: The collocation methods are widely used for solving practice problems in engineering due to its easiness of implementation and high-order accuracy; see, e.g., Douglas and Dupont in [2], Fernandes and Fairweather in [3], and Lu [4][5][6] for linear elliptic and parabolic problems with constant coefficients, Bialecki and Cai in [7] for orthogonal space collocation schemes for elliptic boundary value problems with variable coefficients in two space variable, and Bialecki and Fernandes in [8,9] for orthogonal spline collocation LM, ADI methods, and ADI Crank-Nicolson methods for linear parabolic problems with variable coefficients, respectively. The mathematical controlling model for two-phase incompressible flow in porous media is a strongly nonlinear coupling system of partial differential equations of two different types.…”