We propose a fractional order model for human papillomavirus (HPV) dynamics, including the effects of vaccination and public health education on developing cervical cancer. First, we discuss the general structure of Caputo fractional derivatives and integrals. Next, we define the fractional HPV model using Caputo derivatives. The model equilibrium quantities, with their stability, are discussed based on the magnitude of the reproduction number. We compute and simulate numerical solutions of the presented fractional model using the Adams-Bashforth-Moulton scheme. Meanwhile, real data sourced from reports from the World Health Organization is used to establish the parameters and compute the basic reproduction number. We present figures of state variables for different fractional orders and the classical integer order. The impacts of vaccination and public health education are discussed through numerical simulations. From the results, we observe that an increase in both vaccination rates and public health education increases the quality of life, and thus, reduces disease burden and suffering in communities. The results also confirm that modeling HPV transmission dynamics using fractional derivatives includes history effects in the model, making the model further insightful and appropriate for studying HPV dynamics.
Keywords: Adams-Bashforth-Moulton method; Caputo fractional derivative; Routh-Hurwitz criterion; cervical cancer; human papillomavirus; mean value theorem.