A mathematical model, composed of two non-linear differential equations that describe the population dynamics of CD4 T cells in the human immune system, as well as viral HIV particles, is proposed. The invariance region is determined, classical equilibria stability analysis is performed using the basic reproduction number, and numerical simulations are carried out, in order to illustrate stability results. Later, the model is modified with a delay term, which describes the time that cells require for immunological activation. This generates a two-dimensional integro-differential system, which is transformed into a system with three ordinary differential equations, via auxiliary variable use. For the new model, equilibrium points are determined, their local stability is examined, and results are studied by way of numerical simulation.