Abstract
In this study, an integer-order rabies model is converted into the fractional-order epidemic model. To this end, the Caputo fractional-order derivatives are plugged in place of the classical derivatives. The positivity and boundedness of the fractional-order mathematical model is investigated by applying Laplace transformation and its inversion. To study the qualitative behavior of the non-integer rabies model, two steady states and the basic reproductive number of the underlying model are worked out. The local and global stability is investigated at both the steady states of the fractional-order epidemic model. After analytic treatment, a structure-preserving numerical template is constructed to numerically solve the fractional-order epidemic model. Moreover, the positivity, boundedness and symmetry of the numerical scheme are examined. Lastly, numerical experiment and simulations are accomplished to substantiate the significant traits of the projected numerical design. Consequences of the study are highlighted in the closing section.