A Structure Preserving Numerical Method for Solution of Stochastic Epidemic Model of Smoking Dynamics