When approximating multidimensional partial differential equations, the values of grid functions from neighboring layers are taken previous time layer or approximation. As a result, along with approximation discrepancy, an additional discrepancy numerical solution is formed. To reduce this when solving stationary elliptic equation, parabolization carried out and resulting equation solved by met...