RESUMEN
Public health education is pivotal in the management and control of infectious and non-infectious diseases. This manuscript presents and analyses a nonlinear fractional model of tungiasis dynamics with the impact of public health education for the first time. The human population is split into five classes depending on their disease status. The infected population is split into two subgroups; infected but unaware and infected but aware. The model focuses on the impacts of public health education, contact and treatment contact on tungiasis transmission dynamics. Notably, public health education is important for containing as well as reducing disease outbreaks in communities. The Caputo fractional derivative is utilised in defining the model governing equations. Model equilibrium points existence and stability are investigated using simple matrix algebra. Model analysis shows that tungiasis is contained when the reproduction number is less than unity. Otherwise, if it is greater than unity, the disease persists and spread in the population. The generalised Adams-Bashforth-Moulton approach is utilised in solving the derived tungiasis model numerically. The impacts of public health education, treatment and contact rate on overall disease dynamics are discussed through numerical simulations. From the simulations, we see that for given fractional order, public health education and treatment increase the quality of life plus reduce equilibrium numbers of tungiasis-infected individuals. We observe that population classes converge quicker to their steady states when α is increased. Thus, we can conclude that the derivative order α captures the role of experience or knowledge that individuals have on the disease's history.