Modeling and Mathematical Analysis of Labor Force Evolution

In this paper, we propose a delayed differential system to model labor force (occupied labor force and unemployed) evolution.,e mathematical analysis of our model focuses on the local behavior of the labor force around a positive equilibrium position; the existence of a branch of periodic solutions bifurcated from the positive equilibrium is then analyzed according to the Hopf bifurcation theorem. Finally, we performed numerical simulations to illustrate our theoretical results.


Introduction
Mathematical modeling has gained considerable importance in recent decades in all fields, particularly in the macroeconomic field.is branch of applied mathematics is interested in the representation, with enough loyalty, of the real phenomena by abstract mathematical objects (variables, parameters, and equations) accessible to analysis and calculation.Its main objective is the exploitation of the conditions of yesterday and/or today to predict those of tomorrow (for example, we can use ordinary differential equations or delayed differential equations).
is work is an opportunity to take stock of the evolution of the labor force over time to explain the reasons for the fluctuations.To do this, we propose delayed differential equations to model labor force evolution (knowing the past and the present, foresee the future).
e study of the labor force evolution (occupied labor force and unemployed) finds all its relevance in the economic and mathematical researches for two reasons, the first relating to the role of labor force in the national economy and the second linking up to the fluctuations in labor force, that occur at different times, with different durations and amplitudes [1][2][3][4][5].ese fluctuations are natural responses of population to the exogenous and/or endogenous causes: technological changes, school failure, vocational or technical training failure, or noncompatibility of the training with the labor market.
e article is organized as follows.In Section 2, we formulate our model of the labor force evolution.In Section 3, local stability of the positive equilibrium and existence of Hopf bifurcation are given.Numerical simulations are performed, in Section 4, to illustrate our theoretical results.Finally, a summary of the results and some suggestions for future research are presented in Section 5.

Modeling with Differential Equations
For the estimation of labor force evolution, one needs to know the employment rate c, the rate of job loss σ, the mortality rate μ, the maximum population growth rate ρ, the carrying capacity N c , and the time lag needed to contribute in the reproductive process of a new individual looking for a job, τ. erefore, our labor force evolution model is schematized in Figure 1.In this scheme, the active population is divided into unemployed individuals (compartment U) and occupied labor force individuals (compartment L).
e arrows specify the proportion of individuals passing from one compartment to another, and the arrows passing through the small circle indicate the process of reproduction that gives the number of new people looking for a job.e temporal evolution of this scheme is governed by the following differential equations: where the variable L and U can be interpreted as the number of occupied labor force and the number of unemployed, respectively.e initial conditions for system (1) are where e first equation models the evolution of occupied labor force population.In this equation, c is the employment rate.
us, the term cU can be interpreted as the number of new employees among people looking for a job.e term σL corresponds to the number of employees who have lost their job, and the term μL counts the number of deaths in the occupied labor force compartment.
e second equation models the evolution of unemployed population.In this equation, the function ρ(1 − (L + U/N c ))L indicates the number of new people looking for a job.It reflects both maximum growth of job seekers when the unemployment rate is very low and very low labor force growth when the unemployment rate is very high.In other words, population growth is limited by the sharp increase in the number of unemployed.

Equilibria, Local Stability, and Existence of Hopf Bifurcation
3.1.Equilibria.In the following proposition, we prove the existence of the two equilibria.
Proposition 1. System (1) has the following equilibria.(1) has two equilibria: the trivial equilibrium E 0 � (0, 0) and the positive equilibrium Proof.To find equilibria, we consider the following system: It is easy to verify that (0, 0) is a solution of the system (3), thus E 0 � (0, 0) is a trivial equilibrium of the system (1).In the following, we prove that E * is the unique positive equilibrium of the system (1).
(L, U) is a positive solution of the system (3) if and only if It is clear that if (cρ/μ(c + σ + μ)) > 1, then the second equation of ( 4) has a unique positive solution ( e characteristic equation associated to system (5) takes the following form: where To give our results, we recall the Routh-Hurwitz criterion and the result of Kuang [6] for the characteristic equations of the second degree.(1) If there is no such imaginary solution of the characteristic equation ( 6), then the stability of the equilibrium does not change for any τ ≥ 0. (2) If there is only one imaginary solution with positive imaginary part of the characteristic equation ( 6), an unstable equilibrium never becomes stable for any τ ≥ 0. If the equilibrium is asymptotically stable for τ � 0, then there exists a critical value, τ 0 of time delay, such that the equilibrium is uniformly asymptotically stable for τ ≤ τ 0 , and it becomes unstable for τ > τ 0 .(3) If there are two imaginary roots with positive imaginary part, iω + and iω − , such that ω + > ω − > 0, then the stability of the equilibrium can change a finite number of time at most as τ is increased, and eventually it becomes unstable.
According to Lemma 1 (Routh-Hurwitz criterion for a second-degree polynomial) and Lemma 2, we have the following results.Proposition 2. For τ � 0, the equilibrium (L * , U * ) is locally asymptotically stable.

􏽱
}i, with positive imaginary part.Applying Lemma 3.3, we conclude that the statement of our proposition is true.
Next, we establish sufficient conditions for the local existence of Hopf bifurcation at the positive equilibrium.□ Theorem 1. Suppose one of the following situations: Then, a Hopf bifurcation of periodic solutions of system (1) occurs at (L * , U * ) when τ � τ 0 .
Proof.For the proof of this theorem, we verify the conditions of Hopf bifurcation theorem (see, for example [7]).From Proposition 4, the characteristic equation ( 6) has a pair Modelling and Simulation in Engineering of imaginary roots ±iω 0 at τ � τ 0 .It is easy to show that this root is simple.
For the transversality condition, we have After some calculations, we get If one of the two hypotheses is concludes the proof.From eorem 1, the following form of Hopf bifurcation theorem can now be stated (see eorem 2.7., page 291 in [8]).
In view of Figure 2, we observe an oscillatory regime on the labor market: the number of occupied labor force and the number of unemployed oscillate and converge towards equilibrium.is shows that time delay, in the reproduction process, is a key parameter for the existence of a branch of bifurcated periodic solutions from the positive equilibrium and therefore of the existence of oscillations when the delay crosses certain critical values.

Conclusion
In this paper, we presented a study on the dynamics of a model of the evolution of the employed labor force and the unemployed population.We have shown the existence of a positive equilibrium E * , and we have studied the analytic behavior around this point.Specially, we have shown the existence of a branch of periodic solutions bifurcated from the positive equilibrium E * .e existence of the limit cycles, for certain positive values of the parameters of the model, supports the phenomenon of coexistence of the two populations: active occupied and unemployed. is coexistence can be explained by the fact that job seekers occur relatively quickly because of population growth.
As a result of this paper, we plan to study the direction of the bifurcating periodic solutions and resonant codimension-two bifurcation (see, for example [9]).Also, we plan to develop this work by introducing the economic growth effect.e resulting model determines the conditions of interaction between economic growth and labor market evolution.

Data Availability
e parameter values can be found at https://en.wikipedia.org/wiki/Demographics_of_Morocco.

□ 3 . 2 .
Stability and Hopf Bifurcation.By analyzing the characteristic equation associated to (1), we prove the existence of Hopf bifurcation around the equilibrium E * .Consider the change of variables x � L − L * and y � U − U * .en by linearizing system (1) around (L * , U * ), we have dx dt � cy −(σ + μ)x,

Figure 2 :
Figure 2: Examples of times series obtained by numerical integration of labor force evolution model (1).