  The intersection numbers of p-spin curves are computed through correlation
functions of Gaussian ensembles of random matrices in an external matrix
source. The p-dependence of intersection numbers is determined as polynomial in
p; the large p behavior is also considered. The analytic continuation of
intersection numbers to negative values of p is discussed in relation to
SL(2,R)/U(1) black hole sigma model.
