An efficient method for solving difference systems for elliptic differential equations

A method for solving systems of linear algebraic equations arising in connection with the approximation of boundary value problems for elliptic partial differential equations is proposed. This method belongs to the class of conjugate directions method applied to a preliminary transformed system of equations. A model example is used to explain the idea underlying this method and to investigate it. Results of numerical experiments that confirm the method’s efficiency are discussed.