We report density-functional studies of several small molecules (H 2 , N 2 , CO, H 2 O, and CH 4 ) within the Krieger-Li-Iafrate (KLI) approximation to the exact Kohn-Sham local exchange potential, using a threedimensional real-space finite-difference pseudopotential method. It is found that exchange-only KLI leads to markedly improved eigenvalue spectra compared to those obtained within the standard local-density approximation (LDA), the generalized gradient approximation (GGA), and the Hartree-Fock (HF) method. For structural properties, exchange-only KLI results are close to the corresponding HF values. We find that the addition of LDA or GGA correlation energy functionals to the KLI exact exchange energy functional does not lead to systematic improvements.