extended Raviart–Thomas mixed finite element method. In this method, the numerical
solution of this equation is obtained using triangular meshes. Also, for discretization in time
direction, we use an implicit finite difference scheme. In addition, error estimation and
stability analysis of both methods are shown. Finally, some numerical examples are
considered to confirm the theoretical results. Copyright© 2017 John Wiley & Sons, Ltd.