In this paper, we discuss the patterns of the transmission of how the Human papillomavirus (HPV) spreads in the public. We analyze the behaviors of infectious disease in SIR, SEIR, and Multistrian models. In terms of this disease, HPV has many subtypes, some are high-risk, and some are low-risk. Unfortunately, many people in some areas are not aware of it. That will always cause many delays for the disease and aggravations. The harm caused by HPV is different for different groups of people, people will take various measures to handle the situation according to the severity. We divided them into 3 cases. That is case1-negative for people who may have caught HPV; case 2-negative people who may have caught HPV, but excluding exposures; case 3-negative for people who is sure to catch HPV. After that, by computing the disease-free equilibrium (DFE) to get the value for calculating reproduction number
R
0
to estimate the stability for every case. Then the way of getting the derivatives and partitioned matrix and computing reproduction number with the measure of the Jacobian matrix were used for deciding the response for the epidemic eventually.