Next Article in Journal
A Novel Decomposition-Based Multi-Objective Evolutionary Algorithm with Dual-Population and Adaptive Weight Strategy
Previous Article in Journal
The Frobenius Number for Jacobsthal Triples Associated with Number of Solutions
Previous Article in Special Issue
Multidimensional Markovian BSDEs with Jumps and Continuous Generators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spectral Treatment of High-Order Emden–Fowler Equations Based on Modified Chebyshev Polynomials

by
Waleed Mohamed Abd-Elhameed
1,*,
Mohamed Salem Al-Harbi
2,
Amr Kamel Amin
3 and
Hany M. Ahmed
4
1
Department of Mathematics, Faculty of Science, Cairo University, Giza 12613, Egypt
2
Department of Mathematics, College of Science, University of Jeddah, Jeddah 23218, Saudi Arabia
3
Department of Basic Sciences, Adham University College, Umm AL-Qura University, Makkah 28653, Saudi Arabia
4
Department of Mathematics, Faculty of Technology and Education, Helwan University, Cairo 12612, Egypt
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(2), 99; https://doi.org/10.3390/axioms12020099
Submission received: 19 November 2022 / Revised: 7 January 2023 / Accepted: 9 January 2023 / Published: 17 January 2023
(This article belongs to the Special Issue Higher Order Differential Equations)

Abstract

:
This paper is devoted to proposing numerical algorithms based on the use of the tau and collocation procedures, two widely used spectral approaches for the numerical treatment of the initial high-order linear and non-linear equations of the singular type, especially those of the high-order Emden–Fowler type. The class of modified Chebyshev polynomials of the third-kind is constructed. This class of polynomials generalizes the class of the third-kind Chebyshev polynomials. A new formula that expresses the first-order derivative of the modified Chebyshev polynomials in terms of their original modified polynomials is established. The establishment of this essential formula is based on reducing a certain terminating hypergeometric function of the type 5 F 4 ( 1 ) . The development of our suggested numerical algorithms begins with the extraction of a new operational derivative matrix from this derivative formula. Expansion’s convergence study is performed in detail. Some illustrative examples of linear and non-linear Emden–Flower-type equations of different orders are displayed. Our proposed algorithms are compared with some other methods in the literature. This confirms the accuracy and high efficiency of our presented algorithms.

1. Introduction

Numerous studies have focused on the solutions of the Lane–Emden high-order initial value problems, which can be found in many fields, including astrophysics and stellar structure [1,2,3]. Multiple techniques were utilized for solving such problems, including the Haar wavelet collocation method [4], the Legendre wavelets spectral method [5], the Chebyshev operational matrix algorithm [6], the Galerkin operational matrix method [7], the Adomian decomposition method [8] and its modified method [9,10] and the cubic B-spline method [11]. See, for instance, [12,13,14,15,16] for some other articles exploring various Lane–Emden–Fowler equations.
The field of numerical analysis benefits greatly from the use of spectral methods. They play important roles in a number of applied fields, including fluid dynamics and engineering (see, for instance [17,18,19,20,21]). These methods provide accurate solutions to several types of differential and integral equations. The core idea behind the employment of these methods is based on the selection of suitable linear combinations of different special functions that are often orthogonal polynomials (see, for example, [22,23,24]).
Spectral methods can be divided into three main categories. The methods are collocation, tau, and Galerkin. The residual must be orthogonal to the test basis functions in the Galerkin technique in order for it to work. We select trial basis functions that satisfy the specified initial or boundary conditions. This method was successfully utilized for several differential equations [25,26,27,28]. The tau method is more flexible than the Galerkin method because there is no need to choose the basis functions such that they satisfy the given underlying conditions, and it also enforces the residual to be orthogonal to a set of test functions. These test functions may be chosen differently from the set of trial functions; see, for example, [29,30,31]. Out of all the spectral techniques, the collocation method is by far the most widely used. This is due to its capability to handle various differential equations regulated by various initial/boundary conditions [32,33,34,35].
Chebyshev polynomials have essential parts in the scope of special functions and approximation theory. These polynomials are essential for numerical analysis in general and for solving differential equations numerically in particular. The four families of Chebyshev polynomials, which are themselves special polynomials of Jacobi polynomials, are widely used in numerous contributions to the field of numerical analysis. For example, in [36], the authors employed the first kind of Chebyshev polynomials together with the radial basis functions to handle the space-time fractional Klein-Gordon equation. The orthogonal projection method is applied in [37] based on utilizing the Chebyshev polynomials of the second kind. The Chebyshev polynomials of the third kind are employed in [38] to treat some types of integral and integro-differential equations. Recently, modified Chebyshev polynomials of the second kind were introduced and utilized in [6] to treat a certain type of Emden–Fowler type equation. These polynomials are constructed in the sense that they satisfy the underlying initial and boundary conditions.
In the field of numerical analysis, the use of the operational matrices of derivatives and integrals has been highlighted. This utilization provides efficient algorithms to obtain accurate approximate solutions of various types of differential and integral equations using flexible computations. The idea to extract an operational matrix of derivatives is based on selecting suitable basis functions in terms of celebrated special functions and expressing the first derivative of these basis functions in terms of their original types (see for instance, [39,40]).
In this research, four techniques for handling singular high-order type equations—in particular, Emden–Fowler equations regulated by initial conditions—are presented. These algorithms are based on the use of appropriate spectral approaches. Building two types of operational matrices of derivatives of the third kind of the shifted Chebyshev polynomials and a type of modified one is the main idea for presenting the suggested methods.
The manuscript is structured as shown below. The shifted Chebyshev polynomials of the third kind and their modified counterparts are discussed in Section 2. Section 3 is limited to developing a new operational matrix of modified third-kind Chebyshev polynomials’ derivatives. A new Galerkin operational matrix of derivatives is used in Section 4 to handle the singular high-order type equations regulated by initial conditions. The use of tau and collocation methods to solve two numerical approaches for singular high-order Emden–Fowler-type equations governed by initial value problems is examined in Section 5. The convergence study of the two suggested expansions is addressed in Section 6. The examples in Section 7 are all corroborated by comparisons with several other methods from the literature. Finally, Section 8 presents our findings and discussions.

2. An Overview on the Shifted Third-Kind Chebyshev Polynomials and Their Modified Types

This section concentrates on presenting some elementary properties of the third-kind Chebyshev polynomials and their shifted types. Furthermore, new kinds of orthogonal polynomials, which we call “modified third-kind Chebyshev polynomials” will be introduced.

2.1. An Overview on the Third-Kind Chebyshev Polynomials

The orthogonal Chebyshev polynomials of the third kind V j ( x ) , of degree j, can be defined on [ 1 , 1 ] as
V j ( x ) = cos j + 1 2 θ cos θ 2 , where x = cos ( θ ) , θ [ 0 , π ] .
The recurrence relation shown below can be used to build these polynomials
V j ( x ) = 2 x V j 1 ( x ) V j 2 ( x ) , j 2 ,
with the initial values:
V 0 ( x ) = 1 , V 1 ( x ) = 2 x 1 .
It is easy to extend the polynomials that are defined in (1) to be defined on a general interval [ a , b ] .
The orthogonal Chebyshev polynomials of the third kind can be defined on [ a , b ] as
V j * ( x ) = cos j + 1 2 ϑ cos ϑ 2 , where cos ( ϑ ) = 2 x a b b a , ϑ [ 0 , π ] .
So, in this respect, we have
V j * ( x ) = V j 2 x a b b a , x [ a , b ] .
The orthogonality relation for these polynomials is
a b w 0 ( x ) V i * ( x ) V j * ( x ) d x = π ( b a ) 2 , i = j , 0 , i j ,
where w 0 ( x ) = x a b x .

2.2. Introducing Modified Third-Kind Chebyshev Polynomials

For our current goals, it is quite helpful to present the family of polynomials ϕ n , j ( x ) j 0 defined as:
ϕ n , j ( x ) = ( x a ) n V j * ( x ) .
The main characteristic of the new introduced polynomials ϕ n , j ( x ) is that they fulfil the n initial conditions
D r ϕ n , j ( x ) x = a = 0 , r = 0 , 1 , , n 1 .
Remark 1.
Since, for n = 0 , we have
ϕ 0 , j ( x ) = V j * ( x ) ,
so, the polynomials ϕ n , j ( x ) are generalizations of the orthogonal polynomials V j * ( x ) .
According to the weight function: w n ( x ) = 1 ( x a ) 2 n x a b x , the polynomials ϕ n , j ( x ) are orthogonal on [ a , b ] in the sense that:
a b w n ( x ) ϕ n , i ( x ) ϕ n , j ( x ) d x = π ( b a ) 2 , i = j , 0 , i j .

3. Operational Matrix of Derivatives of the Modified Third-Kind Chebyshev Polynomials

In this section, we will establish the operational matrix of derivatives of the third-kind modified Chebyshev polynomials ϕ n , i ( x ) , n = 0 , 1 , 2 , . This will lead to the development of a new formula that expresses the first-order derivative of ϕ n , i ( x ) in terms of the polynomials themselves. The desired operational derivatives matrix can be determined from this formula. It is necessary to start with the following lemma.
Lemma 1.
For all non-negative integers p and i, the following reduction formula is valid:
5 F 4 p , i p , i p + 1 2 , 2 i p + 1 , i + n p + 1 2 i 2 p , i p + 1 , i p + 3 2 , i + n p 1 = ( 1 ) p ( 2 i 2 p + 1 ) ! 2 ( 1 + 2 i ) ( i + n p ) ( 2 i p ) ! ( 2 i p + 2 n ( 1 + p ) ) p ! , p e v e n , ( 1 2 n ) ( p + 1 ) ! , p o d d .
Proof. 
First, if we assume that
S p , i , n = 5 F 4 p , i p , i p + 1 2 , 2 i p + 1 , i + n p + 1 2 i 2 p , i p + 1 , i p + 3 2 , i + n p 1 ,
then due to Zeilberger’s algorithm (see, [41]), it can be demonstrated that S p , i , n given by (9) satisfies the following recurrence relation of order four:
( 2 i p + 1 ) ( 2 i p + 2 ) ( 2 i p + 3 ) ( i + n p + 3 ) ( p 2 ) ( p 1 ) p S p 3 , i , n 2 ( 2 i 2 p + 7 ) ( i p + 3 ) ( 2 i p + 1 ) ( 2 i p + 2 ) ( i + n p + 2 ) ( p 1 ) p S p 2 , i , n 4 ( 2 i 2 p + 5 ) ( 2 i 2 p + 7 ) ( i p + 2 ) ( i p + 3 ) ( 2 i p + 1 ) ( i + n p + 1 ) p S p 1 , i , n + 8 ( 2 i 2 p + 3 ) ( 2 i 2 p + 5 ) ( 2 i 2 p + 7 ) ( i p + 1 ) ( i p + 1 ) ( i p + 3 ) ( i + n p ) S p , i , n = 0 ,
with the following initials:
S 0 , i , n = 1 , S 1 , i , n = 2 n 1 ( 2 i + 1 ) ( i + n 1 ) , S 2 , i , n = i + 3 n 1 ( i 1 ) ( 2 i + 1 ) ( i + n 2 ) , S 3 , i , n = 6 ( 2 n 1 ) ( i 2 ) ( 2 i 3 ) ( 2 i + 1 ) ( i + n 3 ) .
One can exactly solve the recurrence relation (10) to obtain:
S p , i , n = ( 1 ) p ( 2 i 2 p + 1 ) ! 2 ( 1 + 2 i ) ( i + n p ) ( 2 i p ) ! ( 2 i p + 2 n ( 1 + p ) ) p ! , p e v e n , ( 1 2 n ) ( p + 1 ) ! , p o d d .
This proves Lemma 1. □
Theorem 1.
Let i be any positive integer. In terms of their original polynomials, the first-order derivative of the polynomials ϕ n , i ( x ) can be stated directly as
D ϕ n , i ( x ) = j = 0 i 1 θ i , j ( n ) ϕ n , j ( x ) + ϵ n , i ( x ) ,
where the coefficients θ i , j ( n ) and ϵ n , i ( x ) are given in the form of
θ i , j ( n ) = 2 ( 2 n 1 ) ( j i ) b a , i > j , ( i + j ) e v e n , 2 ( ( 2 n + 1 ) i ( 2 n 1 ) j + 1 ) b a , i > j , ( i + j ) o d d , 0 , o t h e r w i s e ,
and
ϵ n , i ( x ) = ( 1 ) i ( 2 i + 1 ) n ( x a ) n 1 .
Proof. 
Without loss of generality, we will show the validity of (11) for [ a , b ] [ 0 , 1 ] , that is, we will prove the validity of the following identity:
D ϕ n , i ( x ) = 2 ( 1 2 n ) p = 0 p o d d i 1 ( p + 1 ) ϕ n , i p 1 ( x ) + 2 p = 0 p e v e n i 1 i ( 2 n + 1 ) + ( 1 2 n ) ( i p 1 ) + 1 ϕ n , i p 1 ( x ) + ( 1 ) i ( 2 i + 1 ) n x n 1 .
Now, to prove (14), we start with the analytic form of the shifted third-kind Chebyshev polynomials given by ([42])
ϕ 0 , i ( x ) = ( 2 i + 1 ) r = 0 i ( 1 ) i + r 4 r ( i + r ) ! ( i r ) ! ( 2 r + 1 ) ! x r ,
to obtain the following formula:
D ϕ n , i ( x ) = ( 2 i + 1 ) x n r = 1 i ( 1 ) i + r 4 r ( i + r ) ! ( n + r ) ( i r ) ! ( 2 r + 1 ) ! x r 1 .
Equation (16) can be written alternatively in the form:
D ϕ n , i ( x ) = ( 1 ) i ( 2 i + 1 ) n x n 1 + ( 2 i + 1 ) r = 1 i ( 1 ) i + r 4 r ( i + r ) ! ( n + r ) ( i r ) ! ( 2 r + 1 ) ! x n + r 1 .
The inversion formula of the shifted polynomials V * ( x ) = ϕ 0 , ( x ) is given by [42]
x r = ( 2 r + 1 ) ! 4 r = 0 r 1 ( r ) ! ( + r + 1 ) ! ϕ 0 , ( x ) .
Inserting Formula (18) into (17) yields the following formula:
D ϕ n , i ( x ) = ( 1 ) i ( 2 i + 1 ) n x n 1 + 2 ( 2 i + 1 ) r = 1 i ( 1 ) i + r ( n + r ) ( i + r ) ! r ( 2 r + 1 ) ( i r ) ! × = 0 r 1 1 ( r 1 ) ! ( r + ) ! ϕ n , ( x ) .
Using some lengthy algebra computations, Formula (19) converts into the following equivalent formula:
D ϕ n , i ( x ) = 4 ( 2 i + 1 ) p = 0 i 1 ( 1 ) p ( i + n p ) ( 2 i p ) ! p ! ( 2 i 2 p + 1 ) ! × 5 F 4 p , i p , i p + 1 2 , 2 i p + 1 , i + n p + 1 2 i 2 p , i p + 1 , i p + 3 2 , i + n p 1 ϕ n , i p 1 ( x ) + ( 1 ) i ( 2 i + 1 ) n x n 1 .
Based on Lemma 1, the hypergeometric function that appears in (20) can be reduced as in Equation (8). This enables one to write (20) in the form
D ϕ n , i ( x ) = p = 0 i 1 G p , i ϕ n , i p 1 ( x ) + ( 1 ) i ( 2 i + 1 ) n x n 1 ,
where G p , i is given by
G p , i = 2 ( 1 2 n ) ( p + 1 ) , p e v e n , 2 ( 2 i + 2 n p + 2 n p ) , p o d d .
Now, if x in (21) is replaced by x a b a , then Formula (11) can be obtained.
Theorem 1 is now proved. □
Remark 2.
It is to be noted that the modified polynomials ϕ n , i ( x ) that are defined in (6) generalize the celebrated third kind Chebyshev polynomials. The following corollary presents the corresponding formula to the special case n = 0 .
Corollary 1.
The first-order derivative of the third-kind Chebyshev polynomials V i ( x ) can be explicitly represented in the form
D V i * ( x ) = 2 p = 0 p o d d i 1 ( p + 1 ) V i p 1 * ( x ) + p = 0 p e v e n i 1 ( 2 i p ) V i p 1 * ( x ) .
Proof. 
Setting n = 0 in Formula (11) yields the following formula
D V i * ( x ) = j = 0 i 1 θ i , j ( 0 ) V i * ( x ) ,
which can be written alternatively as in (23). □
Remark 3.
Formula (23) is a special case of the formula that is obtained in Ref. [27].
Now, assume that the function u ( x ) L w n 2 [ a , b ] can be expanded as
u ( x ) = i = 0 a i ϕ n , i ( x ) ,
where
a i = 2 π ( b a ) a b w n ( x ) u ( x ) ϕ n , i ( x ) d x .
Let the series in Equation (25) have the following approximation:
u ( x ) u N ( x ) = i = 0 N a i ϕ n , i ( x ) = A T Φ n ( x ) ,
where
A T = a 0 , a 1 , , a N , Φ n ( x ) = [ ϕ n , 0 ( x ) , ϕ n , 1 ( x ) , , ϕ n , N ( x ) ] T , Φ 0 ( x ) = [ V 0 * ( x ) , V 1 * ( x ) , , V N * ( x ) ] T .
Now, the following corollary gives the operational matrix of derivatives of Φ n ( x ) . It is a direct consequence of Theorem 1.
Corollary 2.
It is possible to calculate the first derivative of the given vector Φ n ( x ) using the following formula:
d Φ n ( x ) d x = G n Φ n ( x ) + ϵ n ( x ) , n 0 , G 0 Φ 0 ( x ) , n = 0 ,
where ϵ n ( x ) = ϵ n , 0 ( x ) , ϵ n , 1 ( x ) , , ϵ n , N ( x ) T a n d G n = θ i j ( n ) 0 i , j N .
For instance, if N = 6 and n = 4 , we obtain
G n = 4 b a 0 0 0 0 0 0 0 5 0 0 0 0 0 0 7 6 0 0 0 0 0 14 7 7 0 0 0 0 14 15 7 8 0 0 0 23 14 16 7 9 0 0 21 24 14 17 7 10 0 7 × 7 .
Based on Corollary 2, the high-order derivatives of the vector Φ n can be calculated. This outcome is shown by the corollary that follows.
Corollary 3.
For a positive integer m, the following formula can be used to calculate the mth-derivative of the vector Φ n ( x ) :
d m Φ n ( x ) d x m = G n m Φ n ( x ) + η n ( m ) ( x ) , n 0 , G 0 m Φ 0 ( x ) , n = 0 ,
with η n ( m ) ( x ) = k = 0 m 1 G n k ϵ n ( m k 1 ) ( x ) .
Proof. 
Repeated application of (29) yields (31). □

4. A Galerkin Operational Matrix Method for High-Order Singular Type Equations with Initial Conditions

This section provides a numerical solution to the higher-order Emden–Fowler equations using the Galerkin operational matrix of derivatives of the basis functions defined in (6). Consider the Emden–Fowler equation shown below, which has the form ([14]):
u ( p + q ) ( x ) + r = 1 p p r 1 x r j = 1 r ( k j + 1 ) u ( p + q r ) ( x ) + f ( x ) g ( u ) = 0 , 0 x 1 ,
subject to the non-homogeneous initial conditions:
u ( i ) ( 0 ) = α i , i = 0 , 1 , , q 1 , and u ( i ) ( 0 ) = 0 , i = q , q + 1 , , p + q 1 .
An important step in deriving the proposed algorithm is transforming Equation (32) with respect to non-homogeneous initial conditions (33) into the corresponding homogeneous initial conditions equation. To do so, we use the following transformation:
u ˜ ( x ) = u ( x ) i = 0 q 1 α i i ! x i ,
that transforms Equation (32), subject to the non-homogeneous initial conditions (33), into the following modified high-order singular equation:
u ˜ ( p + q ) ( x ) + r = 1 p p r 1 x r j = 1 r ( k j + 1 ) u ˜ ( p + q r ) ( x ) + f ( x ) g ( u ˜ ) = 0 ,
governed by the homogeneous initial conditions
u ˜ ( i ) ( 0 ) = 0 , i = 0 , 1 , , p + q 1 .
Take note of the fact that the single point at x = 0 appears p times at most as x , x 2 , , x p with corresponding distinct shape factors. In addition, (32) reduces to the generalized Lane–Emden-type equations of higher orders for f ( x ) = 1 .
The general equation in (35) involves many specific equations that can be studied separately. We refer here to the following special cases:
  • The second-order Emden–Fowler type equations that can be obtained from Equation (35) by selecting p + q = 2 , that is, we have the following only option for the two parameters p and q: p = 1 , q = 1 .
  • The third-order Emden–Fowler type equations that can be obtained from Equation (35) by selecting p + q = 3 , that is, we have the following two options for the two parameters p and q: p = 2 , q = 1 and p = 1 , q = 2 .
  • The fourth-order Emden–Fowler type equations that can be obtained from Equation (35) by selecting p + q = 4 , p , q 1 , that is, we have the following three options for the two parameters p and q: p = 3 , q = 1 , p = 2 , q = 2 and p = 1 , q = 3 .
Remark 4.
For Emden–Fowler-type equations of nth-order in Equation (32), we should select p + q = n , which leads to the ( n 1 ) choices for the two parameters p and q. They are p = n 1 , q = 1 , p = n 2 , q = 2 , , p = 1 , q = n 1 .
To obtain approximations of the solutions to Equation (35)–(36), we suggest two spectral approaches in the current section: the shifted Chebyshev third-kind Galerkin operational matrix method (SC3GOMM) and the shifted Chebyshev third-kind collocation operational matrix method (SC3COMM). To be specific, if Equation (35) is linear, SC3GOMM is used, and SC3COMM is utilized if it is non-linear.

4.1. SC3GOMM for Handling High-Order Emden–Flower-type Equations

Here, we offer a numerical algorithm for solving Equations (35) and (36). We can describe the approximate solution u ¯ N ( x ) and its high-order derivatives in matrix forms thanks to Formulas (29) and (31). This allows one to write the residual of Equation (35) as
R n , N ( x ) = A T G n p + q Φ n ( x ) + η n ( p + q ) ( x ) + r = 1 p p r 1 x r j = 1 r ( k j + 1 ) A T G n p + q r Φ n ( x ) + η n ( p + q r ) ( x ) + f ( x ) g ( A T Φ n ( x ) ) = 0 ,
with n = p + q .
and as a result, the relation produced by the Galerkin approach is as follows:
( R n , N ( x ) , ϕ n , i ( x ) ) w n 2 = 0 1 w 0 ( x ) R n , N ( x ) ϕ 0 , i ( x ) d x = 0 , i = 0 , 1 , , N , n = p + q .
Then, Equation (38) yields a system of ( N + 1 ) linear/nonlinear equations whose solution can be found by an appropriate solver.

4.2. SC3COMM for Handling High-Order Emden–Flower-type Equations

The shifted Chebyshev third-kind collocation operational matrix approach that we provide in this part is based on using the collocation method to numerically solve Equations (35) and (36). The ( N + 1 ) zeros of V N + 1 * ( x ) are our choice for the collocation points. We denote them by x i , i = 0 , 1 , , N , so we have
R n , N ( x i ) = 0 , i = 0 , 1 , . . . , N , n = p + q .
The approximate solution to the system of equations resulting from the application of the relation (39) can be obtained by solving a system of ( N + 1 ) linear or nonlinear algebraic equations for the unknown coefficients a i using any suitable solver.

5. Tau and Collocation Operational Matrix Methods for Treating High-Order Emden–Flower-Type with Initial Conditions

For solving (32)–(33), we implement two numerical methods in this part. To achieve this goal, the operational matrix of derivatives of the shifted Chebyshev polynomials of the third kind is utilized. We will propose two approaches, shifted Chebyshev third-kind tau method (SC3TM) and shifted Chebyshev third-kind collocation method (SC3CM).

5.1. SC3TM for Handling High-Order Singular Type Equations

In this section, the numerical solution to Equation (32) with the initial conditions (33) is implemented and examined in this section. Our numerical method is based on applying the spectral tau spectral method. The two Formulas (29) and (31) (for the case n = 0 ) are used to write the residual in the form (37).
By using the tau technique, it is possible to derive the following linear/nonlinear equations of dimension N ( p + q ) + 1 in the unknown vector A:
( x p R 0 , N ( x ) , V i * ( x ) ) w 0 = 0 1 x p w 0 ( x ) R 0 , N ( x ) V i * ( x ) d x = 0 , i = 0 , 1 , . . . , N ( p + q ) ,
as well as, the initial conditions (33) lead to the following ( p + q ) equations:
G 0 i Φ 0 ( 0 ) = α i , i = 0 , 1 , , q 1 , and G 0 i Φ 0 ( 0 ) = 0 , i = q , q + 1 , , p + q 1 .
The system of ( N + 1 ) linear/nonlinear equations that make up the set of equations in (40) and (41) can be solved by any appropriate numerical solver to obtain the approximate solution u N ( x ) .

5.2. SC3CM for Handling High-Order Singular Type Equations

In this part, we demonstrate how to use the collocation method to numerically solve Equations (32) with (33). We enforce the residual R 0 , N ( x ) , which is defined by Equation (37), to vanish at the collocation points, which we choose to be the ( N ( p + q ) + 1 ) zeros of the shifted polynomials V N ( p + q ) + 1 * ( x ) on ( 0 , 1 ) . In other words, we have
R 0 , N ( x i ) = 0 , i = 0 , 1 , , N ( p + q ) .
Equations (41) and (42) constitute ( N + 1 ) linear or nonlinear algebraic equations that can be solved by suitable numerical solvers to obtain u N ( x ) .
Remark 5.
The SC3GOMM and SC3TM methods are better suited if Equation (32) is linear because of the anticipated difficulties in the integrals in (38) and (40), respectively, while the SC3COMM and SC3CM methods can be efficiently applied if Equation (32) is linear or nonlinear.
Remark 6.
The choice of basis in Section 4 guarantees the satisfaction of the initial conditions (36), while the choice of basis in Section 5 implies that the initial conditions are set as constraints. This distinguishes the difference between the Galerkin method and other spectral methods.

6. Convergence and Error Analysis

The examination of error analysis for the two proposed expansions is presented in this part. In this context, the following two theorems are presented and demonstrated.
Theorem 2.
Assume that f ( x ) L w q [ 0 , 1 ] with | f ( q ) ( x ) | M , M is a positive constant. Assume that f ( x ) has the following expansion:
f ( x ) = i = 0 c i ϕ 0 , i ( x ) .
The coefficients c i fulfill the following inequality
| c i | < M 2 q 1 1 ( i q + 1 ) q , i > q 1 ,
and the series converges uniformly to f ( x ) .
Proof. 
The orthogonality relation (5) clearly indicates that
c i = 2 π 0 1 x 1 x f ( x ) ϕ 0 , i d x .
In (44), setting 2 x 1 = cos θ yields
c i = 2 π 0 π f 1 + cos θ 2 cos ( θ / 2 ) cos ( ( i + 1 / 2 ) ) θ d θ .
If the right-hand side of (45) is integrated by parts q times, then we obtain
c i = 1 2 q 1 π 0 π f ( q ) 1 + cos θ 2 R i ( θ ) d θ ,
where R i ( θ ) is a trigonometric polynomial in cos θ and sin θ . So, using the two well-known simple inequalities: | cos θ | 1 , | sin θ | 1 and with a little computation, we can achieve our goal. □
Theorem 3.
Let y N ( x ) = i = 0 N c i ϕ 0 , i ( x ) be an approximate solution to (32). The following inequality is satisfied by the global error under the same hypotheses of Theorem 2:
y y N w < M 2 q 3 2 1 N q 1 2 .
Proof. 
The orthogonality property of V i * ( x ) , i 0 immediately yields
y y N w 0 2 = π 2 i = N + 1 c i 2 .
If we apply Theorem 2, then Identity (47) leads to the following inequality:
y y N w 0 2 < M 2 π 2 2 q 1 i = N + 1 1 ( i q + 1 ) 2 q ,
and as a result, the following inequality also holds: (see, [43])
y y N w 0 2 < M 2 π 2 2 q 1 N 1 ( x q + 1 ) 2 q d x .
Finally, the inequality (48) implies that
y y N w 0 < M 2 q 3 2 1 N q 1 2 .
This proves Theorem 3. □
Remark 7.
The results of two Theorems 2 and 3 also hold if f ( x ) is expanded as f ( x ) = i = 0 c ¯ n , i ϕ n , i ( x ) , where the coefficients c ¯ n , i are defined by the relation:
c ¯ n , i = 2 π 0 1 1 x 2 n x 1 x f ( x ) ϕ n , i d x .

7. Numerical Results

This section focuses on providing some examples to demonstrate how easily our suggested techniques can be used. Additionally, some comparisons with several other approaches are provided.
Example 1.
Consider the non-linear Emden–Flower-type Equation (32), for the special case of ( p = 2 , q = 1 , k = 3 ) [11,44,45]
u ( x ) + 6 x u ( x ) + 6 x 2 u ( x ) = 6 ( x 6 + 2 x 3 + 10 ) e 3 u ( x ) , 0 x 1 ,
subject to the initial conditions
u ( 0 ) = 0 , u ( 0 ) = 0 , u ( 0 ) = 0 .
The exact solution of (50) is u ( x ) = ln ( 1 + x 3 ) .
Table 1 shows the maximum errors that our method SC3COMM might produce, whereas Table 2 compares our technique to the following methods:
  • B-spline collocation method (BSCM) developed in [11].
  • Varitional itereration method (VIM) in [44].
  • Quartic B-spline method (QBSM) in [45].
Additionally, Figure 1 displays a comparison between the obtained error using the two presented methods, SC3COMM and SC3CM.
Table 1. Maximum absolute error of Example 1.
Table 1. Maximum absolute error of Example 1.
NSC3COMMSC3CM
6 2.04719 × 10 5 1.19833 × 10 3
9 5.30930 × 10 8 2.04719 × 10 5
15 7.59007 × 10 12 2.22160 × 10 9
21 5.52367 × 10 15 4.62410 × 10 13
22 4.54192 × 10 15 9.12330 × 10 14
Table 2. Comparison between different methods of Example 1.
Table 2. Comparison between different methods of Example 1.
SC3COMM
( N = 22 )
SC3CM
( N = 22 )
VIM [44]QBSM [45]BSCM [11]
4.542 × 10 15 2.123 × 10 14 1.10 × 10 1 1.29 × 10 5 4.08 × 10 7
Example 2.
Consider the nonlinear Emden–Flower-type Equation (32), for the special case of ( p = 2 , q = 2 , k = 4 ) [14,15]
u ( 4 ) ( x ) + 8 x u ( 3 ) ( x ) + 12 x 2 u ( 2 ) ( x ) + u m ( x ) = 0 , 0 x 1 ,
subject to the initial conditions
u ( 0 ) = 1 , u ( i ) ( 0 ) = 0 , i = 1 , 2 , 3 .
The exact solution of (51) is u ( x ) = 1 x 4 360 , for m = 0 .
The use of our two suggested techniques S C 3 G O M M and S C 3 C O M M for N = 0 , 1 , 2 , , gives the follwoing approximate solution: u N ˜ ( x ) = i = 0 N a i ϕ i ( x ) , where a 0 = 1 360 and a i = 0 , i = 1 , 2 , , N . In view of (34), we obtain u N ( x ) = 1 x 4 360 , which is the exact solution. Additionally, the application of two proposed methods S C 3 T M and S C 3 C M for N = 4 , 5 , 6 , , gives u N ( x ) = i = 0 N a i V i * ( x ) , where a 0 = 5113 5120 , a 1 = 7 7680 , a 2 = 1 2560 , a 3 = 1 10240 ,   a 4 = 1 92160 , and a i = 0 , i = 5 , 6 , 7 , . . . , N .
Example 3.
Consider the nonlinear Emden–Flower-type Equation (32), for the special case of ( p = 3 , q = 1 , k = 4 ) [14,15]
u ( 4 ) ( x ) + 12 x u ( 3 ) ( x ) + 36 x 2 u ( 2 ) ( x ) + 24 x 3 u ( 1 ) ( x ) + 60 ( 7 18 x 4 + 3 x 8 ) u 9 ( x ) = 0 , 0 x 1 ,
subject to the boundary conditions
u ( 0 ) = 1 , u ( i ) ( 0 ) = 0 , i = 1 , 2 , 3 .
The exact solution of (52) is u ( x ) = 1 x 4 + 1 .
The maximum errors introduced by applying our suggested method SC3COMM for various values of specific parameters N are shown in Table 3, while Table 4 compares the results obtained from our algorithm with the following two methods:
  • The Adomian decomposition method (ADM) that is applied in [14].
  • The reproducing kernel Hilbert space method (RKHSM) that is derived in [15].
Additionally, Figure 2 displays the Log-error for different N to demonstrate the stability of the solutions to Example 3 if the two methods SC3COMM and SC3CM are employed, respectively.
Table 3. Maximum absolute error of Example 3.
Table 3. Maximum absolute error of Example 3.
NSC3COMMSC3CM
6 1.046 × 10 4 2.284 × 10 2
9 1.459 × 10 6 4.187 × 10 4
15 2.740 × 10 10 8.886 × 10 8
21 5.664 × 10 12 1.606 × 10 11
26 6.664 × 10 14 1.012 × 10 11
Table 4. Comparison between different methods of Example 3.
Table 4. Comparison between different methods of Example 3.
SC3COMM
( N = 26 )
SC3CM
( N = 21 )
RKHSM [15]ADM [14]
6.664 × 10 14 1.606 × 10 11 1.23016 × 10 6 5.21075 × 10 5
Example 4.
Consider the following nonlinear Emden–Flower-type Equation (32), for the special case of ( p = 3 , q = 1 , k = 1 ) [14,15,46]
u ( 4 ) ( x ) + 3 x u ( 3 ) ( x ) 96 ( 1 10 x 4 + 5 x 8 ) e 4 u = 0 , 0 x 1 ,
subject to the boundary conditions
u ( i ) ( 0 ) = 0 , i = 0 , 1 , 2 , 3 .
The exact solution of (53) is u ( x ) = ln ( 1 + x 4 ) .
The maximum errors produced by applying our two approaches, S C 3 C O M M and S C 3 C M , are shown in Table 5, and Table 6 compares these two methods to the following methods:
  • The Adomian decomposition method (ADM) that applied in [14].
  • The reproducing kernel Hilbert space method (RKHSM) that derived in [15].
  • The Quartic B-spline method (QBSM) that developed in [46].
Table 5. Maximum absolute errors of Example 4.
Table 5. Maximum absolute errors of Example 4.
NSC3COMMSC3CM
6 9.323 × 10 5 9.030 × 10 2
13 6.421 × 10 9 4.312 × 10 6
17 2.370 × 10 10 6.421 × 10 9
22 1.512 × 10 5 7.179 × 10 13
Table 6. Comparison between different methods of Example 4.
Table 6. Comparison between different methods of Example 4.
SC3COMM
( N = 22 )
SC3CM
( N = 22 )
RKHSM [15]ADM [14]QBSM [46]
1.512 × 10 13 7.179 × 10 11 1.98579 × 10 6 7.56423 × 10 5 6.88 × 10 7
Example 5.
Consider the linear Emden–Flower-type Equation (32), for the special case of ( p = 4 , q = 1 , k = 4 )
u ( 5 ) + 16 x u ( 4 ) + 72 x 2 u ( 3 ) + 96 x 3 u ( 2 ) + 24 x 4 u ( 1 ) u = f ( x ) , 0 x 1 ,
subject to the initial conditions
u ( i ) ( 0 ) = 0 , i = 0 , 1 , 2 , 3 , 4 ,
where f ( x ) = ( 8400 + 9744 x + 3696 x 2 + 592 x 3 + 41 x 4 ) e x . The exact solution of (54) is u ( x ) = x 5 e x .
Table 7 displays the maximum errors that resulted from the application of our four methods SC3COMM, SC3GOMM, SC3CM and SC3TM. Again, Figure 3 presents Log-error for various N to show the stability of solutions of Example 5 if the presented four methods are applied.

8. Results and Discussions

New spectral solutions to singular high-order type equations were presented in this paper. New types of modified Chebyshev polynomials of the third kind and their well-known classical types were used as basis functions. The core of the Galerkin matrix method that was designed in this paper for treating the high-order Emden–Fowler equations is based on the derivation of the first-order derivative of the modified Chebyshev polynomials in Theorem 1 and novel operational matrices of the derivatives of these modified polynomials in Corollary 2. The convergence analysis of the two proposed expansions is carefully investigated in Theorems 2 and 3. Numerical results in Section 7 show the efficiency and applicability of our algorithm. We think that the modified third-kind Chebyshev polynomials proposed in this paper could be used to handle diverse differential equations.

Author Contributions

W.M.A.-E. contributed to conceptualization, methodology, software, validation, formal analysis, investigation, Writing—original draft, Writing—review & editing, and Supervision. M.S.A.-H. contributed to validation and investigation, A.K.A. contributed to validation and funding Acquisition. H.M.A. contributed to conceptualization, methodology, software, formal analysis, investigation, Writing—original draft, Writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

The third author: Amr Kamel Amin ([email protected]) is funded by the Deanship of Scientific Research at Umm Al-Qura University.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work by Grant Code: (22UQU4331287DSR07).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Lane, H.J. On the theoretical temperature of the sun, under the hypothesis of a gaseous mass maintaining its volume by its internal heat, and depending on the laws of gases as known to terrestrial experiment. Am. J. Sci. Arts 1870, 2, 57–74. [Google Scholar] [CrossRef] [Green Version]
  2. Emden, R. Gaskugeln: Anwendungen der mechanischen Wärmetheorie auf kosmologische und meteorologische Probleme; Teubner: Leipzig, Germany, 1907. [Google Scholar]
  3. Fowler, R.H. Further studies of Emden’s and similar differential equations. Q. J. Math. 1931, 1, 259–288. [Google Scholar] [CrossRef]
  4. Singh, R.; Garg, H.; Guleria, V. Haar wavelet collocation method for Lane–Emden equations with Dirichlet,Neumann and Neumann–Robin boundary conditions. J. Compt. Appl. Math. 2019, 346, 150–161. [Google Scholar] [CrossRef]
  5. Dizicheh, A.K.; Salahshour, S.; Ahmadian, A.; Baleanu, D. A novel algorithm based on the Legendre wavelets spectral technique for solving theLane–Emden equations. Appl. Numer. Math. 2020, 153, 443–456. [Google Scholar] [CrossRef]
  6. Abd-Elhameed, W.M.; Ahmed, H.M. Tau and Galerkin operational matrices of derivatives for treating singular and Emden–Fowler third-order-type equations. Int. J. Mod. Phys. 2022, 33, 2250061. [Google Scholar] [CrossRef]
  7. Abd-Elhameed, W.M. New Galerkin operational matrix of derivatives for solving Lane-Emden singular-type equations. Eur. Phys. J. Plus 2015, 130, 1–12. [Google Scholar] [CrossRef]
  8. Wazwaz, A.M.; Rach, R.; Bougoffa, L.; Duan, J.S. Solving the Lane–Emden–Fowler type equations of higher orders by the Adomian decomposition method. CMES Comput. Model. Eng. Sci. 2014, 100, 507–529. [Google Scholar]
  9. Hasan, Y.Q.; Zhu, L.M. Solving singular boundary value problems of higher-order ordinary differential equations by modified Adomian decomposition method. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 2592–2596. [Google Scholar] [CrossRef]
  10. Kim, W.; Chun, C. A modified Adomian decomposition method for solving higher-order singular boundary value problems. Z. Für Naturforschung A 2010, 65, 1093–1100. [Google Scholar] [CrossRef]
  11. Iqbal, M.K.; Abbas, M.; Wasim, I. New cubic B-spline approximation for solving third order Emden–Flower type equations. Appl. Math. Comput. 2018, 331, 319–333. [Google Scholar] [CrossRef]
  12. Roul, P.; Thula, K. A fourth-order B-spline collocation method and its error analysis for Bratu-type and Lane–Emden problems. Int. J. Comput. Math. 2019, 96, 85–104. [Google Scholar] [CrossRef]
  13. Wazwaz, A.M. The variational iteration method for solving systems of third-order Emden-Fowler type equations. J. Math. Chem. 2017, 55, 799–817. [Google Scholar] [CrossRef]
  14. Wazwaz, A.M.; Rach, R.; Duan, J.S. Solving new fourth–order Emden–Fowler-type equations by the Adomian decomposition method. Int. J. Comput. Methods Eng. Sci. 2015, 16, 121–131. [Google Scholar] [CrossRef]
  15. Dezhbord, A.; Lotfi, T.; Mahdiani, K. A numerical approach for solving the high-order nonlinear singular Emden–Fowler type equations. Adv. Differ. Equ. 2018, 2018, 161. [Google Scholar] [CrossRef]
  16. Madduri, H.; Roul, P. A fast-converging iterative scheme for solving a system of Lane–Emden equations arising in catalytic diffusion reactions. J. Math. Chem. 2019, 57, 570–582. [Google Scholar] [CrossRef]
  17. Canuto, C.; Hussaini, M.Y.; Quarteroni, A.; Zang, T.A. Spectral Methods in Fluid Dynamics; Springer: Berlin/Heidelberg, Germany, 1988. [Google Scholar]
  18. Hesthaven, J.; Gottlieb, S.; Gottlieb, D. Spectral Methods for Time-Dependent Problems; Cambridge University Press: Cambridge, UK, 2007; Volume 21. [Google Scholar]
  19. Boyd, J.P. Chebyshev and Fourier Spectral Methods; Courier Corporation: Chelmsford, MA, USA, 2001. [Google Scholar]
  20. Trefethen, L.N. Spectral Methods in MATLAB; SIAM: Philadelphia, PA, USA, 2000; Volume 10. [Google Scholar]
  21. Ahmed, H.M. Numerical solutions of Korteweg-de Vries and Korteweg-de Vries-Burger’s equations in a Bernstein polynomial basis. Mediterr. J. Math. 2019, 16, 102. [Google Scholar] [CrossRef]
  22. Abd-Elhameed, W.M.; Alkenedri, A.M. Spectral solutions of linear and nonlinear BVPs using certain Jacobi polynomials generalizing third-and fourth-kinds of Chebyshev polynomials. CMES Comput. Model. Eng. Sci. 2021, 126, 955–989. [Google Scholar] [CrossRef]
  23. Napoli, A.; Abd-Elhameed, W.M. An innovative harmonic numbers operational matrix method for solving initial value problems. Calcolo 2017, 54, 57–76. [Google Scholar] [CrossRef]
  24. Doha, E.H.; Abd-Elhameed, W.M.; Bhrawy, A.H. New spectral-Galerkin algorithms for direct solution of high even-order differential equations using symmetric generalized Jacobi polynomials. Collect. Math. 2013, 64, 373–394. [Google Scholar] [CrossRef]
  25. Abdelhakem, M.; Alaa-Eldeen, T.; Baleanu, D.; Alshehri, M.G.; El-Kady, M. Approximating real-life BVPs via Chebyshev polynomials’ first derivative pseudo-Galerkin method. Fractal Frac. 2021, 5, 165. [Google Scholar] [CrossRef]
  26. Doha, E.H.; Abd-Elhameed, W.M.; Bassuony, M.A. On using third and fourth kinds Chebyshev operational matrices for solving Lane-Emden type equations. Rom. J. Phys 2015, 60, 281–292. [Google Scholar]
  27. Doha, E.H.; Abd-Elhameed, W.M.; Bassuony, M.A. On the coefficients of differentiated expansions and derivatives of Chebyshev polynomials of the third and fourth kinds. Acta Math. Sci. 2015, 35, 326–338. [Google Scholar] [CrossRef]
  28. Alsuyuti, M.M.; Doha, E.H.; Ezz-Eldien, S.S.; Bayoumi, B.I.; Baleanu, D. Modified Galerkin algorithm for solving multitype fractional differential equations. Math. Methods Appl. Sci. 2019, 42, 1389–1412. [Google Scholar] [CrossRef]
  29. Abd-Elhameed, W.M. Novel expressions for the derivatives of sixth-kind Chebyshev polynomials: Spectral solution of the non-linear one-dimensional Burgers’ equation. Fractal Fract. 2021, 5, 53. [Google Scholar] [CrossRef]
  30. Tu, H.; Wang, Y.; Lan, Q.; Liu, W.; Xiao, W.; Ma, S. A Chebyshev-Tau spectral method for normal modes of underwater sound propagation with a layered marine environment. J. Sound Vib. 2021, 492, 115784. [Google Scholar] [CrossRef]
  31. Mokhtary, P.; Ghoreishi, F.; Srivastava, H.M. The Müntz-Legendre Tau method for fractional differential equations. Appl. Math. Model. 2016, 40, 671–684. [Google Scholar] [CrossRef]
  32. Mohebbi, A. Crank–Nicolson and Legendre spectral collocation methods for a partial integro-differential equation with a singular kernel. J. Comput. Appl. Math. 2019, 349, 197–206. [Google Scholar] [CrossRef]
  33. Çelik, I. Gegenbauer wavelet collocation method for the extended Fisher-Kolmogorov equation in two dimensions. Math. Methods Appl. Sci. 2020, 43, 5615–5628. [Google Scholar] [CrossRef]
  34. Khader, M.M.; Adel, M. Numerical and theoretical treatment based on the compact finite difference and spectral collocation algorithms of the space fractional-order Fisher’s equation. Int. J. Mod. Phys. C 2020, 31, 2050122. [Google Scholar] [CrossRef]
  35. Wang, F.; Zhao, Q.; Chen, Z.; Fan, C.M. Localized Chebyshev collocation method for solving elliptic partial differential equations in arbitrary 2D domains. Appl. Math. Comput. 2021, 397, 125903. [Google Scholar] [CrossRef]
  36. Bansu, H.; Kumar, S. Numerical solution of space-time fractional Klein-Gordon equation by radial basis functions and Chebyshev polynomials. Inter. J. Appl. Comput. Math. 2021, 7, 1–19. [Google Scholar] [CrossRef]
  37. Nouri, K.; Fahimi, M.; Torkzadeh, L.; Baleanu, D. Numerical method for pricing discretely monitored double barrier option by orthogonal projection method. AIMS Math. 2021, 6, 5750–5761. [Google Scholar] [CrossRef]
  38. Sakran, M.R.A. Numerical solutions of integral and integro-differential equations using Chebyshev polynomials of the third kind. Appl. Math. Comp. 2019, 351, 66–82. [Google Scholar] [CrossRef]
  39. Balaji, S.; Hariharan, G. An efficient operational matrix method for the numerical solutions of the fractional Bagley–Torvik equation using wavelets. J. Math. Chem. 2019, 57, 1885–1901. [Google Scholar] [CrossRef]
  40. Loh, J.R.; Phang, C. Numerical solution of Fredholm fractional integro-differential equation with Right-sided Caputo’s derivative using Bernoulli polynomials operational matrix of fractional derivative. Mediterr. J. Math. 2019, 16, 28. [Google Scholar] [CrossRef]
  41. Koepf, W. Hypergeometric Summation, 2nd ed.; Universitext Series; Springer: London, UK, 2014. [Google Scholar]
  42. Andrews, G.E.; Askey, R.; Roy, R. Special Functions; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  43. Stewart, J. Single Variable Essential Calculus: Early Transcendentals; Cengage Learning: Davis Drive Belmont, CA, USA, 2012. [Google Scholar]
  44. Wazwaz, A.M. Solving two Emden-Fowler type equations of third order by the variational iteration method. Appl. Math. Inf. Sci. 2015, 9, 2429. [Google Scholar]
  45. Mishra, H.K.; Saini, S. Quartic B-spline method for solving a singular singularly perturbed third-order boundary value problems. Am. J. Numer. Anal. 2015, 3, 18–24. [Google Scholar]
  46. Iqbal, M.; Abbas, M.; Zafar, B. New Quartic B-spline approximations for numerical solution of fourth order singular boundary value problems. J. Math. 2020, 52, 47–63. [Google Scholar]
Figure 1. Comparison between the errors of Example 1 for N = 22 .
Figure 1. Comparison between the errors of Example 1 for N = 22 .
Axioms 12 00099 g001
Figure 2. Comparison between the errors of Example 3 for N = 3 , 4 , , 26 .
Figure 2. Comparison between the errors of Example 3 for N = 3 , 4 , , 26 .
Axioms 12 00099 g002
Figure 3. Comparison between the errors of Example 5 for N = 3 , 4 , , 17 .
Figure 3. Comparison between the errors of Example 5 for N = 3 , 4 , , 17 .
Axioms 12 00099 g003
Table 7. Maximum absolute error of Example 5.
Table 7. Maximum absolute error of Example 5.
NSC3COMMSC3GOMMSC3CMSC3TM
6 4.373 × 10 8 4.249 × 10 8 1.863 × 10 1 1.910 × 10 1
9 4.362 × 10 13 3.672 × 10 13 9.633 × 10 5 2.019 × 10 4
11 5.450 × 10 16 4.926 × 10 16 4.373 × 10 8 6.944 × 10 6
16 1.671 × 10 16 1.221 × 10 16 9.710 × 10 15 2.163 × 10 11
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abd-Elhameed, W.M.; Al-Harbi, M.S.; Amin, A.K.; M. Ahmed, H. Spectral Treatment of High-Order Emden–Fowler Equations Based on Modified Chebyshev Polynomials. Axioms 2023, 12, 99. https://doi.org/10.3390/axioms12020099

AMA Style

Abd-Elhameed WM, Al-Harbi MS, Amin AK, M. Ahmed H. Spectral Treatment of High-Order Emden–Fowler Equations Based on Modified Chebyshev Polynomials. Axioms. 2023; 12(2):99. https://doi.org/10.3390/axioms12020099

Chicago/Turabian Style

Abd-Elhameed, Waleed Mohamed, Mohamed Salem Al-Harbi, Amr Kamel Amin, and Hany M. Ahmed. 2023. "Spectral Treatment of High-Order Emden–Fowler Equations Based on Modified Chebyshev Polynomials" Axioms 12, no. 2: 99. https://doi.org/10.3390/axioms12020099

APA Style

Abd-Elhameed, W. M., Al-Harbi, M. S., Amin, A. K., & M. Ahmed, H. (2023). Spectral Treatment of High-Order Emden–Fowler Equations Based on Modified Chebyshev Polynomials. Axioms, 12(2), 99. https://doi.org/10.3390/axioms12020099

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop