In this talk, we show a general formula for an any order weak approximation of stochastic differential equations using multidimensional Malliavin weights. Moreover, we give a version of weak approximation by applying Bismut-Elworthy-Li type formula and a computation method for distributions on Wiener space. Numerical examples for SABR stochastic volatility model confirm the validity of the scheme.