After breast cancer, cervical cancer is the second most frequent cancer in women globally. The Human Papillomavirus is the most common cause of cervical cancer. In this paper, we used a nonlinear ordinary differential equation system to build a mathematical model of cervical cancer with six compartments (the number of susceptible women, vaccinations of susceptible women, the infected women with HPV, the number of infected with cervical cancer, treatment individual, and recovery class). The model is examined using the existence of bounded and positive solutions, numerical analysis, sensitivity analysis, and stability analysis of disease-free and endemic equilibrium points as a function of R0 values. The numerical simulations of the system are carried out using the ODE45 subroutine of MATLAB and the results are revealed using graphs and biologically interpreted. Using numerical simulation, applying vaccination and increasing treatment for everyone can help to reduce and control the spread of cervical cancer.