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.