New nonlinear model of population growth

The model of population growth is revised in this paper. A new model is proposed based on the concept of fractional differentiation that uses the generalized Mittag-Leffler function as kernel of differentiation. The new model includes the choice of sexuality. The existence of unique solution is investigated and numerical solution is provided.


Introduction
Researchers within the field of biology and mathematical biology are interested to know whether or not the certain specie will be instinct or not. This study has fascinated many researchers around the world in recent passed years. For instance to control the spread of a given infectious diseases researchers are interested in their reproductive number, that help to know whether or not the disease will be extinct [1][2][3][4][5][6][7]. However if the model is accurate enough they can give reliable predictions, if the predictions show the extinction of a given species, then laws-makers can take some decisions to protect the specie. We can find many examples of this in developed countries, for instance in South Africa, the government gave strict law against the killing of rhinos. In China, we have the protection of the tigers. A global protection of whale in all oceans and Africans elephants that are nowadays consider as rare species [4][5][6][7]. In case of infectious diseases, the aim is to end the spread of the virus that can considered as a specie, in this case also the control can be done via mathematical predictions. It is therefore important that in both cased the mathematical formulas should be able to portray more accurately the dynamic of the specie in time [8][9]. Generally speaking mathematical models allows a better thoughtful of how the complex interfaces and processes work. Indeed exhibiting of dynamic interactions in nature can provide a manageable way of understanding how numbers alter in excess of time or relation to each other. The aim of this paper is to provide a new model that will be able to describe the population growth more accurately.

Model of population growth
Biological population demonstrating is worried with the changes in populace size and age spreading within a population as a significance of collaborations of creatures with the physical setting, with individuals of their own species, and with bacteria of other kinds. The biosphere is full of interfaces that varies from modest to dynamic. Earth's processes affect human life and a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 are momentously stochastic and seem disordered to naked eye. Nevertheless, a embarrassment of patterns can be perceived and are brought forth by using inhabitants demonstrating as an instrument. Population reproductions are employed to control maximum fruitage for agriculturists, to comprehend the changing aspects of biotic annexations, and have plentiful environmental safeguarding insinuations [8][9]. Thomas Malthus was one of the former to demonstrate that, inhabitants evolved with symmetrical configuration despite the fact that envisioning the providence of humankind [8][9]. Nonetheless, Nurgaliev's law is a mathematical equation that portrays the rate of change of proportions of a population at a given time, in terms of the current population size. It is a deterministic conventional discrepancy equation in which the rate change is articulated as a quadratic function of the population size and this equation is given as: In this equation, the size of a population is denoted by n time is measured in years, a is half of the average probability of a birth of male also for females, of a potential arbitrary parents pair within a is year. b is an average probability of a death of a person within a year. The above model has some limitations, the variation of growth of population is an average between two given times, which is not naturally true because the averaging is not the same at the different interval of time. The second problem is that the model does not take into account the choice of sex and also the memory effect. In this work we shall introduce new parameters to the model and also consider the memory effect induces by the fractional differentiation based on the Mittag-Leffler function.

Fractional differentiation
The topic of fractional differentiation is one of the hot topics nowadays in almost all the fields of science, technology and engineering due to its wide applicability and also its ability of model real world problems more accurately than the classical differentiation. The first definition was proposed by Riemann and Liouville is given below as [10][11][12][13][14] RL D a t ½f ðtÞ ¼ Caputo when working with a real world problem later modified this definition, as he was unable to recover the usual initial conditions, then Caputo modified Eq (2) by putting the derivative inside the integral to obtain: Definition (2) and (3) have been intensively used and misused in the last decades in many fields. However, when looking at the definition, we can see that, they are convolution of functions and the power law, or the kernel of transformation is the well-known power law. However, it is clear when looking at the behaviour of some physical phenomena that, not all of them follow the power law. Recently Caputo and Fabrizio suggested a step a head in fractional differentiation when they replaced the power law with exponential decay law as presented below [10][11][12][13][14] CF 0 D a t ½f ðtÞ ¼ And Goufo and Atangana proposed the modified version in several research papers and it is given as follows [10][11][12][13][14] CFR 0 D a t ½f ðtÞ ¼ But their proposition was rejected due to the criteria that need to be satisfied for an operator to be called fractional derivative. However their idea was great because a new kernel was introduced with no singularity. Atangana and Baleanu, to solve the problem in Caputo and Fabrizio operator, they suggested a new kernel based on the generalized Mittag-Leffler function, that is the more suitable function that was introduced to solve some problems of disc of convergence of power law. The function is also considered as the queen of fractional calculus and is more natural than power law. Their definitions are given below as follow [15][16][17][18][19][20]: This new version was applied in the theory of Chaos with great success therefore in this paper, we make use of this fractional differentiation to provide new model of population growth.

New model of population growth
Many physical observed facts are said to follow the power law evolution. The use of fractional differentiation to predict the population growth was investigated before with Caputo power law fractional derivative in the following [20][21][22][23]. More importantly the expansion of mankind of population growth is obviously one of those one of those. However the chose of power law used to model such dynamical system must be chosen with care. The growth of population does take place exponentially as indicated by several classical models, or this does not take place with the trend of power law of x α (the power law population growth can be observed in less developed countries were the rate of birth is very high). Additionally this does not occur with only a fading local memory (the real world situation for fading memory population can be found in developed countries where the rate of birth is very small as time goes on) as in the process of the diffusion within a porous media but rather combine both fading memory and power law. The only fractional operator that can with care and accurately replicate this dynamical system is perhaps the fractional differentiation with generalized Mittag-Leffler kernel. In this paper to accurately include into mathematical formulation the effect of fading memory and also power law, we convert the classical derivative with Atangana-Baleanu fractional derivative, which takes into account the power-law population growth together with fading local memory population growth. In this section, a new model of population growth is suggested using the concept of fractional, in addition to this, a new parameter taking into account the choice of partner will be introduced to well represent the physical investigation into mathematical formulas. Let assume N(t) to be the size of density population in a given period of time. Let p be the probability of an adult to choose a partner with same sex then the following equation is suggested: The new function v(t) is the selection function that a given individual will be convince to choose a partner with same sex. The new model induces also the memory effect due to the fractional differentiation.
We shall first present the equilibrium points of this dynamical system. To obtain them, we assume that the function is time independent therefore Eq (8) is transformed to Therefore the realistic equilibrium point is when the proportion of the rate of death with the difference between the birth contribution and the factor of choice of partner. However if the following inequality holds a − (1 − p)v < 0 then mankind specie will die out. If a = (1 − p)v then in a near future also mankind will vanish. If the quantity is big enough then mankind will survive.

Existence of solution
The conditions within which the new equation admits a positive solution will be discussed in this section. To do this, we consider X = C[a, b] the Banach space of every continuous real functions defined in the closed set [a, b], which bestowed with the sub norm and Z be the shaft defined as: Z = {N 2 X, N(t) ! 0, a t b}. We shall present the following Banach fixed-point theorem that will be used for the existence demonstration.
Definition 1: Let E be a real Banach space with a cone H. H initiates a restricted order in E in the succeeding approach [18] x y ) y À x 2 H: For every x, y 2 E the order interval is defined as ha, bi = {f 2 E: a f b}. A cone K is denoted normal, if one can find a positive constant j such that h, d 2 K, F < h < d ) khk jkdk, where F denotes the zero element of K.
Theorem 1 [18]: Let H be a closed set subspace of a Banach space of D. let G be a contraction mapping with Lipschitz constant g < 1 from H to H. Thus G possesses a fixed-point t Ã in H. In addition, if t 0 is a random point in H and {t n } is a sequence defined by t n+1 = Gt n (n = 0, 1, 2. . .), then for a large number n, t n tends to t Ã in H and dðt n ; t Ã Þ g n ð1À gÞ dðt 1 ; t 0 Þ. We present also some properties of Atangana-Baleanu derivative in Caputo sense. Theorem 2 [19]: Let f(t) 2H 1 (a, b), b > a such that the Atangana-Baleanu fractional derivative exists, then the following relationship holds: Proof: By definition we establish the above relation (10) using the Laplace transform operator as follow: Applying on both side of Eq (12) the Laplace transform, we obtain the following expression The prove Eq (10b) we use another method that consists of solving the following time fractional ordinary differential equation with Atangana-Baleanu derivative in Caputo sense ABC 0 D a t f ðtÞ ¼ uðtÞ ð15Þ With the aim to find the function f(t), to do this we employ again the Laplace transform on both sides we obtain Rearranging we obtain This completes the proof.
Here applying the AB-fractional integral on Eq (8), we obtain the following We will treat the above inequality piece by piece we first start with the integral part.
Henceforth the mapping G is increasing. By the conjecture, we get Gm ! m, Gn n. Henceforth the operator G: hn, mi ! hn, mi is compact within the framework of lemma 2 and continuous in view of lemma 1. Since H is a normal cone of G.

Uniqueness of solution
In this section, we discuss with care the conditions under which the unicity of the solution are obtained. To establish these conditions, we consider evaluating the following quantity.
Therefore if the following condition holds 1À a ABðaÞ J þ ab a ABðaÞGðaþ1Þ J < 1 then, the mapping G is a contraction, which implies it has a fixed-point, thus, the new model admits a unique positive solution.

Numerical solution via forward-corrector method
The recent development of fractional differentiation based on the Mittag-Leffler function has induced a new type of Volterra fractional differential equations. As presented earlier, the fractional integral calculus associated o the new fractional calculus is the set of functions for which the their fractional integral in Atangana and Baleanu sense is an average of the given function and the Riemann-Liouville fractional integral. This new design has therefore opened way to many new studies, for instance what can we do to solve the Volterra version of a given equation numerically. It is well known that the Corrector method is very accurate method to handle Volterra equations, due to the Riemann-Liouville, however, with the new fractional integral one could possible apply also the corrector method in the integral part and apply another numerical method in the other part. In this section, we introduce the Forward-Corrector method to solve our new model. Nevertheless there are two ways to handle numerically fractional differential equation based on the new fractional differentiation. In our case, we could solve our problem directly in its present form or solve its Volterra version. We shall start with the Volterra version.
Therefore according to [18] the fractional variant of the one step Adam-Moulton method for the second part of our equation is given by: In the second part we use the forward approximation as follows Putting Eqs (29) and (28) into Eq (17), we obtain the following numerical approximation: This approach can be used to solve many other fractional differential equations based on the new fractional differentiation.
The second approach to solve our problem is to discretize the Atangana-Baleanu time fractional derivative. Koca and Atangana suggested the numerical approximation of the new derivative as follow [19]: ð32Þ E a;2 ðzÞ ¼ Replacing the above in Eq (8), using also the forward numerical scheme, then the numerical approximation solution of the new model is given as: ðt n À t kþ1 ÞE a;2 À a 1À a ðt n À t kþ1 Þ À Á À ðt n À t k ÞE a;2 À a 1À a ðt n À t k Þ À Á

Numerical simulations
In this section, we present the numerical replication of the model for different values of fractional order using the proposed numerical scheme. The numerical solutions are depicted in

Conclusion
The aim of this work was to suggest a nonlinear fractional differential equation that could be used to describe the density of population growth taking into account real world behaviors. To do this, we introduced a new component that considers the choice of partner. The analysis of existence of positive solution of the new model was examined via the fixed-point theorem. The new model was solved numerically using the modified approach that fit well the new fractional integral. Some numerical simulations were done as function of fractional order.