An accurate characteristics-splitting scheme for numerical solution of MHD equations