Crank-Nicholson difference scheme for the system of nonlinear parabolic equations observing epidemic models with general nonlinear incidence rate

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. © 2021 the Author(s), licensee AIMS Press.