Scaling models of random N x N hermitian matrices and passing to the limit N --> infinity leads to integral operators whose Fredholm determinants describe the statistics of the spacing of the eigenvalues of hermitian matrices of large order. For the Gaussian Unitary Ensemble, and for many others as well, the kernel one obtains by scaling in the ''bulk'' of the spectrum is the ''sine kernel'' sin pi(x-y)/pi(x-y). Rescaling the GUE at the ''edge'' of the spectrum leads to the kernel Ai(x)Ai'(y)-Ai'(x)Ai(y)/x-y , where Ai is the Airy function. In previous work we found several analogies between properties of this ''Airy kernel'' and known properties of the sine kernel: a system of partial differential equations associated with the logarithmic differential of the Fredholm determinant when the underlying domain is a union of intervals; a representation of the Fredholm determinant in terms of a Painleve transcendent in the case of a single interval; and, also in this case, asymptotic expansions for these determinants and related quantities, achieved with the help of a differential operator which commutes with the integral operator. In this paper we show that there are completely analogous properties for a class of kernels which arise when one rescales the Laguerre or Jacobi ensembles at the edge of the spectrum, namely J(alpha)(square-root x) square-root yJ(alpha)'(square-root y) - square-root x J(alpha)'(square-root x) J(alpha) (square-root y)/2(x - y) where J(alpha)(z) is the Bessel function of order alpha. In the cases alpha = -/+ 1/2 these become, after a variable change, the kernels which arise when taking scaling limits in the bulk of the spectrum for the Gaussian orthogonal and symplectic ensembles. In particular, an asymptotic expansion we derive will generalize ones found by Dyson for the Fredholm determinants of these kernels.