Efficient recursive algorithms for functionals based on higher order derivatives of the multivariate Gaussian density