In this paper we propose a method for stability studies of functional differential systems. The idea of our method is to reduce the analysis of an n-dimensional system to one for an (n + m)-dimensional system, where m is a natural number, to obtain stability and then to come back and make conclusions on the stability of the given n-dimensional system. As an example, a model describing testosterone by distributed inputs feedback control is considered. The aim of the regulation is to hold testosterone concentration above an appropriate level. The feedback control with integral term is proposed. We have to increase the testosterone level to the normal one. The control we proposed could destroy the stability of the model. That is why we have to choose the parameters of our distributed control, namely a dosage or intensity of assimilation of a medicine in a human body in such a form that the stability of our system is preserved. Thus the problem of regulation of testosterone level leads us to the stability analysis of the functional differential system describing a connection between the concentrations of hormones (GnRH), (LH), and testosterone (Te). Constructing the system, we discard the connections which seem nonessential. To estimate the effect of these connections is an important problem. We construct the Cauchy matrix of integro-differential system to estimate this influence.