Abstract
This study is an attempt to explain a reliable numerical analysis of a stochastic HIV/AIDS model in a two-sex population considering counselling and antiretroviral therapy (ART). The authors are comparing the solutions of the stochastic and deterministic HIV/AIDS epidemic model. Here, an endeavour has been made to explain the stochastic HIV/AIDS epidemic model is comparatively more pragmatic in contrast with the deterministic HIV/AIDS epidemic model. The effect of threshold number H(*) holds on the stochastic HIV/AIDS epidemic model. If H(*) < 1 then condition helps us to control disease in a two-sex human population while H(*) > 1 explains the persistence of disease in the two-sex human population. Lamentably, numerical methods such as Euler-Maruyama, stochastic Euler, and stochastic Runge-Kutta do not work for large time step sizes. The recommended structure preserving framework of the stochastic non-standard finite difference (SNSFD) scheme conserve all vital characteristics such as positivity, boundedness, and dynamical consistency defined by Mickens. The effectiveness of counselling and ART may control HIV/AIDS in a two-sex population.