Abstract
We discuss a fractional-order SIRD mathematical model of the COVID-19 disease in the sense of Caputo in this article. We compute the basic reproduction number through the next-generation matrix. We derive the stability results based on the basic reproduction number. We prove the results of the solution existence and uniqueness via fixed point theory. We utilize the fractional Adams-Bashforth method for obtaining the approximate solution of the proposed model. We illustrate the obtained numerical results in plots to show the COVID-19 transmission dynamics. Further, we compare our results with some reported real data against confirmed infected and death cases per day for the initial 67 days in Wuhan city.