[BACK]
imagesComputer Modeling in Engineering & Sciences images

DOI: 10.32604/cmes.2021.013603

ARTICLE

Spectral Solutions of Linear and Nonlinear BVPs Using Certain Jacobi Polynomials Generalizing Third- and Fourth-Kinds of Chebyshev Polynomials

W. M. Abd-Elhameed1,2,*and Asmaa M. Alkenedri2

1Department of Mathematics, Faculty of Science, Cairo University, Giza, Egypt
2Department of Mathematics, College of Science, University of Jeddah, Jeddah, Saudi Arabia
*Corresponding Author: W. M. Abd-Elhameed. Email: walee_9@yahoo.com
Received: 12 August 2020; Accepted: 26 November 2020

Abstract: This paper is dedicated to implementing and presenting numerical algorithms for solving some linear and nonlinear even-order two-point boundary value problems. For this purpose, we establish new explicit formulas for the high-order derivatives of certain two classes of Jacobi polynomials in terms of their corresponding Jacobi polynomials. These two classes generalize the two celebrated non-symmetric classes of polynomials, namely, Chebyshev polynomials of third- and fourth-kinds. The idea of the derivation of such formulas is essentially based on making use of the power series representations and inversion formulas of these classes of polynomials. The derived formulas serve in converting the even-order linear differential equations with their boundary conditions into linear systems that can be efficiently solved. Furthermore, and based on the first-order derivatives formula of certain Jacobi polynomials, the operational matrix of derivatives is extracted and employed to present another algorithm to treat both linear and nonlinear two-point boundary value problems based on the application of the collocation method. Convergence analysis of the proposed expansions is investigated. Some numerical examples are included to demonstrate the validity and applicability of the proposed algorithms.

Keywords: Jacobi polynomials; high-order boundary value problems; Galerkin method; collocation method; connection problem; convergence analysis

1  Introduction

It is well-known that the Chebyshev polynomials play vital roles in the scope of mathematical analysis and its applications. The first- and second-kinds are special symmetric polynomials of the Jacobi polynomials, so they are ultraspherical polynomials. These two kinds of Chebyshev polynomials are the most popular kinds, and they are employed extensively in numerical analysis, see for example [14]. There are other four kinds of Chebyshev polynomials. Third- and fourth-kinds are also special kinds of Jacobi polynomials, but they differ from the first- and second-kinds since they are not ultraspherical polynomials. In fact, they are non-symmetric Jacobi polynomials with certain parameters. There are considerable publications devoted to investigating these kinds of polynomials from both theoretical and practical points of view. For some articles in this direction, one can be referred to [57]. Fifth- and sixth-kinds of Chebyshev polynomials are not special cases of the Jacobi polynomials. These polynomials were first established by Jamei in his Ph.D. dissertation [8]. Recently, these kinds of polynomials have been employed for treating numerically some types of differential equations, see [911].

Spectral methods are essential in numerical analysis. These methods have been utilized fruitfully for obtaining numerical solutions of various kinds of differential and integral equations. The philosophy in the application of spectral methods is based on writing the desired numerical solution in terms of certain “basis functions,” which may be expressed as combinations of various orthogonal polynomials. There are three celebrated spectral methods; they are tau, collocation, and Galerkin methods. Galerkin method is based on choosing suitable combinations of orthogonal polynomials satisfying the underlying boundary/initial conditions [12]. Tau method is more flexible than the Galerkin method due to the freedom in choosing the basis functions. Collocation method is the most popular method because of its capability to treat all kinds of differential equations. There are considerable contributions employ this method in different types of differential equations [1315]. For a survey on spectral methods and their applications, one can consult [1621].

The problems of establishing derivatives and integrals formulas of different orthogonal polynomials in terms of their original polynomials are of great interest. These formulas play important parts in obtaining spectral solutions of various differential equations. For example, the authors in [5,6] derived explicit formulas of the high-order derivatives and repeated integrals of Chebyshev polynomials of third- and fourth-kinds, and after that, they utilized the derived formulas for obtaining spectral solutions of some boundary value problems (BVPs).

Boundary value problems have important roles due to their appearance in almost all branches related to applied sciences such as engineering, fluid mechanics, and optimization theory (see, [22]). Some examples concerning applications of BVPs are given with details in [22]. Moreover, the author showed that, at some times, it is more appropriate to transform high-order BVPs into systems of differential equations. Many real-life phenomena can be described by high-order BVPs.

Regarding the even-order BVPs, they arise in numerous problems. The free vibrations analysis of beam structures is governed by a fourth-order ordinary differential equation [23]. The vibration behavior of a ring structures is governed by a sixth-order ordinary differential equation [24]. Eighth-order BVPs arise in torsional vibration of uniform beams [25]. In addition, 10th- and 12th-orders appear in applications. When a uniform magnetic field is applied across the fluid in the same direction as gravity. Ordinary convection and overstability yield a 10th-order and a 12th-order BVPs, respectively [26]. For applications regarding the high-order BVPs appear in hydrodynamic and hydromagnetic stability, fluid dynamics, astronomy, beam and long wave theory [27,28]. For some other applications of these kinds of problems, one can be referred to [2931].

Several studies were performed for the numerical solutions of these types of equations. Some of the algorithms treat these kinds are the perturbation and homotopy perturbation methods in [32], differential transform methods in [26] and spectral methods in [15,3337].

The main objectives of this article can be summarized in the following items:

Deriving high-order derivatives formulas of certain classes of polynomials that generalize Chebyshev polynomials of the third- and fourth-kinds.

Implementing and analyzing an algorithm for the numerical solutions of even-order linear BVPs based on the application of the spectral Galerkin method.

Introducing a collocation algorithm for treating both linear and nonlinear BVPs.

The current paper is organized as follows. In Section 2, we give some relevant properties of Jacobi polynomials. In Section 3, we derive in detail two new formulas which give explicitly the high-order derivatives of Jacobi polynomials whose parameters difference is one. In Section 4, we present a spectral Galerkin algorithm for solving some high even-order linear BVPs with constant coefficients. Section 5 is devoted to presenting another algorithm for treating linear and nonlinear BVPs based on the application of a suitable collocation algorithm. Section 6 is interested in solving the connection problem between certain Jacobi polynomials and the Legendre polynomials. In Section 7, we investigate the convergence and error analysis of the proposed expansion in Section 5. In Section 8, we present some numerical examples including comparisons with some other methods to ascertain the accuracy and efficiency of the proposed algorithms presented in Sections 4 and 5. Finally, some conclusions are reported in Section 9.

2  Some Fundamental Properties of Jacobi Polynomials

The classical Jacobi polynomials images, is a polynomial of degree m which can be defined in hypergeometric form as

images

where

images

It is easy to see the following special values

images

It is desirable to define the normalized orthogonal Jacobi polynomials as [38]

images

which gives images. We have the following six special cases of the normalized Jacobi polynomials.

images

where images, Tm(x), Um(x), Vm(x), Wm(x) and Pm(x) are the ultraspherical, Chebyshev of the first-, second-, third- and fourth-kinds, and the Legendre polynomials, respectively.

The polynomials images may be generated using the following recurrence relation:

images

with the initial values:

images

The orthogonality relation satisfied by images is given by:

images

where images is the well-known Kronecker delta function, and

images

The special values

images

images

will be of important use later.

It is useful to extend the Jacobi polynomials on a general interval [a, b]. In fact, the shifted Jacobi polynomials can be defined on [a, b] by

images

Hence, properties and relations concerned with the Jacobi polynomials images can be easily transformed to give their counterparts for the shifted polynomials images.

For our upcoming computations, we need the orthogonality relation of images on [a, b]. This relation is given explicitly as:

images

where images is given by:

images

3  The imagesth Derivative of images and images

The main objective of this section is to state and prove two theorems which express the high-order derivatives of the shifted Jacobi polynomials images and images in terms of their original polynomials. The basic idea behind the derivation of the desired formulas is built on employing the power form representation and the inversion formula of the Jacobi polynomials images. Now, the following three lemmas are essential in this regard.

Lemma 3.1. The explicit power form of the normalized Jacobi polynomials images, images, is given by

images

where

images

and

images

where [z] represents the largest integer less than or equal to z.

Proof. We proceed by induction on n. Assume that relation (10) is valid for (n −1) and (n −2), and we have to prove the validity of (10) itself. Starting with the recurrence relation (3), (for the case images) in the form

images

where

images

then, the application of the induction hypothesis twice yields

images

Eq. (14) may be written in the form

images

where

images

and images

It is not difficult to see that

images

where an, k, bn, k are given by (11) and (12), respectively. Therefore

images

and

images

and this completes the proof of formula (10).

Lemma 3.2. For all images, the following inversion formula holds

images

where

images

and

images

Proof. We proceed by induction on m. Form = 1, the left-hand side equals the right-hand side. Assume that relation (15) is valid, and we have to show that

images

Multiplying both sides of relation (15) by x, then with the aid of relation (3), (for the case images), and after some rather manipulation, Eq. (16) can be proved.

Lemma 3.3. For every nonnegative integer m, we have

images

Proof. If we set

images

then, with the aid of Zeilberger’s algorithm [39], Sn, q, m satisfies the following first-order recurrence relation:

images

with the initial value

images

The exact solution for this recurrence relation is

images

This proves Lemma 3.3.

Now, we are in a position to state and prove the desired derivatives formula.

Theorem 3.1. The derivatives of the normalized shifted Jacobi polynomials images of any degree and for any order in terms of images are given explicitly by

images

where

images

and

images

Proof. We prove the theorem for images, that is, we prove the following formula:

images

where

images

and

images

If we differentiate Eq. (10) q-times with respect to x, and noting the identity

images

then, we get

images

where

images

images

Making use of relation (15), enables one to write

images

where

images

and

images

Expanding images and images and collecting similar terms, then after some rather manipulation, we get

images

where

images

and

images

Now, it is not difficult to show that

images

and

images

Finally, making use of Lemma 3.3, An, m, q and Bn, m, q take the forms

images

and

images

Replacing x in (22) by images, formula (19) can be obtained. This completes the proof of Theorem 3.1.

Remark 3.1. It is to be noted here that relation (19) may be written in the alternative form

images

where

images

Remark 3.2. As a direct consequence of relation (30), the qth derivative of images can be easily deduced. This result is given in the following theorem.

Theorem 3.2. The derivatives of the normalized Jacobi polynomials images of any degree and for any order in terms of images are given explicitly by

images

where images and images are given as:

images

and images and images are given, respectively, in (20) and (21), or alternatively in the form

images

where

images

and images is given by (31).

As two direct consequences of Theorems 3.1 and 3.2, the formulas for the high-order derivatives of Chebyshev polynomials of third- and fourth-kinds in terms of their corresponding Chebyshev polynomials can be deduced. These results are given in the following two corollaries.

Corollary 3.1. The derivatives of Chebyshev polynomials of third-kind Vn(x) on [ −1, 1] of any degree and for any order in terms of their original polynomials are given by the following explicit formula:

images

Corollary 3.2. The derivatives of Chebyshev polynomials of fourth-kind Wn(x) on [ −1, 1] of any degree and for any order in terms of their original formulas are given by

images

Remark 3.3. The results of Corollaries 3.1 and 3.2 are in complete agreement with those obtained in [6].

Remark 3.4. From now on, we will employ only the polynomials images and their shifted ones images which will be denoted, respectively, by images and images. All corresponding results using the polynomials images and their shifted ones images can be obtained similarly.

4  The Solution of High Even-Order Differential Equations Using images

This section is confined to analyzing in detail a spectral Galerkin algorithm for treating numerically even-order BVPs. We will employ the shifted Jacobi polynomials images.

Now, consider the following even-order BVPs:

images

governed by the boundary conditions:

images

where images, and images are constant coefficients and images is defined as

images

The philosophy of applying the Galerkin method to solve (36)(37) is based on choosing basis functions satisfying the boundary conditions (37), and after that, enforcing the residual of Eq. (36) to be orthogonal to these functions. Now, if we set

images

then, to apply the shifted Jacobi–Galerkin method for solving (36)(37), we have to find images such that

images

where images, and images is the scalar inner product in the weighted space L2w(a, b).

4.1 The Choice of Basis Functions

In this section, we describe how to choose suitable basis functions satisfying the boundary conditions (37) to be able to apply the Galerkin method to (36)(37). The idea for choosing compact combinations satisfying the underlying boundary conditions was first introduced in [40,41]. In these two papers, the author selected suitable basis functions satisfying the standard second- and fourth-order BVPs which are considered as special cases of (36). The basis functions were expressed in terms of Chebyshev polynomials of the first-kind and Legendre polynomials. The authors in [42] generalized the proposed algorithms in [40,41]. They suggested compact combinations of ultraspherical polynomials that generalize the Chebyshev of the first-kind and Legendre combinations to solve the even-order BVPs in one- and two-dimensions. The authors in [43] solved the same types of BVPs but by using certain kinds of non-symmetric Jacobi polynomials which are called Chebyshev polynomials of the third- and fourth-kinds. In this section, we are going to generalize the basis functions that derived in [43].

First, consider the case images, and set

images

where the coefficients images are chosen such that every member of images verifies the following homogeneous boundary conditions:

images

It is not difficult to see that the boundary conditions (40) along with (39) after making use of identities (5) and (6), (for the case images), lead to the following system of equations in the unknowns images:

images

where images and images

The determinant of the above system is different from zero, hence images can be uniquely determined to give

images

and hence the basis functions images take the form:

images

Remark 4.1. If x in (39) is replaced by images, then it is easy to see that the basis functions take the form

images

and from (43), we see that they have the following explicit expression

images

satisfies the following boundary conditions

images

Remark 4.2. If we choose the basis functions in terms of the general parameters normalized Jacobi polynomials images, that is

images

then, the resulting linear system to be solved to find the coefficients Hm, k can not be found in simple forms free of hypergeometric forms except for special choices of their parameters. The choice images in this paper enables one to get the simple form of the basis functions in (43).

Now, the basis functions images are linearly independent, and hence we have

images

Now, denote

images

Then Eq. (38) can be transformed into the following matrix system

images

and the entries of the matrices An and Ri, n, images, are given explicitly in the following theorem.

Theorem 4.1. Let the family of basis functions images be as taken in (44). Setting images, images, images, images. Then

images

and the nonzero elements of the matrices images and R0, n are given explicitly by:

images

images

images

images

images

images

where

images

and images, images and zm, k are defined in (31), (9) and (42), respectively.

Proof. The nonzero elements of the matrices An, images can be obtained with the aid of the derivatives formula (19). Substituting q = 2n in relation (30), yields

images

where images is defined in (31). In virtue of the orthogonality relation (8), one can show for images, that

images

and, accordingly, relation (44) enables one to express the elements ank, j in the form

images

In particular, it is easy to see that the diagonal elements of An can be simplified to give

images

Regrading the elements of the matrices Ri, n, images, we have

images

Using relation (30) and the orthogonality relation (8), ri, nkj can be computed to give formulas (50)(52). Finally, regarding the elements of the matrix R0, n, we have

images

using the orthogonality relation (8), it can be shown that, for j = k + p, images,

images

This ends the proof of Theorem 4.1.

The special case corresponding to the choice: images, images, is of interest. In such case, the system in (47) reduces to a linear system with nonsingular upper triangular matrix for all values of n. The following corollary exhibits this result.

Corollary 4.1. For the case corresponds to images, images, the system in (47) reduces to images, where An is an upper triangular matrix, and its solution can be given explicitly as

images

where ankk and ankj are given, respectively, by (48) and (49).

Remark 4.3. The following even-order BVP

images

governed by the nonhomogeneous boundary conditions:

images

can be transformed to a modified one governed by homogeneous boundary conditions (see [43]).

5  A Collocation Algorithm for Treating Nonlinear BVPs

In this section, we implement and present a numerical algorithm based on a collocation method to solve even-order nonlinear BVPs. The basic idea is based on employing the operational matrix of derivatives of the shifted Jacobi polynomials images. First, the following corollary is essential in the sequel.

Corollary 5.1. The first derivative of the polynomials images can be written explicitly as

images

where

images

Proof. Direct from Eq. (30), setting q = 1.

Now, if we define the following vector

images

then, it is easy to see that

images

where M is the operational matrix of derivatives whose entries are defined as:

images

Now, consider the following nonlinear even-order BVPs:

images

subject to the first-kind of boundary conditions:

images

or the second-kind of boundary conditions:

images

Now, consider an approximate solution of the form

images

This approximate solution can be expressed as

images

where

images

and images is defined in (63).

In virtue of Eq. (64), it is easy to express the imagesth-derivative of uN(x). In fact, we can write

images

Based on relation (71), the residual of Eq. (65) can be written as

images

Now, for the sake of obtaining a numerical solution of (65), a collocation method is applied. More precisely, the residual in (72) is enforced to vanish at suitable collocation points. We choose them to be the distinct zeros of the polynomial images, that is we have

images

Furthermore, the boundary conditions (66) yield the following (2n) equations:

images

images

while, the boundary conditions (67) yield the following (2n) equations:

images

images

Merging the (N + 1) equations in (73)(75), or the (N + 1) equations in (73), (76) and (77), we get a nonlinear system in the expansion coefficients ak which can be solved with the aid of Newton’s iterative method. Thus an approximate solution uN(x) can be found.

6  Connection Formula between the Polynomials images and the Shifted Legendre Polynomials

This section is confined to present the connection formula between the polynomials images and the shifted Legendre polynomials on [a, b]. This connection formula will be of important interest in investigating the convergence and error analysis in the subsequent section. To this end, the following theorem which links between two different parameters Jacobi polynomials is needed.

Theorem 6.1. For every nonnegative integer n, the following connection formula holds:

images

Proof. The connection formula between two different parameters classical Jacobi polynomials is given by [44]

images

with

images

Taking into consideration relation (2), the connection formula (78) can be obtained.

Lemma 6.1. Let m and n be two nonnegative integers. The following reduction formula holds

images

Proof. If we set

images

then, with the aid of Zeilberger’s algorithm, it can be shown that the following second-order recurrence relation is satisfied by images

images

with the initial values

images

The recurrence relation (82) can be solved exactly to give

images

and, therefore, the reduction formula (81) can be followed.

Now, we are going to state and prove an important theorem that concerns the connection formula between images and the shifted Legendre polynomials P*k(x).

Theorem 6.2. The connection formula between the shifted normalized Jacobi polynomials images and the shifted Legendre polynomials P*k(x) can be written explicitly in the form

images

where the connection coefficients images are given by

images

Proof. To prove the connection formula (83) with coefficients in (84), we prove first the formula on [ −1, 1], that is, we prove the formula

images

The last relation can be written alternatively in the form

images

Setting images and images in Eq. (78) yields the following formula

images

images

If we make use of Lemma 6.1, then the last formula, after some computations can be transformed into the form in (86), and therefore, relation (85) is proved. Replacing x by images in (85) gives (83). Theorem 6.2 is proved.

7  Study of the Convergence and Error Analysis

In this section, we state and prove some results concerning the convergence analysis of the proposed expansions in Sections 4 and 5. First, we mention some theoretical prerequisites. The following lemmas and theorems are needed. Let us agree on the following notation: For two positive sequences ai, bi, we say that images if there exists a generic constant c such that images.

7.1 Prerequisites

Lemma 7.1. Reference [45] Let Iqk(x) denote the q-times repeated integrals of shifted Legendre polynomials, that is

images

then, Iqk(x) can be expressed in terms of shifted Legendre polynomials by the formula

images

where images is a polynomial of degree at most (q −1), and the coefficients Gm, k, q are given explicitly by

images

Lemma 7.2. Let Iqk(x) be defined as in the above lemma. The following estimate holds

images

Proof. The lemma can be proved with the aid of the repeated integrals formula in Lemma 7.1 and noting the well-known inequality: images.

Lemma 7.3. The connection coefficients in the connection problem (83), which expressed explicitly in (84) satisfy the following estimate

images

Proof. The result can be easily followed from the explicit formula of the connection coefficients in (84).

Lemma 7.4. For all images, the following estimate is valid

images

Proof. From the connection formula (83), and noting that images, we get

images

We will prove (89). We consider the proof if n is even. The case if n is odd can be treated similarly. From (90), we can write

images

with

images

With the aid of Zeilberger’s algorithm, it is not difficult to show that

images

and

images

and accordingly, we get

images

The proof is now complete.

7.2 Convergence Analysis

Now, we are ready to state and prove some lemmas and theorems concerning the convergence of the truncated approximate solutions and the estimate for the truncation error of the proposed expansions in Sections 4 and 5. We first prove our results for the proposed expansion in Section 5. Assume that the exact and approximate expansions are respectively as follows:

images

and

images

From now on, and for more convenience, let us agree on the following notations:

images

Theorem 7.1. If the exact solution in (92) is Cq-function, for some q > 3, images and for images, then the following statements are correct:

  1.  The expansion coefficients an satisfy the following estimate:

images

  2.  The series in (92) converges uniformly.

Proof. Due to (8), the expansion coefficients an are given by

images

It is clear from relation (9) that images, and consequently

images

Thanks to the connection formula (83) that links images with the shifted Legendre polynomials, we have

images

and in virtue of Lemma 7.3, the last inequality gives

images

Now, if we apply integration by parts q-times and along with the following result obtained by Leibniz’s Theorem

images

the following estimate can be obtained

images

Now, the direct application of Lemma 7.2, and after performing some manipulations, yield

images

By the hypothesis of the theorem, precisely, the boundedness of the derivatives u(r), we get

images

which completes the proof of the first-part of the theorem.

For the second part, we have

images

and with the aid of the first result of this theorem along with Lemma 7.4, we get

images

Taking into consideration the condition: images, and applying the integral test, the desired series converges uniformly, and this completes the proof of the second-part of the theorem.

Theorem 7.2. Under the hypotheses of Theorem 7.1, the following truncation error estimate is valid:

images

Proof. Using relations (92) and (93), it is clear that

images

Now, by the result of Theorem 7.1 along with Lemma 7.4 yield the following inequality:

images

Now, based on the asymptotic behaviour of the Riemann–Hurwitz function, we have

images

which completes the proof of the theorem.

Remark 7.1. The main idea to investigate the convergence and error analysis of the expansion that used in Section 4 is built on transforming the basis functions (45) into another suitable form, and after that perform similar procedures to those followed in Subsections 7.1 and 7.2. More precisely, the following theorem is the key for such study. It can be proved through induction.

Theorem 7.3. Let images be the basis functions in (45). The following identity holds for every positive integer n

images

8  Numerical Results

In this section, we give some numerical results obtained by using our two proposed algorithms, namely, shifted Jacobi Galerkin method (SJGM) and the shifted Jacobi collocation method (SJCM) which presented, respectively, in Sections 4 and 5. Furthermore, some comparisons between our two proposed algorithms with some other methods exist in the literature are presented hoping to show the efficiency and accuracy of our proposed methods. In this section, the symbol N refers to the number of retained modes. In addition, u(x) is the exact solution and uN(x) is the corresponding approximate solution and the error is denoted by E where images.

Example 1. Consider the sixth-order linear boundary value problem [46,47]

images

subject to the boundary conditions

images

The exact solution for this problem is: images.

Tab. 1 lists the maximum absolute errors E for various choices of N, c, and images. In Tabs. 2 and 3, comparisons between our proposed SJGM and the two methods developed in [46,47] are displayed. In addition, in Fig. 1, we depict the Log-error plot of Example 1 for the case corresponds to various values of N and images.

Table 1: Maximum absolute error of E for Example 1

images

Table 2: Errors at xi for N = 12 and c = 10 for Example 1

images

Table 3: Errors at xi for N = 12 and c = 1000 for Example 1

images

images

Figure 1: Log-error of Example 1 for c = 1

Remark 8.1. The results in Tab. 1 show that the case corresponding to third-kind Chebyshev polynomials images does not always give the best results.

Example 2. Consider the eighth-order linear boundary value problem [48]

images

subject to the boundary conditions

images

The exact solution for this problem is: u(x) = (1 − x)ex.

Tab. 4 presents the maximum absolute error E for some choices of images and N, while Tab. 5 lists a comparison between our proposed algorithm with the absolute errors obtained by the method proposed in [48]. In Fig. 2, we depict the absolute error plot of Example 2 for the case corresponds to N = 14 and various values images.

Table 4: Maximum absolute error of E for Example 2

images

Table 5: Errors at xi for N = 14 for Example 2

images

images

Figure 2: Absolute error plot of Example 2

Example 3. Consider the eighth-order nonlinear boundary value problem [49,50]

images

subject to the boundary conditions

images

The exact solution for this problem is: images.

Tab. 6 lists the maximum absolute error E for various values of N and images. For the sake of comparison with other methods, we compare the best errors obtained from the application of our algorithm with those obtained by the two methods developed in [49,50]. Tab. 7 displays these results.

Example 4. Consider the twelfth-order nonlinear boundary value problem [51,52]

images

subject to the boundary conditions:

images

The exact solution for this problem is: u(x) = 2 ex.

Tab. 8 presents the maximum absolute error E for some choices of images and N, while Tab. 9 lists a comparison between our proposed algorithm for N = 21 for various values of images with the best errors obtained by the methods developed in [51,52]. In Fig. 3, we depict the absolute error plot of Example 4 for the case corresponds to N = 21 and various values images.

Table 6: Maximum absolute error of E for Example 3

images

Table 7: Best errors for N = 40 for Example 3

images

Table 8: Maximum absolute error of E for Example 4

images

Table 9: Best errors for Example 4

images

images

Figure 3: Absolute error plot of Example 4

9  Conclusions

Some direct algorithms for treating both linear and nonlinear BVPs were analyzed and presented. The proposed solutions are spectral. Two different approaches were presented for solving such problems. In the linear case, a Galerkin approach is followed. The basis functions are expressed in terms of certain Jacobi polynomials which generalize the well-known Chebyshev polynomials of the third- and fourth-kinds. Another collocation approach is followed for the treatment of both linear and nonlinear BVPs. A new operational matrix of derivatives of certain shifted Jacobi polynomials was constructed for this purpose. The presented numerical results show that the expansion of the third-kind of Chebyshev is not always the best in approximation. Moreover, the numerical results show that our algorithms are applicable and of high accuracy.

Acknowledgement: The authors would like to thank the anonymous reviewers for carefully reading the article and also for their constructive and valuable comments which have improved the paper in its present form.

Funding Statement: The authors received no specific funding for this study.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

References

 1.  Ghimire, B. K., Li, X., Chen, C. S., Lamichhane, A. R. (2020). Hybrid Chebyshev polynomial scheme for solving elliptic partial differential equations. Journal of Computational and Applied Mathematics, 364, 112324. DOI 10.1016/j.cam.2019.06.040. [Google Scholar] [CrossRef]

 2.  Abd-Elhameed, W. M., Youssri, Y. H. (2019). Explicit shifted second-kind Chebyshev spectral treatment for fractional Riccati differential equation. Computer Modeling in Engineering & Sciences, 121(3), 1029–1049. DOI 10.32604/cmes.2019.08378. [Google Scholar] [CrossRef]

 3.  Hassani, H., Tenreiro Machado, J. A., Naraghirad, E. (2019). Generalized shifted Chebyshev polynomials for fractional optimal control problems. Communications in Nonlinear Science and Numerical Simulation, 75, 50–61. DOI 10.1016/j.cnsns.2019.03.013. [Google Scholar] [CrossRef]

 4.  Heydari, M. H., Avazzadeh, Z. (2020). Chebyshev–Gauss–Lobatto collocation method for variable-order time fractional generalized Hirota–Satsuma coupled KdV system. Engineering with Computers. DOI 10.1007/s00366-020-01125-5. [Google Scholar] [CrossRef]

 5.  Doha, E. H., Abd-Elhameed, W. M. (2014). On the coefficients of integrated expansions and integrals of Chebyshev polynomials of third and fourth kinds. Bulletin of the Malaysian Mathematical Sciences Society, 37(2), 383–398. [Google Scholar]

 6.  Doha, E. H., Abd-Elhameed, W. M., Bassuony, M. A. (2015). On the coefficients of differentiated expansions and derivatives of Chebyshev polynomials of the third and fourth kinds. Acta Mathematica Scientia, 35(2), 326–338. DOI 10.1016/S0252-9602(15)60004-2. [Google Scholar] [CrossRef]

 7.  Sakran, M. R. A. (2019). Numerical solutions of integral and integro-differential equations using Chebyshev polynomials of the third kind. Applied Mathematics and Computation, 351, 66–82. DOI 10.1016/j.amc.2019.01.030. [Google Scholar] [CrossRef]

 8.  Masjed-Jamei, M. (2006). Some new classes of orthogonal polynomials and special functions: A symmetric generalization of Sturm-Liouville problems and its consequences, Ph.D. Thesis, University of Kassel, Kassel. [Google Scholar]

 9.  Abd-Elhameed, W. M., Youssri, Y. H. (2018). Fifth-kind orthonormal Chebyshev polynomial solutions for fractional differential equations. Computational and Applied Mathematics, 37(3), 2897–2921. DOI 10.1007/s40314-017-0488-z. [Google Scholar] [CrossRef]

10. Abd-Elhameed, W. M., Youssri, Y. H. (2019). Sixth-kind Chebyshev spectral approach for solving fractional differential equations. International Journal of Nonlinear Sciences and Numerical Simulation, 20(2), 191–203. DOI 10.1515/ijnsns-2018-0118. [Google Scholar] [CrossRef]

11. Babaei, A., Jafari, H., Banihashemi, S. (2020). Numerical solution of variable order fractional nonlinear quadratic integro-differential equations based on the sixth-kind Chebyshev collocation method. Journal of Computational and Applied Mathematics, 377, 112908. DOI 10.1016/j.cam.2020.112908. [Google Scholar] [CrossRef]

12. Doha, E. H., Hafez, R. M., Youssri, Y. H. (2019). Shifted Jacobi spectral-Galerkin method for solving hyperbolic partial differential equations. Computers and Mathematics with Applications, 78(3), 889–904. DOI 10.1016/j.camwa.2019.03.011. [Google Scholar] [CrossRef]

13. Lin, J., Reutskiy, S. (2020). A cubic B-spline semi-analytical algorithm for simulation of 3D steady-state convection-diffusion-reaction problems. Applied Mathematics and Computation, 371, 124944. DOI 10.1016/j.amc.2019.124944. [Google Scholar] [CrossRef]

14. Tian, H., Zhang, J., Ju, L. (2020). A spectral collocation method for nonlocal diffusion equations with volume constrained boundary conditions. Applied Mathematics and Computation, 370, 124930. DOI 10.1016/j.amc.2019.124930. [Google Scholar] [CrossRef]

15. Napoli, A., Abd-Elhameed, W. M. (2018). Numerical solution of eighth-order boundary value problems by using Legendre polynomials. International Journal of Computational Methods, 15(2), 1750083. DOI 10.1142/S0219876217500839. [Google Scholar] [CrossRef]

16. Atkinson, K., Chien, D., Hansen, O. (2019). Spectral methods using multivariate polynomials on the unit ball. Boca Raton: CRC Press. [Google Scholar]

17. Canuto, C., Hussaini, M. Y., Quarteroni, A., Zang, T. A. (1988). Spectral methods in fluid dynamics. New York: Springer-Verlag. [Google Scholar]

18. Shizgal, B. (2015). Spectral methods in chemistry and physics: Applications to kinetic theory and quantum mechanics. Dordrecht: Springer. [Google Scholar]

19. Hesthaven, J. S., Gottlieb, S., Gottlieb, D. (2007). Spectral methods for time-dependent problems21. Cambridge: Cambridge University Press. [Google Scholar]

20. Boyd, J. P. (2001). Chebyshev and Fourier spectral methods. North Chelmsford: Courier Corporation. [Google Scholar]

21. Trefethen, L. N. (2000). Spectral methods in MATLAB, 10. Philadelphia: SIAM. [Google Scholar]

22. Agarwal, R. P. (1986). Boundary value problems for higher ordinary differential equations. Singapore: World Scientific. [Google Scholar]

23. Tomar, S. (2020). A computationally efficient iterative scheme for solving fourth-order boundary value problems. International Journal of Applied and Computational Mathematics, 6(4), 1–16. 111. DOI 10.1007/s40819-020-00864-9. [Google Scholar] [CrossRef]

24. Baldwin, P. (1987). Asymptotic estimates of the eigenvalues of a sixth-order boundary-value problem obtained by using global phase-integral methods. Philosophical Transactions of the Royal Society of London, Series A, Mathematical and Physical Sciences, 322(1566), 281–305. DOI 10.1098/rsta.1987.0051. [Google Scholar] [CrossRef]

25. Bishop, R. E. D., Cannon, S. M., Miao, S. (1989). On coupled bending and torsional vibration of uniform beams. Journal of Sound and Vibration, 131(3), 457–464. DOI 10.1016/0022-460X(89)91005-5. [Google Scholar] [CrossRef]

26. Islam, S. U., Haq, S., Ali, J. (2009). Numerical solution of special 12th-order boundary value problems using differential transform method. Communications in Nonlinear Science and Numerical Simulation, 14(4), 1132–1138. DOI 10.1016/j.cnsns.2008.02.012. [Google Scholar] [CrossRef]

27. Chandrasekhar, S. (2013). Hydrodynamic and hydromagnetic stability. Dover Publications, INC. New York: Courier Corporation. [Google Scholar]

28. Davies, A. R., Karageorghis, A., Phillips, T. N. (1988). Spectral Galerkin methods for the primary two-point boundary value problem in modelling viscoelastic flows. International Journal of Numerical Methods in Engineering, 26(3), 647–662. DOI 10.1002/nme.1620260309. [Google Scholar] [CrossRef]

29. Paliwal, D. N., Pande, A. (1999). Orthotropic cylindrical pressure vessels under line load. International Journal of Pressure Vessels and Piping, 76(7), 455–459. DOI 10.1016/S0308-0161(99)00010-1. [Google Scholar] [CrossRef]

30. Siddiqi, S. S., Twizell, E. H. (1996). Spline solutions of linear eighth-order boundary-value problems. Computer Methods in Applied Mechanics and Engineering, 131(3), 309–325. DOI 10.1016/0045-7825(96)88162-X. [Google Scholar] [CrossRef]

31. Twizell, E. H., Boutayeb, A., Djidjeli, K. (1994). Numerical methods for eighth-, tenth- and twelfth-order eigenvalue problems arising in thermal instability. Advances in Computational Mathematics, 2(4), 407–436. DOI 10.1007/BF02521607. [Google Scholar] [CrossRef]

32. He, J. H. (2003). Homotopy perturbation method: A new nonlinear analytical technique. Applied Mathematics and Computation, 135(1), 73–79. DOI 10.1016/S0096-3003(01)00312-5. [Google Scholar] [CrossRef]

33. Costabile, F. A., Napoli, A. (2015). Collocation for high order differential equations with two-points Hermite boundary conditions. Applied Numerical Mathematics, 87, 157–167. DOI 10.1016/j.apnum.2014.09.008. [Google Scholar] [CrossRef]

34. Doha, E. H., Bhrawy, A. H., Saker, M. A. (2011). Integrals of Bernstein polynomials: An application for the solution of high even-order differential equations. Applied Mathematics Letters, 24(4), 559–565. DOI 10.1016/j.aml.2010.11.013. [Google Scholar] [CrossRef]

35. Costabile, F. A., Napoli, A. (2016). A new spectral method for a class of linear boundary value problems. Journal of Computational and Applied Mathematics, 292, 329–341. DOI 10.1016/j.cam.2015.07.015. [Google Scholar] [CrossRef]

36. Abd-Elhameed, W. M. (2014). On solving linear and nonlinear sixth-order two point boundary value problems via an elegant harmonic numbers operational matrix of derivatives. Computer Modeling in Engineering & Sciences, 101(3), 159–185. DOI 10.3970/cmes.2014.101.159. [Google Scholar] [CrossRef]

37. Doha, E. H., Abd-Elhameed, W. M., Bhrawy, A. H. (2013). New spectral- Galerkin algorithms for direct solution of high even-order differential equations using symmetric generalized Jacobi polynomials. Collectanea Mathematica, 64(3), 373–394. DOI 10.1007/s13348-012-0067-y. [Google Scholar] [CrossRef]

38. Doha, E. H., Abd-Elhameed, W. M., Ahmed, H. M. (2012). The coefficients of differentiated expansions of double and triple Jacobi polynomials. Bulletin of the Iranian Mathematical Society, 38(3), 739–765. [Google Scholar]

39. Koepf, W. (2014). Hypergeometric summation. London: Springer-Verlag. [Google Scholar]

40. Shen, J. (1994). Efficient spectral-Galerkin method I. Direct solvers of second- and fourth-order equations using Legendre polynomials. SIAM Journal on Scientific Computing, 15(6), 1489–1505. DOI 10.1137/0915089. [Google Scholar] [CrossRef]

41. Shen, J. (1995). Efficient spectral-Galerkin method II. Direct solvers of second-and fourth-order equations using Chebyshev polynomials. SIAM Journal on Scientific Computing, 16(1), 74–87. DOI 10.1137/0916006. [Google Scholar] [CrossRef]

42. Doha, E. H., Abd-Elhameed, W. M., Bhrawy, A. H. (2009). Efficient spectral ultraspherical-Galerkin algorithms for the direct solution of 2nd-order linear differential equations. Applied Mathematical Modelling, 33(4), 1982–1996. DOI 10.1016/j.apm.2008.05.005. [Google Scholar] [CrossRef]

43. Doha, E. H., Abd-Elhameed, W. M., Bassuony, M. A. (2013). New algorithms for solving high even-order differential equations using third and fourth Chebyshev–Galerkin methods. Journal of Computational Physics, 236, 563–579. DOI 10.1016/j.jcp.2012.11.009. [Google Scholar] [CrossRef]

44. Andrews, G. E., Askey, R., Roy, R. (1999). Special functions, 71. Cambridge: Cambridge University Press. [Google Scholar]

45. Doha, E. H. (2002). On the coefficients of integrated expansions and integrals of ultraspherical polynomials and their applications for solving differential equations. Journal of Computational and Applied Mathematics, 139(2), 275–298. DOI 10.1016/S0377-0427(01)00420-4. [Google Scholar] [CrossRef]

46. Sohaib, M., Haq, S., Mukhtar, S., Khan, I. (2018). Numerical solution of sixth-order boundary-value problems using Legendre wavelet collocation method. Results in Physics, 8, 1204–1208. DOI 10.1016/j.rinp.2018.01.065. [Google Scholar] [CrossRef]

47. Noor, M. A., Mohyud-Din, S. T. (2008). Homotopy perturbation method for solving sixth-order boundary value problems. Computers & Mathematics with Applications, 55(12), 2953–2972. DOI 10.1016/j.camwa.2007.11.026. [Google Scholar] [CrossRef]

48. Golbabai, A., Javidi, M. (2007). Application of homotopy perturbation method for solving eighth-order boundary value problems. Applied Mathematics and Computation, 191(2), 334–346. DOI 10.1016/j.amc.2007.02.091. [Google Scholar] [CrossRef]

49. El-Gamel, M., Abdrabou, A. (2019). Sinc-Galerkin solution to eighth-order boundary value problems. SeMA Journal, 76(2), 249–270. DOI 10.1007/s40324-018-0172-2. [Google Scholar] [CrossRef]

50. Ballem, S., Viswanadham, K. K. N. S. (2015). Numerical solution of eighth order boundary value problems by Galerkin method with septic B-splines. Procedia Engineering, 127, 1370–1377. DOI 10.1016/j.proeng.2015.11.496. [Google Scholar] [CrossRef]

51. Noor, M. A., Mohyud-Din, S. T. (2010). Variational iteration method for solving twelfth-order boundary-value problems using He’s polynomials. Computational Mathematics and Modeling, 21(2), 239–251. DOI 10.1007/s10598-010-9068-4. [Google Scholar] [CrossRef]

52. Noor, M. A., Mohyud-Din, S. T. (2008). Solution of twelfth-order boundary value problems by variational iteration technique. Journal of Applied Mathematics and Computing, 28(1–2), 123–131. DOI 10.1007/s12190-008-0081-0. [Google Scholar] [CrossRef]

images This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.