In this work, we study second order Crank-Nicholson difference scheme (DS) for the approximate solution of problem (1). The existence and uniqueness of the theorem on a bounded solution of Crank-Nicholson DS uniformly with respect to time step τ is proved. In practice, theoretical results are presented on four systems of nonlinear parabolic equations to explain how it works on one and multidimensional problems. Numerical results are provided.
Keywords: Bounded solution; epidemic models; realization in computer; system of nonlinear parabolic equations.