Other Routines

`MakeCircleCoord`	Compute coordinates of a circle placed at a given latitude and longitude.
`MakeEllipseCoord`	Compute coordinates of an ellipse placed at a given latitude and longitude.
`RandomN`	Return a pseudo uniform random deviate between 0 and 1 using the algorithm of Park and Miller with a Marsaglia shift sequence.
`RandomGaussian`	Return a pseudo Gaussian deviate of zero mean and unit variance.
`PreGLQ`	Calculate the weights and nodes used in integrating a function by Gauss-Legendre quadrature.
`EigValVecSym`	Compute the eigenvalues and eigenvectors of a real symmetric matrix.
`EigValVecSymTri`	Compute the eigenvalues and eigenvectors of a real symmetric tridiagonal matrix.
`EigValSym`	Compute the eigenvalues of a real symmetric matrix.
`Wigner3j`	Compute the Wigner-3j symbols for all allowable values of J.
`DHaj`	Compute the latitudinal weights used in the Driscoll and Healy (1994) spherical harmonic transform.