  We model X-ray luminosity functions (XLF) of accreting neutron stars and
black holes in $10^{35} \leq L_X \leq 10^{41}$ erg/s range in star-forming
galaxies and galaxies with the initial star formation burst. XLFs are obtained
by combining a fast generation of compact object+normal star population using
the binary population synthesis code BSE and calculation of the subsequent
detailed binary evolution by the MESA code. XLFs in the galaxies of both types
is broadly reproduced using the standard assumptions of the binary star
evolution.
