The Chebyshev-collocation method is briefly presented in a manner favorable for solving ODEs in geophysics depending on geocentric distance and singular at the Earth center. The main advantage of this method is that it is not needed to treat the ODEs particularly at the Earth center as is traditionally done when using the Runge-Kutta method. Comparisons of spheroidal free oscillation periods and load Love numbers with earlier works show very close agreement. ¿ 2001 American Geophysical Union