Building on an insight due to Avramidi, we provide a system of transport equations for determining key fundamental bi-tensors, including derivatives of the world-function, σ (x, x ), the square root of the Van Vleck determinant, ∆ 1/2 (x, x ), and the tail-term, V (x, x ), appearing in the Hadamard form of the Green function. These bi-tensors are central to a broad range of problems from radiation reaction to quantum field theory in curved spacetime and quantum gravity. Their transport equations may be used either in a semi-recursive approach to determining their covariant Taylor series expansions, or as the basis of numerical calculations. To illustrate the power of the semi-recursive approach, we present an implementation in Mathematica which computes very high order covariant series expansions of these objects. Using this code, a moderate laptop can, for example, calculate the coincidence limit [a7(x, x)] and V (x, x ) to order (σ a ) 20 in a matter of minutes. Results may be output in either a compact notation or in xTensor form. In a second application of the approach, we present a scheme for numerically integrating the transport equations as a system of coupled ordinary differential equations. As an example application of the scheme, we integrate along null geodesics to solve for V (x, x ) in Nariai and Schwarzschild spacetimes. * adrian.ottewill@ucd.ie † barry.wardell@aei.mpg.de 1 The Hadamard and DeWitt coefficients also appear in the literature under several other guises. They may be called DeWitt, Gilkey, Minakshisundaram, Schwinger or Seeley coefficients, or any combination thereof (yielding acronyms such as DWSC, DWSG and HDMS).In the coincidence limit, it has been proposed that they be called Hadamard-Minakshisundaram-DeWitt (HaMiDeW) [8] coefficients. For the remainder of this paper, we will refer to them as either DeWitt (for the coefficients a k A B ) or Hadamard (for the coefficients Vr A B ) coefficients.arXiv:0906.0005v3 [gr-qc]