Non-Standard Finite Difference Schemes for Investigating Stability of a Mathematical Model of Virus Therapy for Cancer

In this paper, a special case of finite difference method called non-standard finite difference (NSFD) method was studied to compute the numerical solutions of the nonlinear mathematical model of the interaction between tumor cells and oncolytic viruses. The global stability of the equilibrium points of the discrete model is investigated by using the Lyapunov stability theorem. Some conditions were gained for the local asymptotical stability of the equilibrium points of the system. Finally, numerical simulations are carried out to illustrate the main theoretical results. The discrete system is dynamically consistent with its continuous model, it preserves essential properties, such as positivity, boundedness of the solution, stability properties of the equilibrium points