An algorithm for the numerical solution of a system of two parabolic equations of a special form has been presented, and its stability has been proved. Estimates of the numerical solutions obtained have been made. Conditions guaranteeing the applicability of the sweep method are determined, and the time step at which the iterative method has the inner convergence has been estimated. An example of the application of the algorithm in soil science has been presented.