Abstract: | A chaos-free numerical method will be developed for the solution of a system of non-linear initial-value problems (IVP's) associated with the transmission dynamics of two HIV subtypes. It will be shown that integration of this four-dimensional system with some standard numerical integrators like the fourth-order Runge-Kutta algorithm leads to scheme-dependent numerical instabilities. |