  In this work we consider a mathematical model of the optical Costas loop. The
pull-in range of the model is estimated by analytical and numerical methods.
Difficulties of numerical analysis, related to the existence of so-called
hidden oscillations in the phase space, are discussed.
