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.

We have obtained the first estimates of the masses of the components of the Her X-1/HZ Her X-ray binary system taking into account non-LTE effects in the formation of the H_gamma absorption line: mx=1.8Msun and mv=2.5Msun. These mass estimates were made in a Roche model based on the observed radial-velocity curve of the optical star, HZ Her. The masses for the X-ray pulsar and optical star obtained for an LTE model lie are mx=0.85\pm0.15Msun and mv=1.87\pm0.13Msun. These mass estimates for the components of Her X-1/HZ Her derived from the radial-velocity curve should be considered tentative. Further mass estimates from high-precision observations of the orbital variability of the absorption profiles in a non-LTE model for the atmosphere of the optical component should be made.

We interpret the observed radial-velocity curve of the optical star in the low-mass X-ray binary 2S 0921-630 using a Roche model, taking into account the X-ray heating of the optical star and screening of X-rays coming from the relativistic object by the accretion disk. Consequences of possible anisotropy of the X-ray radiation are considered.We obtain relations between the masses of the optical and compact (X-ray) components, mv and mx, for orbital inclinations i=60, 75, 90 degrees. Including X-ray heating enabled us to reduce the compact object's mass by near 0.5-1Msun, compared to the case with no heating. Based on the K0III spectral type of the optical component (with a probable mass of mv=2.9Msun, we concluded that mx=2.45-2.55Msun (for i=75-90 degrees). If the K0III star has lost a substantial part of its mass as a result of mass exchange, as in the V404 Cyg and GRS 1905+105 systems, and its mass is $m_v=0.65-0.75Msun, the compact object's mass is close to the standard mass of a neutron star, mx=1.4Msun (for i=75-90 degrees). Thus, it is probable that the X-ray source in the 2S 0921-630 binary is an accreting neutron star.

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.

We have analyzed the observed radial-velocity curve for the X-ray binary M33 X-7 in a Roche model. We have analyzed the dependence between the component masses and the degree of filling of the optical star's Roche lobe to obtain the ratio of the masses of the optical star and compact object. For the most probable mass of the optical star, $m_v=70M_{\odot}$, the mass of the compact object is $m_x=15.55\pm3.20 M_{\odot}$. It has been shown that black holes with masses of $m_x=15M_{\odot}$ and even higher can form in binaries. We present characteristic evolutionary tracks for binary systems passing through an evolutionary stage with properties similar to M33 X-7 - type objects. According to population-synthesis analyses, such binaries should be present in galaxies with masses of at least $10^{11}M_{\odot}$. The present number of such systems in M33 should be of the order of unity. We have also studied the evolutionary status of the X-ray binary IC 10 X-1 with a Wolf-Rayet component, which may contain a massive black hole. The final stages of the evolution of the M33 X-7 and IC 10 X-1 systems should be accompanied by the radiation of gravitational waves.

Using "Scenario Machine" we have carried out population synthesis of radio pulsar with black hole binaries (BH+Psr) in context of the most wide assumptions about star mass loss during evolution, binary stars mass ratio distribution, kick velocity and envelope mass lost during collapse. Our purpose is to display that under any suppositional parameters of evolution scenario BH+Psr population have to be abundant in Galaxy. It is shown that in the all models including models evolved by Heger et al. (2002), Woosley et al. (2002), Heger et al. (2003) expected number of the black holes paired with radio pulsars is sufficient enough to discover such systems within the next few years.