There are numerous contagious diseases caused by pathogenic microorganisms, including bacteria, viruses, fungi, and parasites, that have the propensity to culminate in fatal consequences. A communicable disease is an illness caused by a contagion agent or its toxins and spread directly or indirectly to a susceptible animal or human host by an infected person, animal, vector, or immaterial environment. Human immunodeficiency virus (HIV) infection, hepatitis A, B, and C, and measles are all examples of communicable diseases. Acquired immunodeficiency syndrome (AIDS) is a communicable disease caused by HIV infection that has become the most severe issue facing humanity. The research work in this paper is to numerically explore a mathematical model and demonstrate the dynamics of HIV/AIDS disease transmission using a continuous Galerkin–Petrov time discretization of a higher-order scheme, specifically the cGP(2)-scheme. Depict a graphical and tabular comparison between the outcomes of the mentioned scheme and those obtained through other classical schemes that exist in the literature. Further, a comparison is performed relative to the well-known fourth-order Ruge–Kutta (RK4) method with different step sizes. By contrast, the suggested approach provided more accurate results with a larger step size than RK4 with a smaller step size. After validation and confirmation of the suggested scheme and code, we implement the method to the extended model by introducing a treatment rate and show the impact of various non-linear source terms for the generation of new cells. We also determined the basic reproduction number and use the Routh-Hurwitz criterion to assess the stability of disease-free and unique endemic equilibrium states of the HIV model.