The numerical evaluation of coherent-state path-integral representations for partition functions and other quantities in equilibrium quantum statistical mechanics is discussed. Several coherent-state path-integral schemes are introduced, which differ from each other by the order of approximation and by the operator ordering employed in the high-temperature approximation of the density operator. Simple one-dimensional systems are used to test these schemes. For the harmonic oscillator, finite-dimensional approximations to the coherent-state path integral are calculated analytically and compared to each other and to the real-space path integral. For anharmonic systems, integrations must be approximated by quadrature formulas. This leads to a matrix multiplication scheme which is tested for the double-well potential. The results obtained are accurate from zero temperatures way up into the high-temperature regime where quantum effects become negligible. This is a significant advantage over traditional real-space path integral schemes which break down at low temperatures.