In the current work, the deterministic hepatitis B virus epidemic (DHBVE) model and the stochastic hepatitis B virus epidemic (SHBVE) model are two nonlinear mathematical models that serve as the framework to illustrate and predict the dynamic virus behavior of hepatitis B. We employ an approximation based on the outcomes of the deterministic model to solve the stochastic model numerically. Euler-Maruyama method is employed to investigate the SHBVE model, whereas an explicit Runge-Kutta method is exploited to calculate the solution to the DHBVE model. Finally, comparisons between the DHBVE and SHBVE models' frameworks are presented.