We have developed a universal approach to compute accurately the brightness of eclipsing binary systems during the transit of a planet in front of the stellar disk. This approach is uniform for all values of the system parameters and applicable to most limb-darkening laws used in astrophysics. In the cases of linear and quadratic limb-darkening laws we obtained analytical expressions for the light curve and its derivatives in terms of elementary functions, elliptic integrals and piecewise-defined function of one variable. In the cases of logarithmic and square root laws of limb darkening the flux and its derivatives were expressed in terms of integrals which can be efficiently computed using Gaussian quadrature formula, taking into account singularities of the integrand.

High-accuracy multicolor light curves of the binary system HD 189733, which contains an exoplanet, are analyzed. We have determined the radii of the star and the planet in the binary system, as well as the orbital inclination. The limb-darkening coefficients of the stellar disk were obtained in 10 filters in the wavelength interval 5500-10500 AA. The uncertainties of the fitted parameters were estimated using the differential-correction and confidence-area methods. The wavelength dependence of the limb-darkening coefficients is compared to the corresponding theoretical function for a model thin stellar atmosphere. We confirm the wavelength dependence of the exoplanet's radius found by Pont et al. The exoplanet radius increases with decreasing wavelength, which seems to argue for the presence of an atmosphere around the planet.