<abstract><p>The present work implements the continuous Galerkin-Petrov method (cGP(2)-method) to compute an approximate solution of the model for HIV infection of $ \text{CD4}^{+} $ T-cells. We discuss and analyse the influence of different clinical parameters on the model. The work also depicts graphically that how the level of $ \text{CD4}^{+} $ T-cells varies with respect to the emerging parameters in the model. Simultaneously, the model is solved using the fourth-order Runge Kutta (RK4) method. Finally, the validity and reliability of the proposed scheme are verified by comparing the numerical and graphical results with those obtained through the RK4 method. A numerical comparison between the results of the cGP (2) method and the RK4 method reveals that the proposed technique is a promising tool for the approximate solution of non-linear systems of differential equations. The present study highlights the accuracy and efficiency of the proposed schemes as in comparison to the other traditional schemes, for example, the Laplace adomian decomposition method (LADM), variational iteration method (VIM), homotopy analysis method (HAM), homotopy perturbation method (HAPM), etc. In this study, two different versions of the HIV model are considered. In the first one, the supply of new $ \text{CD4}^{+} $ T-cells from the thymus is constant, while in the second, we consider the production of these cells as a monotonically decreasing function of viral load. The experiments show that the lateral model provides more reasonable predictions than the former model.</p></abstract>