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

<abstract><p>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 $ \tau $ 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.</p></abstract>