Research Article  Open Access
A. H. Bhrawy, A. S. Alofi, S. I. ElSoubhy, "An Extension of the LegendreGalerkin Method for Solving SixthOrder Differential Equations with Variable Polynomial Coefficients", Mathematical Problems in Engineering, vol. 2012, Article ID 896575, 13 pages, 2012. https://doi.org/10.1155/2012/896575
An Extension of the LegendreGalerkin Method for Solving SixthOrder Differential Equations with Variable Polynomial Coefficients
Abstract
We extend the application of LegendreGalerkin algorithms for sixthorder elliptic problems with constant coefficients to sixthorder elliptic equations with variable polynomial coefficients. The complexities of the algorithm are O(N) operations for a onedimensional domain with () unknowns. An efficient and accurate direct solution for algorithms based on the LegendreGalerkin approximations developed for the twodimensional sixthorder elliptic equations with variable coefficients relies upon a tensor product process. The proposed LegendreGalerkin method for solving variable coefficients problem is more efficient than pseudospectral method. Numerical examples are considered aiming to demonstrate the validity and applicability of the proposed techniques.
1. Introduction
Spectral methods are preferable in numerical solutions of ordinary and partial differential equations due to its highorder accuracy whenever it works [1, 2]. Recently, renewed interest in the Galerkin technique has been prompted by the decisive work of Shen [3], where new Legendre polynomial bases for which the matrices systems are sparse are introduced. We introduce a generalization of Shen's basis to numerically solve the sixthorder differential equations with variable polynomial coefficients.
Sixthorder boundaryvalue problems arise in astrophysics; the narrow convecting layers bounded by stable layers, which are believed to surround Atype stars, may be modeled by sixthorder boundaryvalue problems [4, 5]. Further discussion of the sixthorder boundaryvalue problems is given in [6]. The literature of numerical analysis contains little work on the solution of the sixthorder boundaryvalue problems [4, 5, 7, 8]. Theorems that list conditions for the existence and uniqueness of solutions of such problems are thoroughly discussed in [9], but no numerical methods are contained therein.
From the numerical point of view, Shen [3], Doha and Bhrawy [10–12], and Doha et al. [13] have constructed efficient spectralGalerkin algorithms using compact combinations of orthogonal polynomials for solving elliptic equations of the second and fourth order with constant coefficients in various situations. Recently, the authors in [14, 15] and [16] have developed efficient Jacobi dualPetrovGalerkin and JacobiGauss collocation methods for solving some oddorder differential equations. Moreover, the Bernstein polynomials have been applied for the numerical solution of high evenorder differential equations (see, [17, 18]).
For sixthorder differential equations, Twizell and Boutayeb [5] developed finitedifference methods of order two, four, six, and eight for solving such problems. Siddiqi and Twizell [7] used sixthdegree splines, where spline values at the mid knots of the interpolation interval and the corresponding values of the even order derivatives were related through consistency relations. A sixthdegree Bspline functions is used to construct an approximate solution for sixthorder boundaryvalue problems (see [19]). Moreover, Septic spline solutions of sixthorder boundary value problems are introduced in [20]. ElGamel et al. [8] proposed SincGalerkin method for the solutions of sixthorder boundaryvalue problems. In fact, the decomposition and modified domain decomposition methods to investigate solution of the sixthorder boundaryvalue problems are introduced in [21]. Recently, Bhrawy [22] developed a spectral LegendreGalerkin method for solving sixthorder boundaryvalue problems with constant coefficients. In this work, we introduce an efficient direct solution algorithm to generalize the work in [3, 22].
The main aim of this paper is to extend the application of LegendreGalerkin method (LGM) to solve sixthorder elliptic differential equations with variable coefficients by using the expansion coefficients of the moments of the Legendre polynomials and their highorder derivatives. We present appropriate basis functions for the LegendreGalerkin method applied to these equations. This leads to discrete systems with sparse matrices that can be efficiently inverted. The complexities of the algorithm is operations for a onedimensional domain with unknowns. The direct solution algorithms developed for the homogeneous problem in twodimensions with constant and variable coefficients rely upon a tensor product process. Numerical results indicating the high accuracy and effectiveness of these algorithms are presented.
This paper is organized as follows. In the next section, we discuss an algorithm for solving the onedimensional sixthorder elliptic equations with variable polynomial coefficients. In Section 3, we extend our results of Sections 2 to the twodimensional sixthorder equations with variable polynomial coefficients. In Section 4, we present two numerical examples to exhibit the accuracy and efficiency of the proposed numerical algorithms. Also a conclusion is given in Section 5.
2. OneDimensional SixthOrder Equations with Polynomial Coefficients
We first introduce some basic notation which will be used in the sequel. We denote by the th degree Legendre polynomial, and we set where denotes thorder differentiation of with respect to .
We recall that the satisfy the orthogonality relation We recall also that is a polynomial of degree and therefore . The following relation (the th derivative of ) will be needed for our main results (see Doha [23]) where
Some other useful relations are
In this section, we are interested in using the LegendreGalerkin method to solve the variable polynomial coefficients sixthorder differential equation in the form: subject to where , are constants and , are given polynomials. Moreover, is a given source function. Without loss of generality, we suppose that , and where , , and are positive integers.
2.1. Basis of Functions
The problem of approximating solutions of ordinary or partial differential equations by Galerkin approximation involves the projection onto the span of some appropriate set of basis functions, typically arising as the eigenfunctions of a singular SturmLiouville problem. The members of the basis may satisfy automatically the boundary conditions imposed on the problem. As suggested in [3, 10–12], one should choose compact combinations of orthogonal polynomials as basis functions to minimize the bandwidth and condition number of the resulting system. As a general rule, for onedimensional sixthorder differential equations with six boundary conditions, one can choose the basis functions of expansion of the form We will choose the coefficients such that verifies the boundary conditions (2.7). Making use of (2.5) and (2.8), hence can be uniquely determined to obtain where . Now, substitution of (2.9) into (2.8) yields It is obvious that are linearly independent. Therefore by dimension argument we have
2.2. Treatment of Variable Polynomial Coefficients
A more general situation which often arises in the numerical solution of differential equations with polynomial coefficients by using the Legendre Galerkin method is the evaluation of the expansion coefficients of the moments of highorder derivatives of infinitely differentiable functions. The formula of Legendre coefficients of the moments of one single Legendre polynomials of any degree is with , where For more details about the above formula, the reader is referred to Doha [23]. This formula can be used to facilitate greatly the setting up of the algebraic systems to be obtained by applying the LGM for solving differential equations with polynomial coefficients of any order. The following lemma is very important and needed in what follows.
Lemma 2.1. We have, for arbitrary constants , where , , and are as defined in (2.8), (2.9), and (2.13), respectively.
Proof. Immediately obtained from relations (2.3), (2.8), and (2.12), the LegendreGalerkin approximation to (2.6)(2.7) is, to find such that
where is the scalar product in and is the inner product associated with the LegendreGaussLobatto quadrature. It is clear that if we take as defined in (2.8) and , then we find that (2.15) is equivalent to
Hence, by setting
where
then the matrix system associated with (2.16) becomes
where the nonzero elements of the matrices , , , and are given explicitly in the following theorem.
Theorem 2.2. If we take as defined in (2.8), and if we denote , , and then the nonzero elements , , for , are given by
Proof. The proof of this theorem is rather lengthy, but it is not difficult once Lemma 2.1 is applied.
From Theorem 2.2, we see that is a band matrix with an upper bandwidth of , lower bandwidth of , and an overall bandwidth . The sparse matrices , , and have bandwidths of , , and , respectively.
In general, the expense of calculating an LU factorization of an dense matrix is operations, and the expense of solving , provided that the factorization is known, is . However, in the case that a banded has bandwidth of , we need just operations to factorize and operations to solve a linear system. In the case of , , the square matrix has bandwidth of . We need just operations to factorize and operations to solve the linear system (2.19). If this represents a very substantial saving. Notice also that the system (2.19) reduces to a diagonal system for and , .
2.3. Constant Coefficients
In the special case, (, i.e., the sixthorder differential equation with constant coefficients), the corresponding matrix system becomes , where ; .
Corollary 2.3. If then the nonzero elements , , , for , are given as follows:
Note that the results of Corollary 2.3 can be obtained immediately as a special case from Theorem 2.2. For more details see [22].
It is worthy to note here that if , , then the nonzero elements of the matrix are given by (2.21) and the solution of the linear system is given explicitly by .
Obviously , , and are symmetric positive definite matrices. Furthermore, is a diagonal matrix, can be split into two tridiagonal submatrices, can be split into two pentadiagonal submatrices, and can be split into two sparse submatrices with bandwidth of 7. Therefore, the system can be efficiently solved. More precisely for odd, , . Hence system of order can be decoupled into two separate systems of order and , respectively. In this way one needs to solve two systems of order instead of one of order , which leads to substantial savings. Moreover, in the case of , , we can form explicitly the LU factorization, that is, . The special structure of and allows us to obtain the solution in operations.
Remark 2.4. If the boundary conditions are nonhomogeneous, one can split the solution into the sum of a lowdegree polynomial which satisfies the nonhomogeneous boundary conditions plus a sum over the basis functions that satisfy the equivalent homogeneous boundary conditions.
3. TwoDimensional SixthOrder Equations with Polynomial Coefficients
In this section, we extend the results of Section 2 to deal with the twodimensional sixthorder differential equations with variable polynomial coefficients: subject to the boundary conditions where , the differential operator is the wellknown Laplacian defined by , , and is a given source function. Moreover, and , are given polynomials. Without loss of generality, we suppose that , , , , , , and where , , , , , and are positive integers.
The LegendreGalerkin approximation to (3.1)(3.2) is, to find such that It is clear that if we take as defined in (2.8), then
We denote Taking in (3.3) for , then one can observe that (3.3) is equivalent to the following equation: which may be written in the matrix form where are the matrices defined in Theorem 2.2.
The direct solution algorithm here developed for the sixthorder elliptic differential equation in two dimensions relies upon a tensor product process, which is defined as follows. Let P and R be two matrices of size and , respectively. Their tensor product is a matrix of size .
We can also rewrite (3.7) in the following form using the Kronecker matrix algebra (See, Graham [24]): where and are and written in a column vector, that is, and denotes the tensor product of matrices, that is, . In brief, the solution of (3.1) subject to (3.2) can be summarized in Algorithm 1.
4. Numerical Results
We report on two numerical examples by using the algorithms presented in the previous sections. It is worthy to mention that the pure spectralGalerkin technique is rarely used in practice, since for a general righthand side function we are not able to compute exactly its representation by Legendre polynomials. In fact, the socalled pseudospectral method is used to treat the righthand side; that is, we replace by (polynomial interpolation over the set of GaussLobatto points); see for instance; Funaro [25].
Example 4.1. Here we present some numerical results for a onedimensional sixthorder equation with polynomial coefficients. We only consider the case , , and ; that is, we consider the subject to where is chosen such that the exact solution of (4.1) is .
For , we have , the vectors of unknowns is the solution of the following system where the nonzero elements of the matrices , , , , , , and can be evaluated explicitly from Theorem 2.2. Table 1 lists the maximum pointwise error of , using the LGM with various choices of and . Numerical results of this problem show that the Legendre Galerkin method converge exponentially.

In order to examine the algorithm proposed in Section 3, we will consider a problem for a twodimensional sixthorder elliptic differential equation with constant coefficients.
Example 4.2. Consider the twodimensional sixthorder equation subject to the boundary conditions where is chosen such that the exact solution of (4.4)(4.5) is
In Table 2, we list the maximum pointwise error of by the LGM with two choices of , , and various choices of . The results indicate that the spectral accuracy is achieved and that the effect of roundoff errors is very limited.

5. Concluding Remarks
We have presented stable and efficient spectral Galerkin method using Legendre polynomials as basis functions for sixthorder elliptic differential equations in one and two dimensions. We concentrated on applying our algorithms to solve variable polynomials coefficients differential equations by using the expansion coefficients of the moments of the Legendre polynomials and their highorder derivatives. Numerical results are presented which exhibit the high accuracy of the proposed algorithms.
Acknowledgments
The authors are very grateful to the referees for carefully reading the paper and for their comments and suggestions which have improved the paper. This paper was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah. The authors, therefore, acknowledge with thanks DSR technical and financial support.
References
 C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Mechanics, Springer, New York, NY, USA, 2006.
 C. Bernardi and Y. Maday, Approximations Spectrales de Problèmes aux Limites Elliptiques, vol. 10, Springer, Paris, France, 1992.
 J. Shen, “Efficient spectralGalerkin method. I. Direct solvers of second and fourthorder equations using Legendre polynomials,” SIAM Journal on Scientific Computing, vol. 15, no. 6, pp. 1489–1505, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Boutayeb and E. Twizell, “Numerical methods for the solution of special sixthorder boundary value problems,” International Journal of Computer Mathematics, vol. 45, pp. 207–233, 1992. View at: Google Scholar
 E. H. Twizell and A. Boutayeb, “Numerical methods for the solution of special and general sixthorder boundary value problems, with applications to Bénard layer eigenvalue problems,” Proceedings of the Royal Society London Series A, vol. 431, no. 1883, pp. 433–450, 1990. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. Baldwin, “Asymptotic estimates of the eigenvalues of a sixthorder boundaryvalue problem obtained by using global phaseintegral methods,” Philosophical Transactions of the Royal Society of London Series A, vol. 322, no. 1566, pp. 281–305, 1987. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. S. Siddiqi and E. H. Twizell, “Spline solutions of linear sixthorder boundaryvalue problems,” International Journal of Computer Mathematics, vol. 60, no. 34, pp. 295–304, 1996. View at: Google Scholar
 M. ElGamel, J. R. Cannon, and A. I. Zayed, “SincGalerkin method for solving linear sixthorder boundaryvalue problems,” Mathematics of Computation, vol. 73, no. 247, pp. 1325–1343, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. P. Agarwal, Boundary Value Problems for Higher Order Differential Equations, World Scientific Publishing, Singapore, 1986. View at: Zentralblatt MATH
 E. H. Doha and A. H. Bhrawy, “Efficient spectralGalerkin algorithms for direct solution for secondorder differential equations using Jacobi polynomials,” Numerical Algorithms, vol. 42, no. 2, pp. 137–164, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. H. Doha and A. H. Bhrawy, “Efficient spectralGalerkin algorithms for direct solution of fourthorder differential equations using Jacobi polynomials,” Applied Numerical Mathematics, vol. 58, no. 8, pp. 1224–1244, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. H. Doha and A. H. Bhrawy, “A Jacobi spectral Galerkin method for the integrated forms of fourthorder elliptic differential equations,” Numerical Methods for Partial Differential Equations, vol. 25, no. 3, pp. 712–739, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. H. Doha, A. H. Bhrawy, and W. M. AbdElhameed, “Jacobi spectral Galerkin method for elliptic Neumann problems,” Numerical Algorithms, vol. 50, no. 1, pp. 67–91, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “A JacobiJacobi dualPetrovGalerkin method for third and fifthorder differential equations,” Mathematical and Computer Modelling, vol. 53, no. 910, pp. 1820–1832, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “A Jacobi dualPetrovGalerkin method for solving some oddorder ordinary differential equations,” Abstract and Applied Analysis, vol. 2011, Article ID 947230, 21 pages, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. H. Bhrawy and W. M. AbdElhameed, “New algorithm for the numerical solutions of nonlinear thirdorder differential equations using JacobiGauss collocation method,” Mathematical Problems in Engineering, vol. 2011, Article ID 837218, 14 pages, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. H. Doha, A. H. Bhrawy, and M. A. Saker, “Integrals of Bernstein polynomials: an application for the solution of high evenorder differential equations,” Applied Mathematics Letters, vol. 24, no. 4, pp. 559–565, 2011. View at: Publisher Site  Google Scholar
 E. H. Doha, A. H. Bhrawy, and M. A. Saker, “On the derivatives of Bernstein polynomials: an application for the solution of high evenorder differential equations,” Boundary Value Problems, vol. 2011, Article ID 829543, 16 pages, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. B. Loghmani and M. Ahmadinia, “Numerical solution of sixth order boundary value problems with sixth degree Bspline functions,” Applied Mathematics and Computation, vol. 186, no. 2, pp. 992–999, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. S. Siddiqi and G. Akram, “Septic spline solutions of sixthorder boundary value problems,” Journal of Computational and Applied Mathematics, vol. 215, no. 1, pp. 288–301, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Wazwaz, “The numerical solution of sixthorder boundary value problems by the modified decomposition method,” Applied Mathematics and Computation, vol. 118, no. 23, pp. 311–325, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. H. Bhrawy, “LegendreGalerkin method for sixthorder boundary value problems,” Journal of the Egyptian Mathematical Society, vol. 17, no. 2, pp. 173–188, 2009. View at: Google Scholar  Zentralblatt MATH
 E. H. Doha, “On the construction of recurrence relations for the expansion and connection coefficients in series of Jacobi polynomials,” Journal of Physics A, vol. 37, no. 3, pp. 657–675, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Graham, Kronecker Products and Matrix Calculus: With Applications, Ellis Horwood Ltd., England, UK, 1981.
 D. Funaro, Polynomial Approximation of Differential Equations, vol. 8 of Lecture Notes in Physics, Springer, Berlin, Germany, 1992.
Copyright
Copyright © 2012 A. H. Bhrawy et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.