In this article, for the linearized difference scheme that approximates the Dirichlet problem for the homogeneous multidimensional quasi-linear parabolic equation with unbounded nonlinearity, two-sided point-wise estimates of the solution are established which are fully consistent with the same estimates for the differential problem. It is interesting to note that the proved two-sided estimates do not depend on diffusion coefficient. The direct application of such estimates is the proof of the convergence of the considered difference scheme in the grid norm L 2 . An example of the calculation by the Crank–Nicolson difference scheme is given, showing that the violation of the consistency conditions of differential and difference estimates leads to non-monotonic numerical solutions.