Generating functions for intersection numbers on moduli spaces of curves

Using the connection between intersection theory on the Deligne-Mumford spaces and the edge scaling of the GUE matrix model (see math.CO/9903176, math.AG/0101147), we express the n-point functions for the intersection numbers as n-dimensional error-function-type integrals and also give a derivation of Witten's KdV equations using the higher Fay identities of Adler, Shiota, and van Moerbeke.