In this paper, an investigation and analysis of a nonlinear fractional-order SIR epidemic model with Crowley–Martin type functional response and Holling type-II treatment rate are established along the memory. The existence and stability of the equilibrium points are investigated. The sufficient conditions for the persistence of the disease are provided. First, a threshold value, [Formula: see text], is obtained which determines the stability of equilibria, then model equilibria are determined and their stability analysis is considered by using fractional Routh-Hurwitz stability criterion and fractional La-Salle invariant principle. The fractional derivative is taken in Caputo sense and the numerical solution of the model is obtained by L1 scheme which involves the memory trace that can capture and integrate all past activity. Meanwhile, by using Lyapunov functional approach, the global dynamics of the endemic equilibrium point is discussed. Further, some numerical simulations are performed to illustrate the effectiveness of the theoretical results obtained. The outcome of the study reveals that the applied L1 scheme is computationally very strong and effective to analyze fractional-order differential equations arising in disease dynamics. The results show that order of the fractional derivative has a significant effect on the dynamic process. Also, from the results, it is obvious that the memory effect is zero for [Formula: see text]. When the fractional-order [Formula: see text] is decreased from [Formula: see text] the memory trace nonlinearly increases from [Formula: see text], and its dynamics strongly depends on time. The memory effect points out the difference between the derivatives of the fractional-order and integer order.
Global dynamics of a novel deterministic and stochastic SIR epidemic model with vertical transmission and media coverage