A weak Simpson method for a class of stochastic differential equations and numerical stability results