We develop in this work a numerical method for stochastic differential equations (SDEs) with weak second order accuracy based on Gaussian mixture. Unlike the conventional higher schemes SDEs It\^o-Taylor expansion and iterated It\^o integrals, proposed scheme approximates probability measure $\mu(X^{n+1}|X^n=x_n)$ by mixture of Gaussians. The solution at next time step $X^{n+1}$ is then drawn f...