The main aim of this contribution is to construct a numerical scheme for solving stochastic time-dependent partial differential equations (PDEs). This has the advantage problems with positive solutions. provides conditions obtaining solutions, which existing Euler–Maruyama method cannot do. In addition, it more accurate than non-standard finite difference (NSFD) method. Theoretically, suggested...