Mathematical modelling of unemployment as the effect of COVID-19 pandemic in middle-income countries.

WHO (World Health Organization) has stamped Coronavirus disease 2019 (COVID-19) as a pandemic. This disease has affected most of the population in every aspect. One of those main aspects is the economical burden on the middle class and poor people. To overcome these issues, there should be multi-stage precautions, and awareness was required apart from these, there should be a multi-level posterior effort on facing the challenges that were much needed. In that way, to join the dealing of facing the posterior challenges of this disease, we try to help the society in our way of modelling approach, so that we developed a mathematical model to identify the complexities of the virus and its actual phenomena along with its economical oriented effects and perspectives. During the time of the COVID-19 pandemic situation, millions of employees lost their jobs irrespective of age consideration. More youngsters in many countries have been looking for a new jobs even after facing so much of struggles as the effect of this pandemic. In these facts, a system of ordinary differential equations (ODEs) model is constructed and analysed the outbreak of COVID-19 and its effects on employment. The basic reproduction number R 0 has determined using the next-generation matrix technique. If R 0 ≤ 1 disease-free equilibrium is globally asymptotically stable. If R 0 ≥ 1 , the endemic equilibrium is global asymptotic stability is achieved. To show the analytic findings, numerical simulations were performed. We hope that this work is suitable to control unemployment for all kinds of income populations.

Introduction

In pandemic situations from Dec 2019 to Dec 2021 worldwide, many developed and developing countries face many ways of struggle came to control it. However, not clear till now, so many difficulties economically and socially. We took the issue of unemployment in middle-income countries. Many of the country execute the given categories of lockdown for the response to COVID-19 (e.g., school, business, and public park closures, as well as mobility limitations). Following the lockdown phase, additional focused efforts like contact tracing and isolation of possibly affected persons would be implemented, aided by testing.

Unemployment rates for young employees, which were already high, suddenly increased during the pandemic situation of COVID-19. The total unemployment rate among young working-age groups 16–24 moved from 8.4 to 24.4% the year from 2019 to 2020. Unemployment among those aged 25 and up increased from 2.8 to 11.3%. In the year 2020, the rate of unemployment had even higher for young workers (Black at 29.6%, Hispanic at 27.5%, and Asian American/Pacific Islander (AAPI) at 29.7%). The financial results of the COVID-19 economic system on younger employees may also persist for years [1].

This unique aspect of the COVID-19 crisis makes it easy to distinguish layoffs from job loss by looking at the impact of activity status on mental health. The main difference is that it reveals whether future employment aspirations give psychological protection against the loss of current wages and activities. Mental health and nervous system problems are the third leading cause of illness burden in South Africa, behind HIV and other infectious diseases [2].

The effect of COVID-19 cases and fatalities are backbone of deflation on economic and social changes in European countries and gives much impact on European economies also on the pandemic situation. The descriptive data suggest that COVID-19 casing and casualties are more common in Italy and the United Kingdom than in other European countries [4]. Furthermore, the progress of creating new vacancies takes time that is outstanding to the continual growth of the population, which is a slim chance that the number of jobless people will alter when the current vacancies are. In this instance, too, freshly created vacancies will be insufficient, resulting in an unmanageable situation. As a result, new opportunities must have produced as quickly as possible, as our findings suggest that delaying the creation of new vacancies destabilises the system [5].

Personal protection, therapy when individuals are diagnosed by treatment when the individuals are late-diagnosed, effective spraying of the environment, and cleaning probably contaminated surfaces are all part of the control set that can assist lower the virus’s number [6]. The SARS-CoV-2 virus is changed with time, which is the standard characteristic of any virus. Generally, when a virus transforms into a mutant, its traits remain the same or have a minimum influence [7]. The comparative analysis of the COVID-19 for pandemic situations has concentrated on the first and second waves. It describes the growth of the virus, death rates, and rate of recovery of the people. To suggest control of the virus takes prevention actions for COVID-19. The statistical tool ANOVA is used to compare the fractal interpolation function [8].

India has implemented public health and social measures based on local epidemiology and risk assessments; similarly, a full-day lockdown, limited mobility, and educational facility closures. In terms of the total infected population, it has a lower COVID-19 death rate [9]. Controlling individuals from travelling and rigidly implementing personal precautions at the same time can considerably minimise the COVID-19 outbreak. The overall infected population plot compares to the actual infected population of India’s COVID-19 outbreak’s second wave [10]. In addition, people throughout the country have been suffering from a shortage of medical resources. Mathematically, the unique coronavirus dynamics have been modelled to analyse its dynamics. However, the fractional-order COVID-19 model with control techniques is the subject of this paper. The COVID-19 model investigates using the Caputo fractional derivative, besides quarantine and hospitalised compartments [11].

This work is developing as follows: Section (2) describes the mathematical model and analysis that comprises the subsections as model assumption and determination, equilibria, stability analysis, and the global stability of our model. Section (3) describes the numerical results of the simulation of and the parameter effect. It simulates the design in Python to show the figures of the deterministic model and project the real-world data from WHO/UN WOMEN/ UNDP and Economic Policy Institute collections. At the end of the section (4), given our result discussion and conclusion.

Mathematical models: formulation and analysis

Model assumption with determination

The total population separates into four subclasses. Based on what we are considered the poor people. Let P(t), I(t), U(t), and R(t) be the fraction of the poor population, the COVID- 19 infected individuals, unemployment population, and the recovered population, respectively, at time t.

This model includes the following assumptions: The flow of people from one class to the next has displayed in the diagram below:

Let be the inflow rate of individuals in a class, also show the natural death rate per capita in each division. The COVID-19-infected peoples increase through touch already affected by COVID-19 people, with poor populations are taking standard mass action incidence.

Let be the rate of transmission of COVID-19 infection, so that denotes the infected people’s contact with uninfected people and they are travel without any protection the individual.

We further assume that poor people are quitting individual quarantine at a rate .

be the inflow of COVID-19 severely infected death values.

is the rate of recovery from COVID-19 infection. COVID-19 infects people who have been taking individual quarantine or some medical treatment due to recovering from the COVID-19 infection.

rate of people the got jobless/unemployment. During the COVID-19 lockdown period, many people have stayed at home. Therefore, they have no job due to the social distancing and people communication stoping purpose closed all the Corporates/Industry/Hotels/Tourism. Daily labour is also affected in this way.

rate of people goes to recover. Some of the Corporate supported the employees working from home at the time few industries also like followed this way and created new vacancies to give opportunities for jobless people, those who got the job would move to R.

Diagram of the model

Hence, we developed the following model:where with as initial conditions.

As N is constant, we use following transformation to get:where

with as initial conditions. Here, note that .

Equation (2) has following two equilibriums:

(i) Trivial equilibrium

(ii) Endemic equilibrium = and is the basic reproduction number of system (2).

Disease-free steady state exists always but endemic equilibrium exists only if . For finding the Basic reproduction number, we use the next-generation matrix method. Using this method, we obtain two non-negative matrices and evaluated at , such thatFor definition and evaluation of matrices, and refer to [15]. Now, the spectral radius of (), which corresponds to basic reproduction number, is

Stability analysis

Stability analysis of disease-free equilibrium

Theorem 1

Disease-free steady state is locally asymptotically stable for .

The matrix of model Eq. (2)The eigenvalues of this matrix are given by the following roots:We can show that all eigenvalues of are negative. The disease-free equilibrium is locally asymptotically stable if [14].

Stability analysis of endemic equilibrium

Theorem 2

Endemic equilibrium , whenever it exists, is locally asymptotically stable when [15, 16].

Proof

The matrix of the system 2One eigenvalue of this matrix is - and the other three eigenvalues are the roots of the following cubic equation:whereBy Routh–Hurwitz criteria, roots of the cubic equation will have negative real parts if . Hence, the equilibrium point is locally asymptotically stable provided that .

Global stability

Global stability of disease-free equilibrium

Theorem 3

The disease-free equilibrium is globally asymptotically stable when .

Consider Lyapunov function asThe derivative of V along solutions of system (2) is given byIf , then . Using Lyapunov LaSalle theorem [3], we conclude that is globally asymptotically stable for .

Global stability of endemic equilibrium

Theorem 4

The endemic equilibrium is globally asymptotically stable when .

we consider the following Lyapunov function:This is a positive definite function in 1 and its derivative along solutions of system (2) is given byClearly, if , . Using Lyapunov LaSalle theorem, we conclude that is globally asymptotically stable whenever [17].

Numerical results and discussion

We used numerical simulation to verify our analytical conclusions under this section. Outcomes that validate are presented here with rationale.

The system (2) is simulated for various set of parameters satisfying the condition of local asymptotic stability of equilibria . The fourth-order Runge–Kutta technique is used to numerically integrate the dynamic behaviour of the model for reaching equilibrium of the system. We took the parameter set , and the parameters have been considered to have a reasonable value.

We evaluated the equilibrium stability with respect to the initial conditions to make sure it was stable (Fig. 2a–e). Figure 2a, b draws that the stability of endemic equilibrium point to different initial values on poor population and COVID-19-infective populations. And how they come to move other can stables at a point. Figure 2c, d depicts the stability of endemic equilibrium point for different initial values on both poor people and unemployment population. Figure 2e depicts the stability of endemic equilibrium point for different initial values on all compartment dynamics behaviour at time t.

Fig. 2

Stability analysis of Eq. (2) and 3D Plot of

Figure 2f, g explores the relation between with other parameters. In Fig. 2f, we to look about the relations between and . In Fig. 2g, we look about the combined effect of , and .

The pandemic situation has low-, middle-, and high-income country; poor populations have unemployment during the lockdown periods. In this aspect embrace, a mathematical model finds in the data. We are trying to our results along with data. In this view, we have absorbed the data on the COVID-19 situation.

In Fig. 4a, b, the rate of unemployment has been assimilated by a pooled average of April, May, and June data in the years 2019 and 2020. The category of racial and ethnic is entirely incompatible. Hispanic refers to Hispanic/Latinx of any race through the White, Black, and AAPI refers to non-Hispanic white, non-Hispanic Blacks, and non-Hispanic Asian American/Pacific Islanders, in orderly [18] [Source:  Economic Policy Institute Current Population Survey. http://microdata.org].

Fig. 4

Unemployment rate in the spring of 2019 and 2020 young workers in the COVID-19 labour market

In 2020, 8.8% of worldwide working hours will be lost, equating to 255 million full-time employment, compared to the fourth quarter of 2019. Even though overall working-hour losses are largish reduction. It indicates the fact that adjustments occurred through reduced per week workdays in the lodging and dining services sector fell by more than 20% [12].

India’s unemployment rate was over 8% in December 2021, up from 7% the previous month. While the unemployment rate had decreased dramatically since peaking in April 2020, the outbreak of new coronavirus strains, combined with periodic lockdowns, resulted in a shifting pattern of unemployment gripping the country in 2021 [13].

Conclusion

Globally, the impact of COVID-19 is well noticed the numbers of the terms of death or terms of people who got unemployed. Particularity, this affected much in middle-income countries. We may consider the middle-income countries from Asia and Europe. We tried to convey that unemployment reduces if we reduce the effectiveness of COVID-19 spread.

Our numerical simulation indicates controlling the number of new confirmed cases and infections by lowering the transmissive. We also suggest that the SOPs for controlling infections results were helping the jobless people expected much reduced/decreased in the future.

The regain of the parameter shows that people who get treatment due to the COVID-19 infection will have been recovering. We have shown this fact in numerical simulation also.

The parameter (infected to unemployed ) is showed some increase due to the effect of lockdown and its consequences. Figure 3c, d shows the (poor people infected rate) which is compared for = 0.5, 0.3 decreasing the unemployment rate, and also decreasing at the same time the peoples are moving from one place to another place certainly.

Fig. 3

Plot in how unemployment rate decreasing and parameter effects

Therefore, they will have new job recruitment are available at that time to improve the lifestyle of jobless to employment the way new vacancies originated. During the COVID-19 pandemic, governments and officials attempted to alleviate unemployment, and the labour market improved by expanding job openings. Over the COVID-19 recession, governments also revised their job reallocation strategies. There is one proposal from general employees to the government that the government has to support an “unemployment insurance policy’ which may introduce which might result in a dramatic drop in labour demand [4].

Table 1

Description of parameters

Parameter	Description
Λ	Recruitment rate
μ	Natural death rate
μ_1	Death rate due to COVID-19 infection and with out immunity
β	Interaction rate of infected and poor people compartments
γ	Rate of progression from COVID-19 infection to recovery
α	Rate of people moves to unemployment
δ	Flow rate of infected to unemployment
σ	Rate at which the population moves to recovery from unemployment

Table 2

The unemployment rate during the COVID-19 pandemic time from Jan 2020 to Dec 2021

Month →
Year ↓	Jan	Feb	Mar	Apr	May	Jun	Jul	Aug	Sep	Oct	Nov	Dec
2020 (%)	7.22	7.76	8.75	23.52	21.73	10.18	7.4	8.35	6.68	7.02	6.5	9.06
2021 (%)	6.53	6.89	6.5	7.97	11.84	9.17	6.96	8.32	6.86	7.74	6.97	7.91