Abstract
In this work, a mathematical model consisting of a compartmentalized coupled nonlinear system of fractional order differential equations describing the transmission dynamics of COVID-19 is studied. The fractional derivative is taken in the Atangana-Baleanu-Caputo sense. The basic dynamic properties of the fractional model such as invariant region, existence of equilibrium points as well as basic reproduction number are briefly discussed. Qualitative results on the existence and uniqueness of solutions via a fixed point argument as well as stability of the model solutions in the sense of Ulam-Hyers are furnished. Furthermore, the model is fitted to the COVID-19 data circulated by Nigeria Centre for Disease Control and the two-step Adams-Bashforth method incorporating the noninteger order parameter is used to obtain an iterative scheme from which numerical results for the model can be generated. Numerical simulations for the proposed model using Adams-Bashforth iterative scheme are presented to describe the behaviors at distinct values of the fractional index parameter for of each of the system state variables. It was shown numerically that the value of fractional index parameter has a significant effect on the transmission behavior of the disease however, the infected population (the exposed, the asymptomatic infectious, the symptomatic infectious) shrinks with time when the basic reproduction number is less than one irrespective of the value of fractional index parameter.