Polynomial Least Squares Method for Fractional Lane – Emden Equations

This paper applies the Polynomial Least Squares Method (PLSM) to the case of fractional Lane-Emden differential equations. PLSM offers an analytical approximate polynomial solution in a straightforward way. A comparison with previously obtained results proves how accurate the method is.


Introduction
The equation analyzed in this article was published at the end of the 19th century by Jonathan Homer Lane in [1] and at the beginning of the 20th century explored in detail by Robert Emden in [2].In the decades that followed, the equation Lane-Emden raised the interest of many researchers who used different methods to determine numerical or analytical solution for the equation.
We remark that finding accurate approximate solutions to the problem (1) and ( 2) is usually not a simple task due to the presence of a singularity in 0.
Depending on the values of the constants and functions involved in (1) and ( 2), there are several particular types of equations with important practical applications, e.g., thermionic currents, gravitational potential of the degenerate white-dwarf stars or isothermal gas spheres.Many studies have been devoted to finding solutions for it, such as [12][13][14][15][16].
This paper has the following structure: The first section will introduce the Polynomial Least Squares Method (PLSM) [17] which permits determination of analytical approximate polynomial solutions for problems of the type (1) and (2).In the second section we will compare approximate solutions obtained by using PLSM with corresponding approximate solutions obtained in previous studies by means of other methods.

The Polynomial Least Squares Method
Attached to the problem (1) and ( 2), we have in view the operator: We will compute approximate polynomial solutions ỹ(x) of ( 1) and ( 2) on the interval [0, 1], satisfying the conditions: Definition 1.An -approximate polynomial solution of the problem (1) and ( 2) is an approximate polynomial solution ỹ(x) satisfying the relations (4) and (5).

Definition 2.
A weak -approximate polynomial solution of the problem (1) and ( 2) is an approximate polynomial solution ỹ(x) satisfying the relation: together with the initial conditions (5).

Definition 3.
We have in view the polynomials One calls the sequence of polynomials P n (x) convergent to the solution of the problem (1) and Theorem 1.The problem (1) and (2) admits a sequence of weak approximate polynomial solutions.
Proof.We compute a weak -approximate polynomial solution, in the sense of the Definition 2, of the type: where d 0 , d 1 , • • • , d n are constants calculated as follows: (1) We substitute the approximate solution (7) in the Equation ( 1) and obtain the remainder: (2) We attach to problem (1) and ( 2) the real functional: Theorem 2. The sequence of polynomials T n (x) from (10) satisfies the property: Moreover, ∀ > 0, ∃n o ∈ N, n > n 0 such that ∀n ∈ N, n > n 0 it follows that T n (x) is a weak -approximate polynomial solution of the problem (1) and (2).
Proof.Taking into account the way the coefficients of polynomial T n (x) are computed and also the relations ( 8)- (10), the following inequalities are satisfied: It follows that: We obtain: From this we obtain that (∀) > 0, ∃n o ∈ N such that (∀)n ∈ N, n > n 0 .It results that T n (x) is a weak -approximate polynomial solution of the problem (1) and (2).

Remark 1.
As far as the above remark is concerned, in order to find -approximate polynomial solutions of the problem (1) and (2) by using the (PLSM), we will first determine weak approximate polynomial solutions, ỹ(x) following the previously described steps 1 to 4. If |D( ỹ(x))| < , then ỹ(x) is also an -approximate polynomial solution.

Error Estimation
We denote with L the following operator: thus, Equation (1) becomes: where x > 0, y(x) : [0, 1] → R, 1 < α ≤ 2, 0 < β ≤ 1, together with the conditions: With ỹ an approximate solution for the problem ( 15)-( 17), using (for simplicity) the notation: , which means as ỹ satisfies: with initial conditions: We define the error function in the following way: ẽ(x) = y(x) − ỹ(x), where y(x) is the exact solution for the problem (1) and ( 2) and we obtain the differential equation for error function: with the conditions: The problem for the error function becomes so : L( ẽ(x)) = −R(x) or: with ẽ(0) = 0, ẽ (0) = 0. Solving the (23) equation in the same manner as described above, we obtain the approximation ẽ(x), we will be able to determine the absolute maximum error: In this manner we estimate the error without knowing the exact solution of the initial problem (1) and (2).

Application 1
Consider the fractional Lane-Emden equation [3]: β = 1 2 together with the initial conditions: The exact solution of the problem (25) and ( 26) is: An approximate solution of this problem using the collocation method was proposed in [3] and the solution presented absolute errors larger than 10 −8 .
Using the Polynomial Least Squares Method (PLSM) presented in the previous section, we choose an approximate solution ỹ(x) of the type: By using the boundary conditions we compute c0 = 0, c1 = 0 and the approximate solution becomes: ỹ The corresponding functional ( 9) is: In order to find the minimum of this functional we can compute the stationary points by equating to zero its partial derivatives with respect to d 2 and d 3 .We obtain as expected d 2 = −1, d 3 = 1 and it is easy to show (by means of the second derivative) that this stationary point is indeed the minimum.It follows that by using PLSM we are able to find the exact solution of the problem, ỹ(x) = x 3 − x 2 .

Application 2
We consider the fractional Lane-Emden Equation [3]: , α = 3 2 , β = 1 together with the initial conditions: The exact solution of the problem ( 27) and ( 28) is: y Again, an approximate solution of this problem using the collocation method was proposed in [3] and the solution presented absolute errors larger than 10 −8 .
Applying PLSM and using the same steps as in the previous example, we are again able to find the exact solution of the problem: ỹ(x) = −x 3 + x 2 .
The exact solution of the problem is y(x) = x 2 cos(x) ( [5]).
For this problem, by using PLSM we obtain the approximate analytical solution: The absolute error of the approximation, computed as the absolute value of the difference between the exact solution and the approximate one, is presented in Figure 1.

Conclusions
The Polynomial Least Squares Method (PLSM) is considered as a simple, efficient, accurate method for calculating approximate polynomial solutions for Lane-Emden-type fractional differential equations.
The comparison with anterior results highlights the accuracy of the method.At the same time, the fact that the solutions are polynomials of relatively small degree leads to solutions that are not only precise but also that present a very simple expression.
Since the method does not actually depend on a particular equation, it can be easily applied to other types of equations, even strongly nonlinear ones.

Figure 1 .
Figure 1.The absolute maximum error using Polynomial Least Squares Method (PLSM) for Application 3.