The DeWitt expansion of the matrix element M xy = x| exp −[ 1 2 (p − A) 2 + V ]t |y , (p = −i∂) in powers of t can be made in a number of ways. For x = y (the case of interest when doing one-loop calculations) numerous approaches have been employed to determine this expansion to very high order; when x = y (relevant for doing calculations beyond one-loop) there appear to be but two examples of performing the DeWitt expansion. In this paper we compute the off-diagonal elements of the DeWitt expansion coefficients using the Fock-Schwinger gauge. Our technique is based on representing M xy by a quantum mechanical path integral. We also generalize our method to the case of curved space, allowing us to determine the DeWitt expansion ofM xy = x| exp 1t|y by use of normal coordinates. By comparison with results for the DeWitt expansion of this matrix element obtained by the iterative solution of the diffusion equation, the relative merit of different approaches to the representation of M xy as a quantum mechanical path integral can be assessed. Furthermore, the exact dependence ofM xy on some geometric scalars can be determined. In two appendices, we discuss boundary effects in the one-dimensional quantum mechanical path integral, and the curved space generalization of the Fock-Schwinger gauge.