This paper proposes a general approximation method for the solutions to second order parabolic partial differential equations (PDEs) by an extension of Leandre's approach and the Bismut identity in Malliavin calculus. We show two types of its applications, new approximations of derivatives...