Linear, quadratic and square root limb-darkening coefficients have been computed by means of least squares fits to model atmospheres by Kurucz (1991) with the main purpose of facilitating the analysis of light curves of eclipsing binaries. Results are compared with previously published values using the same of different models and various fitting procedures.