Abstract
Stiff stochastic differential equations arise in many applications including in the area of biology. In this paper, we present numerical solution of stochastic differential equations representing the Malthus population model and SIS epidemic model, using the improved implicit Milstein method of order one proposed in [6]. The open source programming language SCILAB is used to perform the numerical simulations. Results show that the method is more accurate and stable compared to the implicit Euler method.