Two-phase, incompressible, immiscible flow in porous media is governed by a coupled system of nonlinear partial differential equations. The pressure equation is elliptic, whereas the concentration equation is parabolic, and both are treated by the collocation scheme. Existence and uniqueness of
solutions of the algorithm are proved. A optimal convergence analysis is given for the method.