This paper shows a higher order discretization scheme for the Bismut-Elworthy-Li formula, the differentiation of diffusion semigroups. A weak approximation type algorithm with Malliavin weights is constructed through the integration by parts on Wiener space and is efficiently implemented by a...