This paper aims to present a nonstandard finite difference (NFSD) method to solve of mathematical model of Hepatitis B disease. The proposed method has important properties such as positivity and stability. Numerical comparisons confirm that the accuracy of our method is better than other existence standard methods such as the four-order Runge-Kutta (RK4) method and the Euler method.