Numerical modelling of the perturbation of HIV-1 during combination anti-retroviral therapy

A competitive, chaos-free, implicit, finite-difference method is developed and used for a novel deterministic model for the perturbation of HIV by combination antiretroviral therapy. The compartmental model monitors the interaction between HIV and CD4⁺ T cells, its principal target and site of replication in vivo, in the presence of reverse transcription inhibitors and protease inhibitors. The model exhibits two steady states, an uninfected (trivial) steady state (with no virus present) and an endemically infected state (with virus and infected T cells present). Stability and bifurcation analyses together with numerical simulations of the resulting dynamical system are reported.