Research Article  Open Access
Ghazala Akram, Hamood Ur Rehman, "Solutions of a Class of Sixth Order Boundary Value Problems Using the Reproducing Kernel Space", Abstract and Applied Analysis, vol. 2013, Article ID 560590, 8 pages, 2013. https://doi.org/10.1155/2013/560590
Solutions of a Class of Sixth Order Boundary Value Problems Using the Reproducing Kernel Space
Abstract
The approximate solution to a class of sixth order boundary value problems is obtained using the reproducing kernel space method. The numerical procedure is applied on linear and nonlinear boundary value problems. The approach provides the solution in terms of a convergent series with easily computable components. The present method is simple from the computational point of view, resulting in speed and accuracy significant improvements in scientific and engineering applications.It was observed that the errors in absolute values are better than compared (Che Hussin and Kiliçman (2011) and, Noor and MahyudDin (2008), Wazwaz (2001), Pandey (2012)).Furthermore, the nonlinear boundary value problem for the integrodifferential equation has been investigated arising in chemical engineering, underground water flow and population dynamics, and other fields of physics and mathematical chemistry. The performance of reproducing kernel functions is shown to be very encouraging by experimental results.
1. Introduction
Boundary value problems arise in engineering, applied mathematics, and several branches of physics and have attracted much attention. However, it is difficult to obtain closedform solutions for boundary value problems, especially for nonlinear problems. In most cases, only approximate solutions (either numerical solutions or analytical solutions) can be expected.
It is well known that a wide class of boundary value problems arise in various branches of pure and applied sciences including astrophysics, structural engineering, optimization, and economics. The literature of the numerical solution of sixth order boundary value problems is sparse. These types of problems generally arise in astrophysics; the narrow convecting layers bounded by stable layers which are believed to surround Atype stars may be modeled by sixth order boundary value problems [1].
Chandrasekhar [2] determined that when an infinite horizontal layer of fluid is heated from below and is under the action of rotation, instability sets in. When this instability is an ordinary convection, the ordinary differential equation is a sixth order. Details of the theorems which listed the conditions for the existence and uniqueness of solutions of sixth order boundary value problems are given in Agarwal [3], but no numerical methods are contained therein.
Higher order boundary value problems arise in the study of astrophysics, hydrodynamic and hydromagnetic stability, fluid dynamics, astronomy, beam and long wave theory, engineering, and applied physics. The boundary value problems of higher order have been investigated due to their mathematical importance and the potential for applications in diversified applied sciences.
Siddiqi and Akram [4, 5] presented a second order method using polynomial and nonpolynomial septic splines for the solution of linear sixth order boundary value problems with boundary conditions at first and second order derivatives. Noor and MohyudDin [6] proposed a reliable approach for solving linear and nonlinear sixth order boundary value problems by homotopy perturbation method. Wazwaz [7] used decomposition and modified domain decomposition methods to investigate the solution of sixth order nonlinear boundary value problems by making a comparison among differential transformation method (DTM) and Adomian decomposition method (ADM) [8]. Simos [9] proposed a new closed NewtonCotes trigonometrically fitted differential method of high algebraic order. Simos [10] used a hybrid twostep method for the numerical solution of the radial Schrdinger equation and related problems with periodic or oscillating solutions. Anastassi and Simos [11] developed an efficient parametric symmetric linear phase fitted fourstep method for the numerical solution of the Schrdinger equation and related oscillatory problems. A reproducing kernel Hilbert space is a useful framework for constructing approximate solutions of boundary value problems [12–15]. In this paper, a reproducing kernel method is used for the solution of a class of sixth order boundary value problems. To the best of our knowledge, the sixth order nonlinear boundary value problem for the integrodifferential equation which arises in chemical engineering, underground water flow and population dynamics, and other fields of physics and mathematical chemistry has not been investigated yet. It is illustrated that our methods can also be used to find the approximate solution of such integrodifferential equations.
Consider the following sixth order twopoint boundary value problem (BVP): where and , are finite real constants and , and are the continuous functions on . Let the differential operator be , and homogenization of the boundary conditions of (1) can be transformed into the following form:
Therefore, to solve (1), it suffices to solve (2).
The rest of this paper is organized as follows. In Section 2, the reproducing kernel function is obtained, and reproducing kernel spaces needed in this paper are defined. After that in Section 3, a solution of (2) with initial boundary value conditions is presented. Finally, numerical examples are discussed to demonstrate the accuracy of the presented method in Section 4.
2. Reproducing Kernel Spaces
(i)The reproducing kernel space is defined by = , are absolutely continuous realvalued functions in , . The inner product in is given by (ii)The reproducing kernel space is defined by is absolutely continuous realvalued function in , . Also, the inner product is given by
Theorem 1. The space is a reproducing kernel Hilbert space. That is, for all and each fixed , there exists such that , and is called the reproducing kernel function of space . The reproducing kernel function in is given by
The coefficients are given in the appendix at the end.
3. The Exact and Approximate Solutions
The solution of (2) is given in the reproducing kernel Hilbert space , and the linear operator is bounded. Choose a countable dense subset in the domain , and let where is reproducing kernel of . Further assume that , where is the conjugate operator of .
Theorem 2. is a complete system of and .
Proof. For each fixed , let (), which implies that
Since is dense in , , which implies that from the existence of .
Using reproducing kernel property, it can be written as
To orthonormalize the sequence in the reproducing kernel space , GramSchmidt process can be used as
Theorem 3. If is dense in and the solution of (2) is unique, for all , the series is convergent in the norm of . If is exact solution; then the solution of (2) has the form
Proof. Since and can be expanded in the form of Fourier series about normal orthogonal system as
Since the space is Hilbert space, so the series is convergent in the norm of . From (9) and (11), it can be written as
If is the exact solution of (2) and , then
The approximate solution of is given by
Remark 4. If (2) is linear, that is, , then the solution can be obtained directly from (14).
If (2) is nonlinear, the approximate solution can be obtained using the following method.
Using (13), the following iterative formula can be constructed as
It can be noted that the approximate solution given in (15) satisfied the boundary conditions of problem (2). In fact, the solution of problem (2) is considered as the fixed point of the following functional under the suitable choice of the initial term :
Theorem 5 (Banach’s fixed point theorem). Assume that X is a Banach space and is a nonlinear function, and assume that for some constants . Then, A has a unique fixed point. Furthermore, the sequence with an arbitrary choice of converges to the fixed point of A.
According to Theorem 5, for the nonlinear mapping, a sufficient condition for convergence of the present iteration method is strictly contraction of A. Furthermore, the sequence in (15) converges to the fixed point of A which is also the solution of problem (2).
The approximate solution can be obtained by taking finitely many terms in the series representation of , given by
4. Numerical Examples
To illustrate the applicability and effectiveness of our method, three numerical examples are constructed. All the numerical computations are performed by using Mathematica 5.2 on a system with Intel Core 2 Quad 2.6 GHz CPU and 3 GB of RAM.
Example 6. Consider the following special sixth order boundary value problem involving a parameter c [6]:
The exact solution of problem (20) is A comparison of the errors in absolute values between the method developed in this paper and that of Noor and MohyudDin [6] is shown in Table 1. For small values of , , and , it is noted that the results obtained by our method are better than the method developed by Noor and MohyudDin [6]. The results obtained by Noor and MohyudDin using Adomian's decomposition method, variation iteration method, and Homotopy perturbation method for large values of are no longer valid as given in [6], but the results obtained by the present method for Example 6 are also valid and accurate for large value of .
Example 7. Consider the following nonlinear boundary value problem of sixth order [6, 7]:
The exact solution of problem (22) is A comparison of the errors in absolute values between the method developed in this paper and that of Wazwaz [7] and Noor and MohyudDin [6] is shown in Table 2 and Figures 1 and 2.
Example 8. Consider the following nonlinear boundary value problem of sixth order:
The exact solution of problem (24) is A comparison of the errors in absolute values between the method developed in this paper and that of Che Hussin and Kiliçman [8] is shown in Table 3 and Figures 3 and 4.

Example 9. Consider the following nonlinear integral boundary value problem of sixth order:
The exact solution of problem (24) is The absolute error (between exact and approximate solutions) and relative error are given in Table 4 and Figures 5 and 6.

Example 10. Consider the following nonlinear boundary value problem of sixth order:
The exact solution of problem (24) is
It is observed that the results obtained from our method which are obtained in the form of maximum absolute error (MAE) are better than the results obtained from the finite difference method [16] as shown in Table 5. It is noted that CPU time (s) for is 0.346, is 0.416, is 0.687, is 1.886, and is 4.366.

5. Conclusion
In this paper, a reproducing kernel space method to find the solution of a class of sixth order boundary value problems is considered. The properties of the reproducing kernel space require no more integral computation for some functions, instead of computing some values of a function at some nodes. This simplification of integral computation not only improves the computational speed but also improves the computational accuracy. It was observed that the errors in absolute values are better than the other developed methods [6–8, 16]. Mathematica software is used for all computational work. The numerical results show that only a few number of iteration steps can be used for numerical purpose with a high degree of accuracy. It is analyzed that our proposed method is well suited for use in higher order boundary value problems as it provides the best solution in a less number of iterations and reduces the computational work. Therefore, the present method is an accurate and reliable analytical technique for sixth order boundary value problems.
Appendix
+ , , + , + + + − , + − + − , + − + , , , , , , , , , + − + , , − + , + + + − , + − − − , + − + , , , , , , , , .References
 J. Toomre, J. R. Zahn, J. Latour, and E. A. Spiegel, “Stellar convection theory ii:single mode study of the second convection zone in atype stars,” The Astrophysical Journal, vol. 207, pp. 545–563, 1976. View at: Publisher Site  Google Scholar
 S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, The International Series of Monographs on Physics, Clarendon Press, Oxford, UK, 1961. View at: MathSciNet
 R. P. Agarwal, Boundary Value Problems for Higher Order Differential Equations, World Scientific Publishing, Singapore, 1986. View at: Zentralblatt MATH  MathSciNet
 G. Akram and S. S. Siddiqi, “Solution of sixth order boundary value problems using nonpolynomial spline technique,” Applied Mathematics and Computation, vol. 181, no. 1, pp. 708–720, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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  MathSciNet
 M. A. Noor and S. T. MohyudDin, “Homotopy perturbation method for solving sixthorder boundary value problems,” Computers & Mathematics with Applications, vol. 55, no. 12, pp. 2953–2972, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A.M. 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  MathSciNet
 C. H. Che Hussin and A. Kiliçman, “On the solutions of nonlinear higherorder boundary value problems by using differential transformation method and Adomian decomposition method,” Mathematical Problems in Engineering, vol. 2011, Article ID 724927, 19 pages, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. E. Simos, “New stable closed NewtonCotes trigonometrically fitted formulae for longtime integration,” Abstract and Applied Analysis, Article ID 182536, 15 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. E. Simos, “Optimizing a hybrid twostep method for the numerical solution of the Schrödinger equation and related problems with respect to phaselag,” Journal of Applied Mathematics, Article ID 420387, 17 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. A. Anastassi and T. E. Simos, “A parametric symmetric linear fourstep method for the efficient integration of the Schrödinger equation and related oscillatory problems,” Journal of Computational and Applied Mathematics, vol. 236, no. 16, pp. 3880–3889, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Akram and H. U. Rehman, “Solution of first order singularly perturbed initial value problem in reproducing kernel Hilbert space,” European Journal of Scientific Research, vol. 53, no. 4, pp. 516–523, 2011. View at: Google Scholar
 G. Akram and H. U. Rehman, “Numerical solution of eighth order boundary value problems in reproducing kernel space,” Numerical Algorithms, vol. 62, no. 3, pp. 527–540, 2013. View at: Publisher Site  Google Scholar
 G. Akram and H. U. Rehman, “Solution of fifth order boundary value problems in the reproducing kernel space,” Middle East Journal of Scientific Research, vol. 10, no. 2, pp. 191–195, 2011. View at: Google Scholar
 H. Yao, “New algorithm for the numerical solution of the integrodifferential equation with an integral boundary condition,” Journal of Mathematical Chemistry, vol. 47, no. 3, pp. 1054–1067, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. K. Pandey, “High order finite difference method for numerical solution of general two point boundary value problems involving sixth order differential equation,” International Journal of Pure and Applied Mathematics, vol. 76, no. 3, pp. 317–326, 2012. View at: Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2013 Ghazala Akram and Hamood Ur Rehman. 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.