MODELLING OF THE ENTERPRISE FUNCTIONING STABILITY USING THE AUTOMATIC CONTROL THEORY APPARATUS

Modern Ukrainian society seeks to continually improve the level and conditions of life, which can only be achieved through economic growth of the country. In turn, economic growth can only be achieved through sustainable and stable development of enterprises – manufacturers of products that meet the growing needs of society. On the other hand, in any economic system, long-term economic growth and development cannot be monotonous MODELLING OF THE ENTERPRISE FUNCTIONING STABILITY USING THE AUTOMATIC CONTROL THEORY APPARATUS


Introduction
Modern Ukrainian society seeks to continually improve the level and conditions of life, which can only be achieved through economic growth of the country.In turn, economic growth can only be achieved through sustainable and stable development of enterprises -manufacturers of products that meet the growing needs of society.On the other hand, in any economic system, long-term economic growth and development cannot be monotonous
Crisis phenomena, instability and uncertainty in the economic environment are aggravated by unthoughtful longterm reforms, fierce competition and outdated technologies so characteristic for Ukraine.As a result, a large number of Ukrainian enterprises are characterized by non-fulfillment or failures in the implementation of plans, losses, loss of efficiency, i. e. instability.This leads to insolvency and bankruptcy of such enterprises, which further exacerbates the problems of the Ukrainian economy.Consequently, in today's economic conditions, the problem of the enterprise functioning stability and approaches to stability evaluation becomes particularly relevant.
The overcoming of the above-mentioned problems depends on the timely identification of the enterprise functioning stability threat, as well as the creation of mechanisms for effective management of its sustainable functioning.This becomes possible with the application of economic and mathematical modeling and logistic approach.The result of this approach application should be minimization of costs and increase of the enterprise efficiency.
Taking into account the above-mentioned, the actual task is the development, research and application of a complex of economic-mathematical models for assessment and analysis of the enterprise functioning stability.This complex will allow timely diagnosis of the enterprise instability and making effective management decisions.

Literature review and problem statement
In modern scientific publications, various approaches to modeling the stability of the enterprise functioning are considered.In the works [1,2], the structure of the model basis of the crises early warning and prevention mechanism, aimed at ensuring the stable functioning and development of industrial-economic systems (IES) on different hierarchy levels is proposed.However, the issues related to the construction of the model complex were highlighted fragmentarily.In a number of sources, applied models of wind power plants stability assessment are considered, which allow forming a system of stability diagnostic indicators on the basis of signs complete and incomplete reduction methods [3,4], carrying out classification of IES stability states on the basis of cluster analysis methods [3], DEA methods [5,6], identifying the class of IES stability on the basis of classification trees [3], models of multiple choice [4], discriminant analysis [5].It should be noted that along with the undoubted advantages, the above works do not adequately cover the issues of the comparative analysis of application effectiveness of various simulation methods.In addition, the terms of the model complex use in the IES operation are not fully considered.
The studies related to solving the task for constructing an integral indicator and predicting the level of the IES stability are sufficiently widely presented [7,8].Thus, the possibilities of using taxonomy methods [9,10], Fourier analysis [9], fractal [11], bifurcation analysis [12] are explored.In the paper [13], the methods of crisis dynamics analysis and modeling are complemented by an approach based on the evaluation of the resonance interaction of cycle formation factors.The effectiveness of using the instrumental methods presented above in the comprehensive stability assessment, diagnostics of time series of integrated stability indicators, forecasting of crises is undoubted.However, it would be advisable to consider more fully the issues related to the implementation of a system approach in assessing the IES operation stability.In the works [14,15], a set of models is proposed, which, based on the theories of fuzzy logic, neural networks, allows predicting the enterprise operation stability class.Possibilities for using neural network models to predict stability indicators are also considered in the papers [16,17].In the work [18], a multi-agent system with a batneural network is used for these purposes.To study the influence of factors in changing conditions, the authors of the work [19] suggest the use of a Bayesian network.The above methods have proven effective in investigating mass phenomena in credit scoring systems, bankruptcy monitoring, market value and investment management.However, the application of such an approach in the financial management system of enterprises with a flexible organizational structure is limited to the requirements for information provision.
Under conditions of crisis dynamics, the studies involving the use of a wide range of formalized and non-formalized methods for substantiating the choice of a strategy ensuring sustainable operation and development of enterprises were widely used.In the paper [20], a cognitive model of enterprise stability is proposed, which allows the formation of various scenarios for the development of a production system under various managerial influences and obstacles.In the works [21,22], it is suggested to use simulation modeling based on the concept of system dynamics to justify the choice of a strategy that ensures the stable operation of the IES under conditions of the dominant threats.To substantiate the functional strategies of IES, in the works [23,24], a multiagent approach is used that allows developing flexible models of system behavior taking into account the model agents individual behavior.Considering the undoubted prospects of the authors' research, it should be noted that in the above-mentioned works the issues of defining the parameters at which stable operation of separate subsystems of the IES as well as the IES in a whole are not sufficiently highlighted.
Summarizing the above analysis, it should be noted that, despite the effectiveness of the approaches proposed by the authors, the issues related to the enterprise functioning stability assessment in the dynamics and the effective management of this process remain insufficiently studied.Proposed in scientific publications  models, which are constructed for solving partial problems, solve specific scientific problems and are not aimed at complex solution.In addition, many methods and approaches for determining stability are too complicated, expensive and time-consuming, which is unjustified.
All this necessitates research on the development of economic and mathematical models of stability assessment based on the use of the classical apparatus of the automatic control theory (ACT) and logistic approach in economic systems [25].This approach opens up the broad possibilities for a formal description of such complex systems as industrial enterprises and the application of well-known stability criteria to them.

The aim and objectives of the study
The conducted studies aimed to develop a complex of economic and mathematical models of enterprise operation stability assessment and analysis, using the logistic approach and the apparatus of the automatic control theory.Such models will allow timely diagnosis of the enterprise instability and making effective management decisions.
To achieve the aim, in the process of research the following tasks were solved: -to develop a structural model and models of the enterprise functional subsystems dynamics; -to construct a generalized transfer function of the enterprise and to determine the coefficients of the indicated function; -to form a scenario model for assessing the enterprise operation stability and to build a simulation model for the scenario implementation.

1. Information base used in the study
Parametrization of the enterprise generalized transfer function and the scenario model formation were carried out using the economic indicators of two dairy complexes: State Enterprise "Experimental Farm "Kutuzivka" (Kharkiv region, Ukraine), Institute of Animal Husbandry of UAAS and LLC Agrofirm "Batkivshchyna" (Mykolaiv region, Ukraine).Among the indicators obtained from the financial reports of the complexes in 2016, those were selected that will be used to determine the parameters of the enterprise generalized transfer function coefficients (Table 1).The determination of the demand level and market conditions was carried out on the basis of the dairy products sales prices and milk production volumes in all categories of farms in accordance with the Kharkiv and Mykolaiv regions in 2016, obtained from the State Statistics Service of Ukraine (Table 2) [26].Thus, the information base for constructing the following set of models is the data of agribusiness enterprise financial reporting and the open statistical base of the State Statistics Service of Ukraine.

2. A methodology for constructing a complex of models for estimating and analyzing the enterprise functioning stability
In the conditions of permanent crisis phenomena in the economy, the assessment of the enterprise functioning stability in the dynamics becomes especially important.To solve this problem and to recreate the system functioning algorithm in time, it will be effective to apply simulation modeling.The construction of the simulation model of the enterprise is proposed through a system of differential equations, which is one of the main tools for studying the stability in the classical apparatus of the ACT.This approach becomes possible due to the fact that the structure of the control scheme in systems of different nature (technical, biological, social, economic) is largely similar.Therefore, mathematical models can be used to describe the functioning of economic systems, which describe processes in technical systems.
It was proven that the most effective direction of the organization of enterprises involved in production and sales activities is the formation of a logistic production and sales system (PSS) [27].The construction of such a system was carried out using a logistic approach and represents a solid organizational and economic structure consisting of the enterprise, suppliers of raw materials, components and products, consumers of finished products, and system of transport and warehouse facilities (Fig. 1, a).
The organizational and functional structure of the PSS is a complex cybernetic system with a plurality of material and information flows, managed by managers (block 1) with the help of feedback system.On these feedbacks, operative information arrives.In order to implement the principle of control by deviations, the system is closed with the help of the main feedback.This allows providing the system with the input signal proportional to the output value.Since this signal must act on the system in such a way that the deviation of the output value resulting from the perturbation is reduced, then the main feedback is negative.Along with the main feedback, to improve the accuracy of control in the system, local feedback that covers one or more links is also applied.
The system consists of the complex itself (block 3), warehouses for raw materials, semifinished products, equipment, materials (block 2) and finished products (block 4), external transport (block 5), consumers (block 6), as well as suppliers of raw materials (block 7).With such a structure, the efficiency of the entire PSS is maximized.
System input θ in (Fig. 1, a) is a prognostic production plan that is tracked and implemented by block 1.Output θ out are finished products in the consumer field (block 6).The corresponding indicators of the prognostic production plan go to blocks 2-7, are compared with the possibilities of implementation and are taken into account in the process of PSS functioning.
A somewhat simplified management process of the PSS can be represented as systematic tracking of the system functioning in time by the block 1 when comparing outputs with the actual performance of the prognostic production plan.This comparison occurs for each local subsystem (blocks 2-7), and for the whole system.
In order to construct the simulation model of the enterprise, each block of the organizational-functional structure of the PSS is presented in the form of the corresponding interconnected dynamic links.As a result, we obtain the structural model of the enterprise in terms of the Laplace transformation (Fig. 1, b).
The determination of the types of dynamic links representing the functional blocks of the system on a structural scheme was carried out by compiling the dynamics equations of all the system elements.In drawing up these equations, the essence of economic phenomena occurring in one or another block was studied in detail, the methods of mathematical description of these phenomena were found, and the influence of all operating factors was taken into account.From the equations, it was possible to obtain the transfer and transition functions of the links, to determine transfer coefficients and time constants of the links, and therefore to establish the types of dynamic links (Table 3) [27].Management of the production and sales system (aperiodic link) Warehouse of raw materials, semi-finished products, equipment, materials (real-differentiating link) Production (aperiodic link) External transport (delay link) e -pτ --6 Consumers (aperiodic link)  After collapse of the whole structural model of the enterprise (Fig. 1, b) in accordance with the given rules, collecting like terms and implementation of the substitution, the generalized transfer function of the PSS will be as follows [27]: ( ), ( )..., ( ) ( ), ( )..., ( )

Table 4
The coefficients of the system generalized transfer function Coefficient Formula for the coefficient calculation The structural model of the PSS is converted into one equivalent block (Fig. 2) with a generalized transfer function (1), which is the basis of the simulation model of the enterprise.The input of such a block receives a prognostic production plan, and at the output, we receive finished products in the field of consumption.
For the local criterion of optimality of the PSS functional subsystems in the management process, the minimization of the discrepancy can be accepted: where ∆θ -the value of discrepancy; θ pl out -the prognostic production plan, UAH; θ a out -the actual execution of the plan, UAH.
To determine the parameters of the enterprise structural model coefficients, the theoretical concepts of the modern economic theory and economic systems management organization were used, as well as the indicators presented in the information base of the research.The resulting parameters of the coefficients are shown in Table 5.
To investigate the enterprise functioning stability, the algebraic Hurwitz criterion and the graph-analytical criterion of Mikhailov were used.The choice of these methods is due to the fact that they are developed and successfully used in the ACT, allow investigating the stability in dynamics, and the process of determination of stability is well-algorithmized.
The algebraic Hurwitz criterion allows establishing if the system is stable or not fast, even without solving the characteristic equation.Thus, according to the Hurwitz criterion, the system (enterprise) will be stable, that is, the roots of the characteristic equation obtained from expression (1): a 6 p 6 +a 5 p 5 +a 4 p 4 +a 3 p 3 +a 2 p 2 +a 1 p+a 0 =0, will have negative valid parts if the determinant of Hurwitz and all its diagonal minors (∆ i ) are positive.Thus, the stability conditions have the form: The graph-analytic criterion of Mikhailov provides more information about the system behaviour.In it, the curve that describes the end of the vector obtained from the characteristic equation of the system (3) is used as a determinant, after the substitution p=jω, where 1 = -j : when changing the frequency ω from -∞ to +∞.Since it is proven that the Mikhailov curve is formed as a vector hodograph and symmetric with respect to the abscissa axis, it is sufficient to change ω from 0 to +∞.The Mikhailov's hodograph is constructed on a complex plane in the coordinates P(ω), Q(ω), pre-splitting the expression F( jω) into two parts -valid P(ω) and imaginary Q(ω).In order for a linear automatic control system having a characteristic equation of the n-th order was stable, it is necessary and sufficient to change ω from 0 to ∞ and the complete change of the vector argument F( jω) was equal to , 2 where n -is the degree of the characteristic equation.
In other words, the Mikhailov's curve must be placed in such a way as to cross over sequentially n quadrants of the plane P(ω)-Q(ω).Such a curve (stable system) is called correct.
The formulation of the criterion is retained when placing the roots on the imaginary axis.The presence of roots on the imaginary axis corresponds to the fact that the system is on the boundary of stability and the Mikhailov curve passes through the origin.
For the characteristic equation ( 5) when changing ω from 0 to ∞, the complete change of the vector argument F( jω) will be equal to 3π.That is, for the stable system the Mikhailov curve should be placed so that it would consistently cross six quadrants in the plane P(ω)-Q(ω).
The problem with the coefficients of the characteristic equation, which tend to infinity, was solved by the method of finitization and reduction to the Cauchy task in the construction of the Mikhailov's hodograph.
Finitization is done by replacing the argument ω=t/(1-t), tÎ[0;1] in the characteristic equation (5).The resulting characteristic polynomial has the form: Determination of the angle of the radius vector rotation around the origin of coordinates does not require knowledge of its modulus.Therefore, without losing information about the argument of a complex number, both parts of the equation can be multiplied by the real function (1-t) n .These actions will change the modulus of the complex number, but will not change its argument.Finally, we have the expression: where R * (t) -the finitized radius vector obtained from F( jω).
As a result, for the stability analysis of an equation with constant coefficients, the determination of the radius vector rotation angle R * (t), tÎ[0;1] can be used instead of considering the vector rotation angle F( jω), which modulus tends to infinity.The normalization of the vector was carried out by multiplying the real and imaginary parts of the equation by the fraction (1+t)/R(t).
The curve, which describes the radius vector R * (t), is built on the integrated plane P(t)-Q(t), pre-splitting the expression R * (t) into two parts -valid P(t) and imaginary Q(t).

Results of the enterprise functioning stability studies
The results of calculations of the generalized transfer function coefficients for enterprises, performed in the Mathcad package are given in Table 6.
For the enterprises under study, according to Table 6, the coefficients of the characteristic equation (3) were found and the Hurwitz determinant of the sixth order compiled according to this equation, as well as its five diagonal minors were calculated (Fig. 3).
As can be seen, in the first case, the determinant of Hurwitz and all of its five diagonal minors are positive (Fig. 3, a), that is, the conditions (4) are fulfilled.In the second case, the determinant of Hurwitz and three diagonal minors are negative (Fig. 3, b), therefore, the conditions (4) are not fulfilled.This means that, at the time of modeling, the functioning of the enterprise SE EF "Kutuzivka" is stable under this production plan, and the functioning of LLC Agrofirm "Batkivshchyna", on the contrary, is unstable.
However, as noted above, algebraic criteria allows making only qualitative judgments about the nature of processes, i.e. establishing if the process is stable or not.Concerning how fast the process is damped, the algebraic criterion of Hurwitz does not give an answer.It may turn out that the system is stable, but processes in the system are damped extremely slowly, and such a system turns out to be practically unsuitable.
The Mikhailov curve, constructed according to the data of the investigated enterprises for the characteristic equation (5) when changing ω from 0 to ∞ is shown in Fig. 4. In fact, it is the simulation model of the investigated enterprises built on the basis of generalized transfer functions.
As can be seen from Fig. 4, Mikhailov's hodograph cannot be obtained in the form of a spiral, and the angle of rotation of the radius vector can only be guessed.This is due to the fact that the coefficients of the characteristic equation ( 5) tend to infinity and have a large spread of values among themselves.As a result, with a fairly large ω, the Mikhailov curve also tends to infinity, and the graph does not display its initial areas.

Table 6
The values of the enterprise generalized transfer function coefficients In this regard, it becomes unclear how the curve behaves: how many quadrants crosses and in what sequence.Consequently, it is impossible to make correct conclusions about stability in this case, and it is necessary to go to the adjustment of the simulation model.
The finitized Mikhailov hodograph for the investigated systems of the sixth order, constructed according to the data of the enterprises, is depicted in Fig. 5.
The finitized Mikhailov hodograph for the enterprise SE EF "Kutuzivka" is an integral spiral, which successively crosses six quadrants of the plane P(t)-Q(t), that is, the curve is correct (Fig. 5, a).This means that, according to the Mikhailov's criterion, the functioning of the investigated enterprise is stable, and the adjustment of the simulation model can be considered complete.a b Fig. 5. Finitized Mikhailov's hodograph for the characteristic equation ( 7) when changing t from 0 to 1 built for the enterprises: a -SE EF "Kutuzivka"; b -LLC Agrofirm "Batkivshchyna" For comparison, an appearance of a formalized Mikhailov hodograph built for the enterprise LLC Agrofirm "Batkivshchyna" is shown alongside (Fig. 5, b).As we see, the curve does not intersect six quadrants consistently, but makes a "loop", so it is wrong.This means that the operation of the mentioned enterprise was unsustainable.

Discussion of the enterprise functioning stability results study and the determination of the stability margin
The analysis of the functioning stability of two enterprises by the algebraic Hurwitz criterion and the graph-analytic Mikhailov criterion was carried out.The first is actually a mathematical model for assessing the system's functioning stability, and the second one is a simulation model.Since simulation results coincide in both cases, we can assume that the simulation model constructed is adequate.
It should be noted that in addition to the above cases of stable and unstable systems, it may turn out that the system is on the verge of stability.In this case, there are two possible options of placement of the Mikhailov curve.One corresponds to the presence of imaginary roots (Fig. 6, a), and the other -to the presence of a real zero root (Fig. 6, b).A characteristic feature of the curve placement is that with an infinitely small deformation (shown by strokes in Fig. 6), the correct Mikhailov curve can be obtained.So, if the curve F( jω) passes near the origin of coordinates, we can say that the system is near the limit of stability, that is, we can determine the degree of stability.
For an unstable system, it is also possible that the curve F( jω) will pass through zero, but in this case, an infinitely small deformation does not significantly change the placement of the curve (Fig. 6, c).
The determination of the stability domain by the equations of the first approximation does not give full confidence that the enterprise functioning will be stable at all values of the parameters belonging to the selected area.The calculation of the areas of stability, taking into account the exact equations, is either impossible or loses its practical value as a result of bulkiness.Therefore, it is necessary to perform calculations on the approximate equations, and then introduce the correction coefficients (stability margin).
The need to determine the stability margin is due to the following circumstances: -when compiling the original equations, only the basic laws of economic theory were taken into account and the secondary factors were rejected; -the initial equation is linearized; -economic indicators, which express the time constants and the gain factors of the units, are usually determined with the error; -the calculation is conducted for typical conditions and the structural model of the PSS.
In reality, it is necessary to take into account the statistical nature of the external factors change and the distribution of parameters in different types of systems.Therefore, the application of such a quantitative characteristic as stability margin in the calculations is a kind of guarantee of stability in real conditions.Consequently, in case of negative factors influence, the system will remain stable for a certain time or will be on the verge of stability.How long the system will remain stable depends on the magnitude and duration of these factors.
The stability margin can be expressed in different ways, depending on the criterion, accepted as the basis of the calculation.Let us consider the order of the stability margin using the Mikhailov's criterion.If the system location at the stability limits is indicated by the fact that the curve passes through the origin, then the stability margin is given as the radius of the circle r with a center at the coordinates origin.Moreover, this circle should not cross the Mikhailov curve.
After analyzing the finitized Mikhailov hodograph for the investigated system of the sixth order depicted in Fig. 5, a, it can be revealed that the radius of the circle with the centre at the origin is not more than 0.45 (Fig. 7).Otherwise, the circle will intersect the curve.This fact reflects the presence of a stability margin equal to 0.45.In addition, to determine the stability margin, algebraic criteria also can be used.The stability condition for the Hurwitz algebraic criterion is expressed by the inequality, on the right side of which there is zero (a n-1 a n-2 -a n-3 a n )>0.In order to specify a stability margin, it is needed to take a small positive number instead of zero and write it down (a n-1 a n-2 -a n-3 a n )>e.
Let us find out the stability margin for the SE EF "Kutuzivka".To do this, we will find a solution of the inequality (a 4 a 5 -a 3 a 6 )>0 provided that the transition coefficients K 1 K 2 K 3 K 4 are unknown.Then we will get: If we will analyze the coefficients of this inequality, it becomes clear that it is impossible to define separately K 1 , K 2 or even a product K 1 K 2 K 3 K 4 .However, each term of the inequality includes the coefficients K 1 and K 2 , so we will solve the inequality with respect to the product K 1 K 2 considering K 3 and K 4 as constants.After the corresponding transformation, we get: 2.771×10 13 (K 1 K 2 ) 2 +2.765×10 25 K 1 K 2 + +1.1519×10 19 > 0 .
( 9 ) The roots of the left side of the inequality are: (K 1 K 2 ) I =1.837×10 6 ; (K 1 K 2 ) II =-9.98×10 11 .(10) The second root is rejected because the transition coefficient cannot be negative, and the first root represents the limit value of the product of the coefficients K 1 K 2 at which the system will still be stable: K lim =K 1 K 2 =1.837×10 6 .
Thus, the inequality 1.837×10 6 >K 1 K 2 >0 indicates the condition of the system stability relative to the product K 1 K 2 .In other words, the ratio of supply, demand for goods, overhead, and demand intensity should not exceed 1.837×10 6 .
Similarly, it is possible to determine the boundary values of all coefficients under the condition of the constancy of others.In turn, the determination of the coefficients limit values when assessing the stability margin allows determining the system stability degree.
Consequently, the proposed complex of models for assessing and analyzing the enterprise operation stability is universal, and with appropriate adaptation, opens up wide opportunities for its use for the PSS study in the conditions of any economy branches in different countries.In particular, a set of models can be used to make managerial decisions in real time in order to ensure stable operation and increase the efficiency of enterprises.

Conclusions
1.A complex of economic-mathematical models for the enterprise functioning stability estimation and analysis is developed.The mentioned complex includes the structural model of the enterprise constructed using the logistic approach and the ACT apparatus, as well as partial models of the enterprise functional subsystems dynamics.
2. Based on the structural model of the enterprise, due to the transition from the differential equations of dynamics of elements of the system to algebraic ones, as provided by the ACT, a generalized transfer function of the enterprise in a market environment is constructed.The coefficients of the generalized transfer function of the system are determined.
3. Based on the performance indicators of two dairy complexes, the parametrization of the enterprise generalized transfer function was made and scenario models were formed.Using the algebraic Hurwitz criterion and the graph-analytic criterion of Mikhailov, mathematical and simulation models for the study of the enterprise functioning stability dynamics, based on the generalized transfer function of the system, were constructed.By analyzing possible scenarios of operation, these models allow us to investigate the effect of changes in economic indicators on sustainability and to assess the margin of stability.
materials and components (aperiodic link)

Fig. 1 .
Fig. 1.Logistic production and sales system (PSS): a -scheme of organizational and functional structure; b -structural model; θ in (p) -system input, prognostic production plan (PPP); W 1 (p), W 2 (p),…, W 7 (p) -transfer functions of PSS blocks; ---(2 1) (3 1) (6||7 1) ( ), ( ),..., ( ) fc fc fc W p W p W p -transfer functions of feedback circuits; θ out (p) -system output, implementation of the plan .. a p a p a e b p b p b p W 1-7 (p) -the generalized transfer function of the PSS; θ out (p) -the output (finished products in the field of consumption); θ in (p) -the input (prognostic production plan); A, B, a 6 , a 5 ,…, a 0 , b 4 , b 3 , b 2 -the coefficients of the generalized transfer function; p -the complex Laplace operator; e -the basis of the natural logarithm; τ -the time of cargo transportation.The values of the generalized transfer function coefficients, representing the corresponding combinations of time constants T 1 , T 2 ,..., T 7 and gain factors K 1 , K 2 ,..., K 7 for links in the main circuit and feedback circuits

Fig. 2 .
Fig. 2. Structural model of PSS in the form of an equivalent block

Fig. 6 .
Placement of Mikhailov's hodographs for a system that is on the verge of stability in the presence of: а -complex root; b -zero root; c -hodograph of an unstable system

Fig. 7 .
Fig. 7. Determination of the stability margin using the Mikhailov's criterion

Table 1
Indicators of the enterprises in 2016

Table 2
The external environment of enterprises in 2016

Table 3
Transfer functions of the structural model parts taking into account the Laplace transformation, their coefficients and time constants

Table 5
Parameters of the enterprise structural model coefficients for 2016