Biochemical systems have important practical applications, in particular to understanding critical intra-cellular processes. Often biochemical kinetic models represent cellular processes as systems of chemical reactions, traditionally modeled by the deterministic reaction rate equations. In the cellular environment, many biological processes are inherently stochastic. The stochastic fluctuations due to the presence of some low molecular populations may have a great impact on the biochemical system behavior. Then, stochastic models are required for an accurate description of the system dynamics. An important stochastic model of biochemical kinetics is the Chemical Langevin Equation. In this work, we provide a numerical method for approximating the solution of the Chemical Langevin Equation, namely the derivative-free Milstein scheme. The method is compared with the widely used strategy for this class of problems, the Milstein method. As opposed to the Milstein scheme, the proposed strategy has the advantage that it does not require the calculation of exact derivatives, while having the same strong order of accuracy as the Milstein scheme. Therefore it may be used for an automatic simulation of the numerical solution of the Chemical Langevin Equation. The tests on several models of practical interest show that our method performs very well.
 D. T. Gillespie, “A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions,” Journal of Computational Physics, Vol. 22, No. 4, 1976, pp. 403-434. doi:10.1016/0021-9991(76)90041-3
 C. V. Rao and A. P. Arkin, “Stochastic Chemical Kinetics and the Quasi-Steady-State Assumption: Application to the Gillespie Algorithm,” Journal of Chemical Physics, Vol. 118, No. 11, 2003, pp. 4999-5010. doi:10.1063/1.1545446
 A. Samant and D. Vlachos, “Overcoming Stiffness in Stochastic Simulation Stemming from Partial Equilibrium: A Multiscale Monte-Carlo Algorithm,” Journal of Chemical Physics, Vol. 123, No. 14, 2005, pp. 144114-144122. doi:10.1063/1.2046628
 S. Ilie, “Variable Time-Stepping in the Pathwise Numerical Solution of the Chemical Langevin Equation,” Journal of Chemical Physics, Vol. 137, No. 23, 2012, pp. 234110-234119. doi:10.1063/1.4771660
 K. Burrage, P. M. Burrage and T. Tian, “Numerical Methods for Strong Solutions of Stochastic Differential Equations: An Overview,” Proceedings of the Royal Society A, Vol. 460, No. 2041, 2004, pp. 373-402. doi:10.1098/rspa.2003.1247
 MATLAB, “The Language of Technical Computing,” 2009. www.mathworks.com.