By the method of finite differences, MDF, give a good numerical approximation of the solution of the partial differential equation modeling the flow of contaminants from an aquifer. For this, we discretize the domain and using appropriate approaches to the partial derivatives of the equation rewrite the problem as an algebraic system of equations. Examine the effectiveness of two finite difference schemes: an explicit method and the Crank-Nicolson method for the problem to be solved by analyzing their consistency, stability and convergence. The problem is computationally implemented using scientific software