Skew orthogonal polynomials arise in the calculation of the n-point distrib
ution function for the eigenvalues of ensembles of random matrices with ort
hogonal or symplectic symmetry. In particular, the distribution functions a
re completely determined by a certain sum involving the skew orthogonal pol
ynomials. In the case that the eigenvalue probability density function invo
lves a classical weight function, explicit formulas for the skew orthogonal
polynomials are given in terms of related orthogonal polynomials, and the
structure is used to give a closed-form expression for the sum. This theory
treates all classical cases on an equal footing, giving formulas applicabl
e at once to the Hermite, Laguerre, and Jacobi cases.