We formulate and analyze an unconditionally stable nonstandard finite difference method for a mathematical model of HIV transmission dynamics. The dynamics of this model are studied using the qualitative theory of dynamical ...