Abstract
We show the convergence of the Milstein scheme in two-dimensional stochastic differential equations using a Matlab implementation for a particular class of invertible stochastic differential equations. In the implementation, we get an order one approximation under a nondegeneracy assumption.