In this paper, a fractional model in the Caputo sense is used to characterize the dynamics of HPV with cervical cancer. Generalized mean value theorem has been used to examine whether the infection model has a unique positive solution. The model has two equilibrium points: the disease-free point and the endemic point. The examination of the system’s local and global stability is provided in terms of the basic reproductive number Rp°. The global stability analysis has been carried out using an appropriate Lyapunov function and the LaSalle invariant principle. The results demonstrate that in the infection model, if Rp°<1, then the solution converges to the disease-free equilibrium, which is both locally and globally asymptotically stable. Whilst Rp°>1, the endemic equilibrium is considered to exist. Simulations are implemented via a finite difference method with Grünwald-Letnikov discretization approach for Caputo derivative operator to define how changes in parameters impact the dynamic behavior of the system using Matlab.