We study the numerical solution of the nonrelativistic Schrödinger equation for two-electron atoms in ground and excited S states using pseudospectral (PS) methods of calculation. The calculation achieves convergence rates for the energy, Cauchy error in the wave function, and variance in local energy that are exponentially fast for all practical purposes. The method requires three separate subdomains to handle the wave function's cusplike behavior near the two-particle coalescences. The use of three subdomains is essential to maintaining exponential convergence and is more computationally efficient than a single subdomain. A comparison of several different treatments of the cusps suggests that the simplest prescription is sufficient. We investigate two alternate methods for handling the semi-infinite domain, one which involves a sequence of truncated versions of the domain and the other which employs an algebraic mapping of the semi-infinite domain to a finite one and imposes no explicit cutoffs on the wave function. The latter prescription proves superior. For many purposes it proves unnecessary to handle the three-particle coalescence in a special way. The presence of logarithmic terms in the exact solution is expected to limit the convergence to being nonexponential but the only clear evidence of that is the rate of convergence of derivatives near the three-particle coalescence point. Higher resolution than achieved in this work will ultimately be needed to see its limiting effect on other measures of error. As developed and applied here the PS method has many virtues: no explicit assumptions need be made about the asymptotic behavior of the wave function near cusps or at large distances, the local energy (Hψ/ψ) is exactly equal to the calculated global energy at all collocation points, local errors go down everywhere with increasing resolution, the effective basis using Chebyshev polynomials is complete and simple, and the method is easily extensible to other bound states. As the number of collocation points grows, the method achieves exponential convergence up to the resolution tested. This study serves as a proof-of-principle of the method for more general two-and possibly three-electron applications.