key: cord-0706272-ccjoq3ed
authors: Kouidere, Abdelfatah; Khajji, Bouchaib; Balatif, Omar; Rachik, Mostafa
title: A multi-age mathematical modeling of the dynamics of population diabetics with effect of lifestyle using optimal control
date: 2021-01-11
journal: J Appl Math Comput
DOI: 10.1007/s12190-020-01474-w
sha: 5dff42fcdfd914c6d3371afbf7755eebc050cb21
doc_id: 706272
cord_uid: ccjoq3ed

Diabetes is a disease which caused by socio-environmental and / or genetic factors. The negative effect of socio-environmental or lifestyle leads a susceptible individual to become a diabetic. On the one hand, social interaction wields a great deal of influence over lifestyle. On the other hand, genetic factors are the main cause of the birth diabetes genetic disorder. Considering these above mentioned factors. In the present paper, we study a discrete age continuous mathematical model that describes the dynamics of diabetics. We highlight the negative impact of socio-environmental on diabetic patients according to age groups. We also suggest an optimal strategy to implement the best campaigns of rising awareness that aims at protecting diabetic patients from the negative impact of a lifestyle that leads them to complications. In addition to psychological treatment and follow-up of diabetic patients with complications, an awareness campaign will also be carried out for people with potential diabetes that aims at educating them about the dangerous of diabetes and its complications. Pontryagin’s maximum principle is used to characterize the optimal controls and the optimality system is solved by an iterative method. The numerical simulation is carried out using MATLAB.

On November 14, 2019, The world celebrated International Diabetes Day on which the World Health Organization (WHO) [1] and the International Diabetes Federation (IDF) [2] have chosen the theme "Family and diabetes" with the slogan "Diabetes is important for the whole family ". The International Diabetes Day has an objective the raising of awareness about this disease and its complications by focusing on prevention, adoption of a healthy lifestyle, early examination, especially for people who are subject to diabetes.

According to the International Diabetes Federation [2] , statistics indicate that more than 463 million people with diabetes in 2019 and that 50 percent of people suffer from diabetes are not aware of being a subject to disease. And these numbers increase exponetially every year. The number of people with diabetes is expected to reach 578 million in 2030 and 700 million in 2045.

Diabetes increases the possibility to be affected by other serious illnesses, especially respiratory illnesses like the SARS-CoV-2 virus. According to the American Diabetes Association (ADA) [42] , people with diabetes are at higher risk for viral infections. According to the Scientific Federation of Diabetics [33] , In response to the current COVID-19 pandemic, governments of many countries have restricted the movements of their citizens, confining them to their family environment. facilities such as, gymnasiums, sports centers and swimming pools are closed. Therefore, the inactivity and the lack of exercising have caused many people to gain excessive weight which will put them at a high of developing diabetes [34] [35] [36] [37] [38] [39] [40] [41] .

Diabetes can be caused by two main factors, genetics and/ or lifestyle. Lifestyle manifested in malnutrition which stem from poverty, and social factors that can result from work place and family problems. These factors inter alia may lead to diabetic complications such as, kidney failure, vision loss and others.

In recent years, a large amount of research has been conducted in relation to diabetes in general and in mathematics particularly. These mathematical models aim at understanding and analyzing the dynamics of diabetics. In a related research work, in 2004 Boutayeb et al. [3] introduced introduced the diabetes complication (DC) model to find out how many changes of diabetics without complications (D) and diabetics with complications (C). In the DC model the number of new patients (incidence) of diabetes is assumed to be constant. In 2007 Derouich et al. [4] proposed a A population model of diabetes and prediabetes using a system of ordinary differential equations. In 2014 Derouich et al. [4] proposed also an optimal control approach modeling the evolution from pre-diabetes to diabetes with and without complications. They showed the existence of an optimal control and then used a numerical implicit finite-difference method to monitor the size of population in each compartment. Their model shows that, using optimal control, the number of diabetics with and without complications can be significantly reduced in a period of 10 years. In 2018 Permatasar et al. [8] proposed Mathematical model to elaborate the prevalence of diabetics has been determined by diabetes complication (DC) model. In the DC model, people with diabetes were classified into two compartments, uncomplicated diabetics (D) and complicated diabetics (C). Diabetes is known as a disease caused by lifestyle and genetic factors. A bad lifestyle leads a susceptible individual to become a diabetic. Bad lifestyle is strongly influenced by risky social interaction. In the other side, a genetic factor is the main cause of the diabetes genetic disorder birth. Consider these both factors, the DC model was developed into a susceptible diabetes complication (SDC) model. Susceptible individuals were involved in the calculation of risky interactions. The SDC model is a first order nonlinear differential equation. The number and the change of individuals in each compartment can be determined from the solution of this model. In this paper, the SDC model is applied to predict changes of diabetics prevalence in the United States. As a result, the SDC model is good enough to predict the prevalence, and Al Helal et al. [32] Insulin injections and exercise scheduling for individuals with diabetes: an optimal control model. In 2019 Kouidere et al. [5] introduced a Discrete Time to the Dynamics of a Population of Diabetics with Highlighting the Impact of Living Environment, and Anarina et al. [18, 19] studied the dynamics of glucose, insulin, and free fatty acids with time delay: the impact of bariatric surgery on type 2 diabetes mellitus. In 2020 Kouidere et al. [6] In order to have realistic model. They proposed to study an optimal control approach with delay in state and control variables in a discrete mathematical model of kouidere et al. [5] , that time with delay represent the measuring the extent of interaction with the means of treatment or awareness campaigns. Also, many researches have focused on this topic and other related topics, and Sweatman et al. [14] introduce a mathematical model of diabetes and dietrelated lipid metabolism, leptin sensitivity, insulin sensitivity and VLDLTG clearance predict pathways to health and type II diabetes, and also Auni Aslah Mat Daud et al. [46] proposed a mathematical model of the population dynamics of DM during pregnancy is developed and analyzed. Four independent variables have been considered, namely the numbers of nonpregnancy nondiabetic women, diabetic nonpregnant women, diabetic pregnant women and diabetic pregnant women with complication. The model is described by a system of ordinary differential equations. The stability of the equilibrium point is analyzed using Routh-Hurwitz criteria. The model is numerically simulated using MATLAB to verify the analytical results. The model has only one nonnegative equilibrium point which is asymptotically stable. The equilibrium solution is further investigated using simple sensitivity analysis. The results of simple sensitivity analysis of the equilibrium solution suggest the key parameters of the model. The equilibrium point of the model indicates the influential parameters that can be controlled to address the issue. Also, many researches have focused on this topic and other related topics [7] [8] [9] 16, [20] [21] [22] .

According th World Health Organization (WHO) [43] , the number of incidences of diabetes is not constant and tends to increase every year. The increase is due to the lifestyle of the world's population. One of the factors influencing lifestyle is social interaction. Social interaction is a significant factor affecting lifestyle changes so as to increase the potential of a healthy-susceptible individual into diabetics [44] . In addition, the incidence of diabetes is also often associated with genetic factors from parents who have a history of diabetes [45] . However, they did not take into account the negative impact of the socio-environmental factors according to age groups, that is the aim of this paper.

The increase in the number of diabetics is related to several factors. Diabetes is caused by unhealthy lifestyles such as, the lack of physical activity, unhealthy eating patterns, and other unhealthy habits. Diabetics who have experienced complications are already aware of their poor heath conditions. In contrast, diabetics without complications were unable to lead a healthy lifestyle, especially in undiagnosed cases. The influence of society on diabetes varies vis-à-vis their ages. the closer people are in age, the more influencial they become towards one another. In addition to the foregoing, the aim of this paper is to shed light on socio-environmental such as behavioral, social, and economic factors; these factors can be controlled in contrast to genetic factors that are difficult to control because diabetes cannot be cured, but one can control the rate of glucose in the blood by virtue of a balanced diet. The extent to which people who do not suffer from diabeties negatively impact diabetics when both groups belong to the same age group. The extent to which people who do not suffer from diabeties negatively impact diabetics when both groups do not belong to different age group.

As previously indicated in discrete mathematical modeling of Kouidere et al., diabetics were divided into two main groups, namely: diabetics due to genetic factors P, diabetics due to behavioral and economic factors in which the entourage plays a crutial role E. furthermore, diabetics without complications D and diabetics with complications C. can be added to the mix. And that was elaborated in detail in the methodology section.

Besides the aforementioned works, we will study a mathematical modeling of the dynamics of a population of diabetics which contains the following additions: * A discrete-age continuous mathematical modeling, * Dividing age into multiple age intervals.

We note, as mentioned above, that most researches about diabetes and its complications focusing on continuous and discrete time models and describing differential equations. Recently, more and more attention has been paid to the study of optimal control. (see [12, [24] [25] [26] [27] [28] [29] [30] [31] …and references cited There).

In Sect. 2 of this paper, we presented the P E DC continuous mathematical model and we gave some basic properties, and we presented the optimal control problem for the proposed model where we gave some results concerning the existence of the optimal controls and we characterized these optimal controls using Pontryagin's maximum principle. In Sect. 3, we proposed mathematical model with multiple age and optimal control where we gave some result concerning the existence of the optimal control and we characterized the optimal controls. We gave numerical simulations, in Sect. 4. Finally, we concluded the paper in Sect. 5.

We considered a mathematical model P E DC that described the dynamics of a population of diabetics with highlighting the effect of lifestyle. by divided the population N into four compartments: P is Number of pre-diabetic peopole through genetics factors (Genetic predisposition), E is represented the Pre-diabetics due to the negative impact of socio-environmental factors on diabetics ( without Genetic predisposition), D is represented the number of diabetics without complications, C is represented the number of diabetics with complications (Complication as Cardiovascular disease, Foot damage, neuropathy, nephropathy and retinopathy). The total population of individuals, N (t) at time t is given as:

The graphical representation of the proposed model is shown in Fig. 1 .

Hence, we presented the P E DC diabetic mathematical model by the following system of differential equations:

with P(0) ≥ 0, D(0) ≥ 0, E(0) ≥ 0, C(0) ≥ 0 are the given initial states.

μ: Natural death rate (Natural mortality is not due to diabetes and its complications). -I 1 : Denote the incidence of pre-diabetes due to genetics factors. -I 2 : Denote the incidence of pre-diabetes due to lifestyle (living environment).

β 1 : The rate of pre-diabetics who has been diabetics without complications.

β 2 : The rate of diabetics whose complications are cured.

β 3 : The rate of diabetics become diabetics with complication through the sudden sock. -γ : The rate of diabetics people probability of developing diabetes through the negative effect of socio-environmental factors on diabetics. -δ: Mortality rate due to complications.

α: The rate of negative effect of socio-environmental factors E(t) on diabetics people D(t)( According to Hill et al. [47] , the interaction between individuals without diabetes E(t) who have unhealthy lifestyles and diabetics without complications D(t) leads to negative health outcomes for those with diabetes).

Proof We have according to the first equation of the system (1) that

The both sides of the last inequality are multiplied by exp (

Integrating this inequality from 0 to t gives:

Then

Similarly, we prove that E(t), D(t) and C(t) are positive for all t ≥ 0.

Proof By adding the equations of system (1), we obtain

where I = I 1 + I 2 that I is the sum of the recruitment rates I 1 and I 2 , and N (0) represents the initial values of the total population. Thus, lim t→∞ sup N (t) = I μ . It implies that the region Ω is a positively invariant set for system ( 1) . So, we only need to consider the dynamics of the system on the set Ω.

The system (1) can be rewritten in the following form:

where

The second term on the right-hand side of (2) satisfies: 2 , and M is some positive constants, independent of the state variables P (t) , E (t) , D (t) and C(t).

then

Thus, it follows that the function ϕ is uniformly Lipschitz continuous. the restriction on P (t) ≥ 0 , E (t) ≥ 0, D (t) ≥ 0 and C(t) ≥ 0, we see that a solution of the system (1) exists [17] .

The strategy of control that has been adopted consists of an awareness program to minimize the negative effect of behavioral factors on diabetics without complication. Our main aim is to minimize the number of people evolving from the stage of pre-diabetes to the stages of diabetes with and without complications. In this model, we included three controls u(t), v(t) and w(t) for t ∈ [0, T ]. The first control is symbolized by u represented treatment, Its role is manifested in treating the complications of diabetes, and overcoming them, which gives an opportunity to control the level of glucose in the blood at a given time t. The second control v represented the awareness program through media and education, by raising awareness of the seriousness of the negative impact of behavioral factors on diabetics without complication D. The third control w represented the awareness program, through raising awareness of non-diabetics about the negative imapct of behavioral, economic and social factors that lead to diabetes at a given time t.

Our aim is to minimize the number of diabetics without complication individuals and the cost of vaccination program. Mathematically, it can be interpreted by optimisation of the objective functional

In other words, we seek the optimals controls u * , v * and w * such that

where U is the set of admissible controls defined by

Theorem 4 Consider the control problem with system (3). There exists an optimal

Proof The existence of the optimal control can be obtained using a result by Fleming and Rishel [16] , checking the following steps :

-The set of controls and corresponding state variables is nonempty. To prove this condition we use a simplified version of an existence result of Boyce and DiPrima ( [15] , Theorem 7.1.1) To prove that the set of controls and the corresponding state variables is nonempty, we will use a simplified version of an existence result [15] . Let

where X 1 , X 2 , X 3 and X 4 form the right-hand side of the system of equations (10) . Let u, v and w for some constants and since all parameters are constants and X 1 , X 2 , X 3 and X 4 are continuous, then F P , F E , F D and F C are also continuous. Additionally, the partial derivatives

. . , 4 are all continuous. therefore, there exists a unique solution (P, E, D, C) that satisfies the initial conditions. therefore, the set of controls and the corresponding state variables is nonempty and condition 1 is satisfied in U .

T ]} is convex and closed by definition (see appendix 1). -All the right hand sides of equations of system are continuous, bounded above by a sum of bounded control and state, and can be written as a linear function of u with coefficients depending on time and state. -The integrand in the objective functional

The state variables being bounded, let

B, ζ 4 = G and ζ = 2 then it follows that:

Then from Fleming and Rishel [16] , we conclude that there exists an optimal control.

In order to derive the necessary condition for optimal control, the Pontryagins maximum principle [13] , given in was used. This principle converts into a problem of minimizing a Hamiltonian H (t) at time t defined by

where f i is the right side of the differential equation of the ith state variable at time t.

Theorem 5 Given the optimals controls (u * , v * , w * ) and the solutions P * , E * , D * and C * of the corresponding state system (3), there exists adjoint variables λ 1 (t), λ 2 (t),λ 3 (t) and λ 4 (t) satisfying:

With the transversality conditions at time T :

Furthermore, for t ∈ [0, T ], the optimal controls u * , v * and w * are given by

Proof We use the Pontryagins maximum principle [13] , for characterized the optimal controls. So, we defined the Hamiltonian H as follows:

At this case it is considered that N is constant. For t ∈ [0, T ], the adjoint equations and transversality conditions can be obtained by using Pontryagin's maximum principle [5, 10, 11, 13, 23] such that

With the transversality conditions [48] at time T : For, t ∈ [0, T ] the optimal controls u, v and w can be solved from the optimality condition

We have

By the bounds in U of the controls, it is easy to obtain u * , v * and w * are given by (16, 17, 18) in the form of system(3).

In this section, the emphasis is on the impact of behavioral factors that can result from the influence that non-diabetics wield over diabetics without complications by age groups. The negative effect of a group of children under the age of 20 on each other behaviorally, or a group between the age of 18 and 50 years on each other, as well as a group over the age of 50. It will also highlight how different age groups affect others, just like adults and young adults. As we divided these age groups k into three main groups: The first category k = 1 is less than 20 years old, the second category k = 2 is between 20 years to 50 years, and the third category k = 3 is over 50 years old.

The negative effect of no diabetics on diabetics within the same age Diabetes is a chronic and deadly disease, it has its negative effects on the individual and society. As previously mentioned, the individual develops diabetes by genetic, behavioral, economic and social factors, according to IDF [2] In this section, we will highlight the negative effect of behavioral factors directly caused by people without diabetes. Influence of age groups of the same age as diabetics 

N . These age groups will have a significant influence on the increase in complications of diabetics for several reasons, three of which are having the same age, having the same way of thinking and lacking the awarness of the dangers that diabetes may present.

In this paragraph, we want to highlight the negative impact of behavioral factors of people who do not suffer from diabetes on diabetics from another age group

for k = 1, 2, 3. Adults with children Fig. 2 .

The model is presented with age group k. Hence, we present the diabetic model by the following system of differential equations: 

The both sides in last inequality are multiplied by exp (

.P k (t)) ≥ 0 Integrating this inequality from 0 to t gives:

Similarly, we prove that E k (t) ≥ 0, D k (t) ≥ 0 and C k (t) ≥ 0.

Proof By collecting system equations (10), we obtain

where I = I 1 + I 2 , and N (0) represents the initial values of the total population. Thus, lim t→∞ sup N (t) = I μ . It implies that the region Ω is a positively invariant set for system ( 10) . So, we only need to consider the dynamics of the system on the set Ω.

so the system (10) can be rewritten in the following form:

where

The second term on the right-hand side of (11) satisfies

where M are some positive constants, independent of the state variables P k (t) , E k (t) , D k (t) and C k (t).

Then ϕ(X 1 ) − ϕ(X 2 ) ≤ V . X 1 − X 2 Thus, it follows that the function ϕ is uniformly Lipschitz continuous. the restriction on P k (t) ≥ 0 , E k (t) ≥ 0, D k (t) ≥ 0 and C k (t) ≥ 0, we see that a solution of the system (11) exists [17] .

We will follow the same strategy above to optimize mathematical modeling with multi category of age. We included three controls u(t), v(t) and w(t) for t ∈ [0, T ], represented awarness program and treatment for diabetics people and sensitization about the negative impact of behavioral, economic and social factors.

The problem is to minimize the objective functional (13) where A > 0, B > 0 and G > 0 are the cost coefficients.they are selected to weigh the relative importance of u(t), v(t) and w(t) at time t. T is the final time.

In other words, we seek the optimals controls u * ,v * and w * such that

where U is the set of admissible controls defined by

Theorem 9 Consider the control problem with system (12) . There exists an optimal control (u * , v * , w * ) ∈ U 3 such that J (u * , v * , w * ) = min

Proof The existence of the optimal control can be obtained using a result by Fleming and Rishel [16] , checking the following steps:

-The set of controls and corresponding state variables is nonempty. To prove this condition we use a simplified version of an existence result of Boyce and DiPrima ( [15] , Theorem 7.

T ]} is convex and closed by definition. -All the right hand sides of equations of system are continuous, bounded above by a sum of bounded control and state, and can be written as a linear function of u, v and w with coefficients depending on time and state.

is clearly convex on U It rest to show that there exists constants ζ 1 , ζ 2 , ζ 3 , ζ 4 > 0, and ζ such that

The state variables being bounded, let

A, ζ 3 = B, ζ 4 = G and ζ = 2 then it follows that:

Then from Fleming and Rishel [16] we conclude that there exists an optimal control.

In order to derive the necessary condition for optimal control, the pontryagins maximum principle [13] , given in was used. This principle converts into a problem of minimizing a Hamiltonian H (t) at time t defined by

where f i is the right side of the differential equation of the ith state variable at time t. Theorem 10 Given the optimals controls (u * , v * , w * ) and the solutions P * k , E * k , D * and C * k of the corresponding state system (12) , there exists adjoint variables λ 1 (t), λ 2 (t), λ 3 (t) and λ 4 (t) satisfying:

With the transversality conditions at time T : λ 1 (T ) = 0, λ 2 (T ) = 0, λ 3 (T ) = 1 and λ 4 (T ) = −1.

Furthermore, for t ∈ [0, T ], the optimals controls u * , v * and w * are given by

Proof We use the Pontryagins maximum principle [13] , for characterized the optimal controls. So, we defined the Hamiltonian H as follows:

For t ∈ [0, T ], the adjoint equations and transversality conditions can be obtained by using Pontryagin's maximum principle [5, 10, 11, 13, 50] such that

For, t ∈ [0, T ] the optimal controls u, v and w can be solved from the optimality condition

By the bounds in U of the controls, it is easy to obtain u * , v * and w * are given by (16)- (18) in the form of system. 

In this section, we investigate and compare numerical results of the following control strategies when applied for reducing the the negative effect of lifestyle on diabetics and miminizing the number complicated diabetics. Te strategies are (i) Strategy 1: treatment with education and awareness program.

(ii) Strategy 2: the awareness program through media and education, by raising awareness of the seriousness of the negative impact of behavioral factors on diabetics without complication.

As explained in [51] [52] [53] [54] , the adjoint system is solved by using the backward in time finite-difference method with terminal conditions λ k (T ) = 0, where T = 120 months and initial conditions P(0) = 6,660,000, D(0) = 10,200,000, E(0) = 10,000,000, C(0) = 5,500,000. The controls u, v and w are considered to be bounded and the weights in the objective functional are estimated to be A = 100, B = 100 and G = 100.. With initial value of controls and the initial condition X (0) = X 0 , the state solutions system is solved forward in time using the finite-difference method of order four. The update of the controls is done using a convex combination of the current and previous controls to obtain the new solution for X and λ. The method continues by using these new updates aiming at finding a fixed point (X , λ, u ) that u = (u, v, w). This iterative process terminates when the last and preceding iterations are negligible close and the last iteration is the solution of the optimal problem. The parameter values in the state system and the objective function are obtained from different literatures and others are estimated depending on the dynamic of diabetic population as shown in Table 1 . In this paper all plots for state variables are in the logarithmic form.

Different simulations can be performed using different parameter values in Table 1 taken from [5] :

Proceeding from Fig. 3 , Table 2 gave the evolution of the number of diabetics without and with complications without control after 120 months.

We note in Fig. 3 and Tables 2 and 1 that the number of diabetics without complications decreased from 10.2 × 10 6 to 4.83 × 10 6 and that the number of diabetics with complications increased from 5.5×10 6 to 14.03 × 10 7 . These changes occured due to several reasons, including the sudden development of diabetes, lack of taking the necessary and appropriate precautions to reduce the disease, in addition to the effect of the behavioral factors of improper eating and unorganized diet and also due to family and work problems. Ordinary people 1.5 × 10 7

Proceeding from Fig. 4 , Table 3 gave the evolution of the number of diabetics without and with complications with control after 120 months. Based on Fig. 4 and Table 3 , we notice an increase in the rate of diabetics without complications from 10.2 × 10 6 to 10.9 × 10 6 . The number of diabetics with complication increased from 5.5×10 6 to 7.24 × 10 6 , and also increase of rate E from 10 6 to 10.9 × 10 6 And by following the medical advice (as ACE inhibitors, acetaminophen, aspirin or ibuprofen [49] ) and the psychological follow-up of diabetics, as happened in many awareness campaigns mediated by the behavior factors of diabetics to avoid the mentioned problems, whether family or private work as possible. 

Proceeding from Fig. 5 , Table 4 gave the evolution of the number of diabetics with multi age without control after 120 months. We noted from Fig. 5 and Table 4 , we wanted to highlight the behavioral factors effecting diabetics when they are with both same and different age groups, so that we can notice the negative impact of these groups which led to a significant decrease in the number of diabetics without complications and at the same time an significant increase in the number of diabetics with complications, and this decrease can be interpreted by the fact that these groups do not take into account the dangers that diabetics can be exposed to by means of their ignorance of the aftermaths that their careless behaviors may have on diabetics.

Proceeding from Fig. 6 , Table 5 gave the evolution of the number of diabetics with multi age without control after 120 months.

In Fig. 6 and Table 5 , compared to Fig. 4 , we observe a slight decrease in the number of diabetics with complications and a noticeable rise in the number of diabetics without complications after making a diagnosis on the reasons leading to this dangerous development. This decrease occured due to awareness campaigns by both the state and specialized federations and NGOs, where all age groups are targeted and sensitized about the risks of the negative impact of the socio-environmental factors on diabetics according to their age. Those campaigns introduced the risks and complications of this disease so as to reduce the development of the disease from dangerous levels that are difficult to control, through a healthy diet, and keeping people as far away as possible from family and work problems as mentioned previously in Fig. 4 .

We can also merge multiple assemblies as (u(t), v(t)) and (u(t), v(t), w(t)) thus get a variety of results.

In this work, we studied a discrete age continuous mathematical model that describes the dynamics of diabetics population. We highlighted the negative impact of the people's lifestyle and the socio-environmental where diabetics get influenced badly by non-diabetics. This negative impact has a direct consequence on the health of diabetics without complications which eventually leads to a rise in the forms of complications. The seriousness of complications can lead directly to death. We also divided diabetes patients into two categories: the first got the disease due to genetic factors and the second due to the negative impact of lifestyle. We have also divided the age groups into three groups; below 18 years old, between 18 and 50 years old, and over 50 years old. The goal is to find out how strong is influence of each group on the other. We have also suggested three optimal strategies to reduce the increase in the number of diabetic patients with complications and to reduce the negative impact of lifestyle on diabetic patients. We presented a new application of the theory of optimal control to reduce the impact of social and economic behavior of age groups on the health of diabetics through systematic differential equations, awareness campaigns as well as psychological treatment and follow-up. We applied the results of the control theory and we managed to obtain the characterizations of the optimal controls. The numerical simulation of the obtained results showed the effectiveness of the proposed control strategies. Fig. 7 The optimal controls u * , v * and w * v * = min v max , max v min , (λ 3 (t) − λ 4 (t))

w * = min w max , max w min , (λ 2 (t) − λ 3 (t))

So, Using these standard optimality arguments, we characterize the control u * (t), v * (t) and w * (t) by:

World Health Organisation: Definition and diagnosis of diabetes mellitus and intermediate hyperglycemia

International Diabetes Federation (IDF). IDF DIABETES ATLAS

A population model of diabetes and prediabetes

Optimal control approach to the dynamics of a population of diabetics

Optimal control strategy for a discrete time to the dynamics of a population of diabetics with highlighting the impact of living environment

Optimal control strategy with multi-delay in state and control variables of a discrete mathematical modeling for the dynamics of diabetic population

A simulation model for the dynamics of a population of diabetics with and without complications using optimal control

A mathematical model for the epidemiology of diabetes mellitus with lifestyle and genetic factors

Mathematical modeling, analysis, and optimal control of abstinence behavior of registration on the electoral lists

On the analysis of a multi-regions discrete SIR epidemic model: an optimal control approach

Using age structure for a multi-stage optimal control model with random switching time

Mathematical Theory of Optimal Processes

Mathematical model of diabetes and lipid metabolism linked to diet, leptin sensitivity, insulin sensitivity and VLDLTG clearance predicts paths to health and type II diabetes

Elementary Differential Equations and Boundary Value Problems

Deterministic and Stochastic Optimal Control

Ordinary Differential Equations

Modeling, the dynamics of glucose, insulin, and free fatty acids with time delay: the impact of bariatric surgery on type 2 diabetes mellitus

A mathematical model of type 1 diabetes involving leptin effects on glucose metabolism

A multi-region discrete time mathematical modeling of the dynamics of Covid-19 virus propagation using optimal control

A new approach for determining multi-objective optimal control of semilinear parabolic problems

A new mathematical modeling with optimal control strategy for the dynamics of population diabetics and its complications with effect of living environment

A mathematical modeling with optimal control strategy of transmission of COVID-19 pandemic virus

A multi-regions discrete-time epidemic model with a travel-blocking vicinity optimal control approach on patches

A multi-regions SEIRS discrete epidemic model with a travel-blocking vicinity optimal control approach on cells

A multi-regions SIS discrete influenza pandemic model with a travel-blocking vicinity optimal control approach on cells

A new epidemic modeling approach: multi-regions discrete-time model with travel-blocking vicinity optimal control strategy

Effect of awareness programs and travel-blocking operations in the control of HIV/AIDS outbreaks: a multi-domains SIR model

On the regional control of a reaction-diffusion system SIR

Age-structured mathematical modeling approach to short-term diffusion of electronic commerce in Spain

A nonlinear dynamic age-structured model of e-commerce in spain: Stability analysis of the equilibrium by delay and stochastic perturbations

Insulin injections and exercise scheduling for individuals with diabetes: an optimal control model

International Diabetes Federation (IDF): IDF diabetes atlas

COVID-19. A challenge for healthcare professionals

COVID-19, an emerging coronavirus infection: current scenario and recent developments-an overview

COVID-19 and diabetes: knowledge in progress

COVID-19 and diabetes management: what should be considered?

COVID-19 infection in Italian people with diabetes: lessons learned for our future (an experience to be used)

Internationat Federation Diabetes: Diabetes Voice "COVID-19 and diabetes

COVID-19 and diabetes: is metformin a friend or foe

Report of a WHO/IDF Consultation

Development and analysis of a mathematical model for the population dynamics of Diabetes Mellitus during pregnancy

Optimal Control Applied to Biological Models. Chapmal Hall/CRC

Optimal control problem for an electoral behavior model

Mathematical modeling and optimal control of malaria

Optimal control techniques on a mathematical model for the dynamics of tungiasis in a community

Optimal Control Applied to Biological Models

Sustainability and optimal control of an exploited prey predator system through provision of alternative food to predator

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

The authors thank the editor and the anonymous reviewers for very helpful suggestions and comments that helped us to improve the paper.

Author contributions All authors contributed equally to the writing of this manuscript. All authors read and approved the final version.Funding There is no funding source for the publication of this paper.

Competing interests All the authors have no competing interests.

The

Proof Take any controls u, v ∈ U and λ ∈ [0, 1]. then 0 ≤ λu

The controls u, v and w in Fig. 7 were utilized to reduce the increase in the number of diabetic with complications by means of treatment and awareness programs. In fact diabetes is incurable, but it can be controlled. We noted that the controls curves were remained at the same level. for t ∈ [0, T ], the optimal controls u * , v * and w * are given by u * = min u min , max u min , (λ 3 (t) − λ 4 (t))