Two phase flow problem in a porous medium is governed by a system of nonlinear equations. One is an elliptic equation for the pressure and the other is a parabolic equation for the concentration of one of the fluids. Since the pressure appears in the concentration only through its velocity field, we choose a mixed finite element to approximate the pressure equation and use the standard Galerkin...