Abstract
The current effort is devoted to investigating and exploring the stochastic nonlinear mathematical pandemic model to describe the dynamics of the novel coronavirus. The model adopts the form of a nonlinear stochastic susceptible-infected-treated-recovered system, and we investigate the stochastic reproduction dynamics, both analytically and numerically. We applied different standard and nonstandard computational numerical methods for the solution of the stochastic system. The design of a nonstandard computation method for the stochastic system is innovative. Unfortunately, standard computation numerical methods are time-dependent and violate the structure properties of models, such as positivity, boundedness, and dynamical consistency of the stochastic system. To that end, convergence analysis of nonstandard computational methods and simulation with a comparison of standard computational methods are presented.