Two efficient methods for solving fractional Lane–Emden equations with conformable fractional derivative

In this paper, we introduce two reliable efficient approximate methods for solving a class of fractional Lane–Emden equations with conformable fractional derivative (CL-M) which are the so-called conformable Homotopy–Adomian decomposition method (CH-A) and conformable residual power series method (CRP). Furthermore, the proposed methods express the solutions of the non-linear cases of the CL-M in terms of fractional convergent series in which its components can be computed in an easy manner. Finally, the results are given by graphs for each case of the CL-M at different values of α in order to demonstrate its accuracy, applicability, and efficiency.


Introduction
The subject of fractional calculus can be considered as a generalization of the classical integer order calculus (derivative and integration), and due its important role in science and engineering, it has been gaining considerable attention of many authors and researchers [1][2][3][4][5]. The reader who is interested with the field of fractional calculus knows that there are many definitions of a fractional order derivative and they have been used to describe several real-life problems in many fields of sciences; the most important and famous ones are Riemann-Liouville, Grunwald-Letnikov, and Caputo derivatives [1,2,[6][7][8][9][10]. Riemann-Liouville and Caputo use the integral in its construction, particularly the Cauchy integral formula with modifications. Hence, occasionally, we need complex computations to get the fractional derivative in the sense of Riemann-Liouville and Caputo. Furthermore, they do not satisfy the non-linear derivative rules as product, quotient, and chain rules. More recently, a new definition of fractional derivatives which is the so-called conformable fractional derivative has been introduced and attracted the attention of many researchers because it contains many characteristics that correspond to the usual derivative, particularly the Leibniz rules [11][12][13][14][15][16]. This definition is very simple and more welcome than other fractional definitions since it has been receiving a lot (2020) 28:42 Page 2 of 11 of attention, many applications and phenomena can be modeled based on the CFDs, and it contains many interesting advantages such as the following: it is a local derivative that simulates the normal derivative because it depends on the limit in its formulation and it generalizes all concepts of ordinary calculus and can solve different fractional differential equations with all cases. In addition to this definition, there is another type of local derivatives called non-conformable fractional derivative, and for this purpose, the authors point out the publications [17][18][19][20][21][22]. In recent years, many authors have handled and studied the Lane-Emden equations because they were used to formulate lots of phenomena in physics and astrophysics [23][24][25][26][27][28]. The aim of this paper is to find the approximate solution of fractional Lane-Emden equations with conformable fractional derivative (CL-M) using the conformable Homotopy-Adomian decomposition method (CH-A) and conformable residual power series method (CRP); both of these methods are effective and easy to use for solving nonlinear CL-M, without linearization or discretization. The benefit of these techniques over the other methods is that they can be performed directly to the given problem by choosing an appropriate value for the initial guess approximations, and they also reduce the difficulty appearing in the computation of the complicated terms [29][30][31][32][33][34][35][36].
This paper is organized as follows: In the "Preliminaries" section, a preliminary introduction of the conformable fractional order derivative is presented. In "The CH-A for solving CL-M" section, we introduced the CH-A for solving CL-M. The basic idea of the CRP technique for solving CL-M was given in the "A CRP for solving CL-M" section. Finally, a conclusion has been drawn.
Definition 1 Given a function y :[ 0, +∞) → R, then the CFD of order α of y is given by:

Definition 2
The conformable integral of order α is defined by: where the integral is the usual Riemann improper integral and α ∈ (0, 1).

Lemma 1
Let α ∈ (0, 1] and f , g be α-differentiable at a point t > 0, then: The CH-A for solving CL-M In this section, we present the CH-A for solving CL-M; the CH-A is a hybrid method in which it combines between the HAM and the Adomian decomposition method (ADM). The calculations involved in this technique are more easy than the standard HAM especially when the non-linear term in the CL-M is decomposed by using the Adomian polynomials. To start the procedure, let us consider the general form of CL-M of order α > 0: where x > 0 and 0 < α ≤ 1.
To solve this problem using CH-A, first, rewrite Eq. (3), as follows: Then, integrate Eq. (4) twice with respect to x, so the general fractional solution of Eq. (3) is given by: where c 1 and c 2 are constants and d α x = x 1−α dx. In order to solve Eq. (5) by means of HAM, we need to seek the auxiliary linear operator: We now define the non-linear operator as: Consequently, the mth-order (m ≥ 1) deformation equations can be expressed using the Adomian polynomials as: where Hence y m = χ m y m−1 + HR m ( y m−1 ) Starting with an initial approximation: So, we have: where m=0 is the set of Adomian polynomials of f (y) which is defined as follows: Finally, the exact solution of Eq. (3) can be calculated by: and the nth-order approximate solution of Eq. (3) is given by: Now, we examine some known and charming cases that have been resulting from Eq. (3) as follows: Eq.(3), then the following equation for x > 0 is called: The CL-M of the first kind.

Case 2 Consider f (y) = e y as in Eq. (3), then the following equation for x > 0 is called:
The CL-M of the second kind.

A CRP for solving CL-M
The implementation of the CRP will be described in this segment in order to get the fractional power series solution of the CL-M represented by Eq. (3) subject to the initial conditions: The solution of Eq. (3) with respect to Eq. (24) is proposed by CRP as a fractional power series expansion about the initial point x = x 0 as follows: where y m is expressed by y m = c m  By truncating the summation into Eq. (27), we get the kth approximate series: For the convergence of the fractional power series Eq. (28), we advise to see [41]. For simplification and before employing the CRP, Eq. (3) can be written in the form: Now, to evaluate c m that appears in the series expansion given by Eq. (28), substituting y k (x) into Eq. (29), we obtain the following kth residual function: Furthermore, CD (k−1)α Res k (x 0 ) = 0, k = 1, 2, 3, ..., k which is coming from lemma 1 (property 3) and this is a basic rule to compute the coefficients c m , m = 2, 3, ..., k.
However, finding c m demands to solve the algebraic expression: In this manner, all the desired coefficients c m will be found. For solving Eq. (3) using CRP, different cases for Eq. (3) will be discussed as follows: Case (i) Consider Eq. (17), subject to the initial conditions y(0) = 1 and y (0) = 0, according to the initial conditions given above, and by using Eq. (30), hence, we have c 0 = 1, c 1 = 0, c 2 = −1 3 , c m = 0 for all m ≥ 3, substituting these values into Eq. (28) then the general fractional series solution in this case can be obtained by: The nature of CRP solution of Eq. (17) Also, the nature of the CRP solution of Eq. (21) is presented in Fig. 3. Case (iii) Consider the problem given by Eq. (23) with respect to the initial conditions y(0) = 0 and y (0) = 0, then after taking the initial conditions in consideration and by employing Eq. (30), we get c 0 = c 1 = 0, c 2 = −1 3 , c 3 = 0, c 4 = 1 5 , c 5 = 0, ...; thus, the fractional power series solution in this case is: The quality of the CRP solution of Eq. (23) is given in Fig. 4.

Conclusions
In this paper, we applied CH-A and CRP to find the approximate analytic solution of some classes of CL-M in terms of infinite fractional power series. The proposed methods introduced an easy way to compute the components of the solution that have been converging rapidly to the exact solution. The results obtained by CH-A and CRP prove that these algorithms are highly effective and convenient in non-linear cases of the CL-M and can be employed to examine a wide class of non-linear fractional mathematical models.