A numerical approach for a model of the precancer lesions caused by the human papillomavirus

We develop two numerical methods to approximate the solutions of a pioneer model of the lesions at the cervical cells caused by the human papillomavirus. Such model is given by a nonlinear advection–diffusion-reaction partial differential equation and the goal of the schemes is to analyze the behaviour of the evolution of infected cells. The developed schemes consist of two explicit non standard finite differences numerical schemes which satisfy positivity conditions. They are based on the subequation method in the context of the non standard scheme methodology. Our approach provides an alternative method to the early diagnosis of the disease and may open up new lines of research.