Modified Auxiliary Equation Method versus Three Nonlinear Fractional Biological Models in Present Explicit Wave Solutions

In this article, we present a modified auxiliary equation method. We harness this modification in three fundamental models in the biological branch of science. These models are the biological population model, equal width model and modified equal width equation. The three models represent the population density occurring as a result of population supply, a lengthy wave propagating in the positive x-direction, and the simulation of one-dimensional wave propagation in nonlinear media with dispersion processes, respectively. We discuss these models in nonlinear fractional partial differential equation formulas. We used the conformable derivative properties to convert them into nonlinear ordinary differential equations with integer order. After adapting, we applied our new modification to these models to obtain solitary solutions of them. We obtained many novel solutions of these models, which serve to understand more about their properties. All obtained solutions were verified by putting them back into the original equations via computer software such as Maple, Mathematica, and Matlab.


Introduction
Since the emergence of humanity, people have taken a great interest in understanding natural phenomena, beginning from fire, lightning, thunder, earthquakes, and volcanoes, ranging all to way to the nano-particle.In the beginning of the 18th century, many scientists demonstrated an interest in studying these phenomena.The development of these studies continued until the middle of the 19th century.Partial differential equations (PDEs) became an essential tool for studying other branches of science, and many phenomena have been explained using PDEs.Many techniques for solving PDEs have been developed by scientists and researchers in their efforts to explain these phenomena.Many methods have been derived to obtain exact and approximate solutions of these models.The most important results in determining explicit solutions of nonlinear partial differential equations (NLPDEs) were derived in [1].These types of methods are best known as continuous symmetry transformation groups [2][3][4][5][6].At present, we use computer software (e.g., Maple, Mathematica) to make this kind of transformation.Since the 20th century, the investigation of the partial differential equation has become an independent field.Partial differential equations have been applied to the study of many phenomena in different fields, such as sound, heat, electrostatics, electrodynamics, fluid dynamics, elasticity, quantum mechanics, solid state physics, fluid mechanics, hydrodynamics, optics, plasma physics, chemical kinetics, biological phenomena, and more.It is known that there exist various analytic solution methods for NPDEs [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21].However, in the general case, there is no central theory for NPDEs.There is no unified method that can be applied to all types of nonlinear PDEs.For these reasons, we have to discover new methods for solving these kinds of models, and cannot just make modifications to methods that were discovered more than 100 years ago.Recently, fractional nonlinear partial differential equations have played an important role in the representation of many phenomena.Fractional rule half differential equations are increasingly chronic in accordance with mannequin problems between fluid flow, pay up and vile areas over the application.Fractional derivatives grant spiffy arms for the story on devotion or ancestral houses concerning various substances then processes.Half-order derivatives and integrals prove according to lie more useful for the method concerning sure electrochemical problems than the first-rate models.Fractional differentiation and integration operators are additionally chronic for extensions about the pervasion then suspense equations.Many researchers tried in conformity with finding out extra properties expecting that form of derivatives.They discovered partially about methods in accordance with changing the nonlinear fractional half differential equations among everyday differential equation together with integer order.For more important points in relation to these types of derivatives, we refer the reader to References [22][23][24][25][26][27][28][29].In this paper, we use conformable derivative properties and apply them to three nonlinear biological models.These models are a fractional biological population model [30-34], a fractional equal width equation [35][36][37][38], and a fractional modified equal width equation [39-42], respectively.
The rest of this research paper is organized as follows: In Section 2, we give the headlines of the modified auxiliary equation method.In Section 3, we employ a novel computation to get the solitary solutions of the fractional biological population model, fractional equal width equation, and fractional modified equal width equation.In Section 4, we study the obtained solutions and their novelty with the previous method.In Section 5, we represent the conclusion of this study.

Fundamental Steps of the New Technique
Consider a nonlinear partial differential equation (NLPDE) with the following form: where F is a polynomial in u(x, t) and its partial derivatives in which the highest-order derivatives and nonlinear terms are involved.In the following, we give the steps of this method: Step 1.Using the wave transformation: Using this transformation on (1), we convert the nonlinear partial differential equation into an ordinary differential equation with the following form: Step 2. Suppose that the solutions of ODE (2) have the following form: where (a i , a) are arbitrary constants that will be determined later, while f (θ) satisfies the following ODE: Step 3. Determine the positive integer N in Equation (3) by balancing the highest-order derivatives and the nonlinear terms.
Step 4. Substituting Equations (3) and (4) into Equation ( 2) and collecting all the terms of the same power (K i f (θ) ) where (i = −N, . . ., N) and equating them to zero, we get a system of algebraic equations that can be solved by Maple or Mathematica to get the values of a i , b i , and (α, β, σ).
Step 5. Substituting these values and the solutions of Equation (4) into Equation (3), we obtain the exact solutions of Equation (1).

Applications
In this part, we apply a novel computational method to two fractional models.
• Fractional biological population model: This model describes population dynamics.It also gives a simple example of how complex interactions and processes work.The model has the following form: where u impersonates the population density and λ(u 2 − s) constitutes the population changes because of deaths and births.

•
Fractional equal width equation: This model is usually used to describe complex physical phenomena in various fields, and has the following formula: where [h, r] are arbitrary constants.

•
Fractional modified equal width equation: This model refers to the replica of one-dimensional wave propagation in nonlinear form with dispersion processes, and has the following formula: where [h, r] are arbitrary constants.
Applying the conformable derivative definition and its properties to Equations ( 5) and (7) in respective order [u(x, we transform the fractional PDE into an integer order ODE in the following order: We balanced the terms in Equations ( 8)-(10) to get the balance value of each of them, obtaining the value of balance equal to one in all of them.According to the novel method, the general solution of Equations ( 8) and (10) is in the following form: while the general solution of Equation ( 9) is in the next form:

Fractional Equal Width Model
Substitute Equation ( 12) and its derivatives into Equation ( 9).Collect all coefficients of the same terms for K f (θ) .Equating them to zero, we obtain a system of algebraic equations.Solving this system via any computer software (e.g., Maple, Mathematica), we obtain: where Consequently, the solitary wave solutions of population density have the following formulae: When β 2 − 4 α σ < 0 and σ = 0, we get: When β 2 − 4 α σ > 0 and σ = 0, we get: Figures 1-3 to give more explanations about these solutions.The accuracy of all obtained solutions was examined by putting them into their original equations using Mathematica 9.