Lyapunov functions for fractional-order systems in biology: Methods and applications

A Fractional-Order Epidemic Model with Quarantine Class and Nonmonotonic Incidence: Modeling and Simulations

In any outbreak of infectious disease, the timely quarantine of infected individuals along with preventive measures strategy are the crucial methods to control new infections in the population. Therefore, this study aims to provide a novel fractional Caputo derivative-based susceptible-infected-quarantined-recovered-susceptible epidemic mathematical model along with a nonmonotonic incidence rate of infection. A new quarantined individual compartment is incorporated into the susceptible-infected-recovered-susceptible compartmental model by dividing the total population into four subpopulations. The nonmonotonic incidence rate of infection is considered as Monod–Haldane functional type to understand the psychological effects in the population. Qualitative analysis of the study shows that the model solutions are well-posed i.e., they are nonnegative and bounded in an attractive region. It is revealed that the model has two equilibria, namely, disease-free (DFE) and endemic (EE). The stability analysis of equilibria is investigated for local as well as global behaviors. Mathematical analysis of the model reveals that DFE is locally asymptotically stable when the basic reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$({R}_{0})$$\end{document}(R0) is lower than one. The basic reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{0}$$\end{document}R0 is computed using the next-generation matrix method. The existence of EE is shown and it is investigated that EE is locally asymptotically stable when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{0}>1$$\end{document}R0>1 under some appropriate conditions. Moreover, the global stability behaviors of DFE and EE are analyzed under some conditions using \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{0}$$\end{document}R0. Finally, some numerical simulations are performed to interpret the theoretical findings.

Monic Chebyshev pseudospectral differentiation matrices for higher-order IVPs and BVPs: applications to certain types of real-life problems

We introduce new differentiation matrices based on the pseudospectral collocation method. Monic Chebyshev polynomials (MCPs) were used as trial functions in differentiation matrices (D-matrices). Those matrices have been used to approximate the solutions of higher-order ordinary differential equations (H-ODEs). Two techniques will be used in this work. The first technique is a direct approximation of the H-ODE. While the second technique depends on transforming the H-ODE into a system of lower order ODEs. We discuss the error analysis of these D-matrices in-depth. Also, the approximation and truncation error convergence have been presented to improve the error analysis. Some numerical test functions and examples are illustrated to show the constructed D-matrices’ efficiency and accuracy.