- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Journal of Applied Mathematics
Volume 2014 (2014), Article ID 135465, 10 pages
A Numerical Iterative Method for Solving Systems of First-Order Periodic Boundary Value Problems
1Applied Science Department, Ajloun College, Al-Balqa Applied University, Ajloun 26816, Jordan
2Department of Mathematics, Al-Balqa Applied University, Salt 19117, Jordan
Received 12 September 2013; Accepted 11 February 2014; Published 25 March 2014
Academic Editor: Hak-Keung Lam
Copyright © 2014 Mohammed AL-Smadi 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.
The objective of this paper is to present a numerical iterative method for solving systems of first-order ordinary differential equations subject to periodic boundary conditions. This iterative technique is based on the use of the reproducing kernel Hilbert space method in which every function satisfies the periodic boundary conditions. The present method is accurate, needs less effort to achieve the results, and is especially developed for nonlinear case. Furthermore, the present method enables us to approximate the solutions and their derivatives at every point of the range of integration. Indeed, three numerical examples are provided to illustrate the effectiveness of the present method. Results obtained show that the numerical scheme is very effective and convenient for solving systems of first-order ordinary differential equations with periodic boundary conditions.
Systems of ordinary differential equations with periodic boundary value conditions, the so-called periodic boundary value problems (BVPs), are well known for their applications in sciences and engineering [1–5]. In this paper, we focus on finding approximate solutions to systems of first-order periodic BVPs, which are a combination of systems of first-order ordinary differential equations and periodic boundary conditions. In fact, accurate and fast numerical solutions of systems of first-order periodic BVPs are of great importance due to their wide applications in scientific and engineering research.
Numerical methods are becoming more and more important in mathematical and engineering applications, simply not only because of the difficulties encountered in finding exact analytical solutions but also because of the ease with which numerical techniques can be used in conjunction with modern high-speed digital computers. A numerical procedure for solving systems of first-order periodic BVPs based on the use of reproducing kernel Hilbert space (RKHS) method is discussed in this work.
Among a substantial number of works dealing with systems of first-order periodic BVPs, we mention [6–10]. The existence of solutions to systems of first-order periodic BVPs has been discussed as described in . In , the authors have discussed some existence and uniqueness results of periodic solutions for first-order periodic differential systems. Also, in  the authors have provided the existence, multiplicity, and nonexistence of positive periodic solutions for systems of first-order periodic BVPs. Furthermore, the existence of periodic solutions for the coupled first-order differential systems of Hamiltonian type is carried out in . Recently, the existence of positive solutions for systems of first-order periodic BVPs is proposed in . For more results on the solvability analysis of solutions for systems of first-order periodic BVPs, we refer the reader to [11–15], and for numerical solvability of different categories of BVPs, one can consult [16–19].
Investigation about systems of first-order periodic BVPs numerically is scarce. In this paper, we utilize a methodical way to solve these types of differential systems. In fact, we provide criteria for finding the approximate and exact solutions to the following system: subject to the periodic boundary conditions where , are unknown functions to be determined, are continuous terms in as , , in which , and , are two reproducing kernel spaces. Here, we assume that (1) subject to the periodic boundary conditions (2) has a unique solution on .
Reproducing kernel theory has important applications in numerical analysis, differential equations, integral equations, probability and statistics, and so forth [20–22]. In the last years, extensive work has been done using RKHS method which provides numerical approximations for linear and nonlinear equations. This method has been implemented in several operator, differential, integral, and integrodifferential equations side by side with their theories. The reader is kindly requested to go through [23–35] in order to know more details about RKHS method, including its history, its modification for use, its applications, and its characteristics.
The rest of the paper is organized as follows. In the next section, two reproducing kernel spaces are described in order to formulate the reproducing kernel functions. In Section 3, some essential results are introduced and a method for the existence of solutions for (1) and (2) is described. In Section 4, we give an iterative method to solve (1) and (2) numerically. Numerical examples are presented in Section 5. Section 6 ends this paper with brief conclusions.
2. Construct of Reproducing Kernel Functions
In this section, two reproducing kernels needed are constructed in order to solve (1) and (2) using RKHS method. Before the construction, we utilize the reproducing kernel concept. Throughout this paper, is the set of complex numbers, , and .
Definition 1 (see ). Let be a nonempty abstract set. A function is a reproducing kernel of the Hilbert space if(1)for each , ,(2)for each and , .
Remark 2. Condition in Definition 1 is called “the reproducing property,” which means that the value of the function at the point is reproducing by the inner product of with . A Hilbert space which possesses a reproducing kernel is called a RKHS.
To solve (1) and (2) using RKHS method, we first define and construct a reproducing kernel space in which every function satisfies the periodic boundary condition . After that, we utilize the reproducing kernel space .
Definition 3. The inner product space is defined as are absolutely continuous real-valued functions on , , and . On the other hand, the inner product and the norm in are defined, respectively, by and , where .
It is easy to see that satisfies all the requirements for the inner product. First, . Second, . Third, . Fourth, , where . It therefore remains only to prove that if and only if . In fact, it is obvious that when , then . On the other hand, if , then by (3), we have ; therefore, and . Then, we can obtain .
Definition 4 (see ). The Hilbert space is called a reproducing kernel if, for each fixed , there exist (simply ) such that for any and .
An important subset of the RKHSs is the RKHSs associated with continuous kernel functions. These spaces have wide applications, including complex analysis, harmonic analysis, quantum mechanics, statistics, and machine learning.
Theorem 5. The Hilbert space is a complete reproducing kernel and its reproducing kernel function can be written as where and , , are unknown coefficients of and will be given in the following proof.
Proof. The proof of the completeness and reproducing property of is similar to the proof in . Now, let us find out the expression form of the reproducing kernel function in the space . Through several integration by parts, we have . Thus, from (3), we can write + + . Since , it follows that ; also since , it follows that . Then But on the other aspect as well, if , , , and , then (5) implies that . Now, for any , if satisfies then . Obviously, is the reproducing kernel function of the space . Next, we give the expression form of the reproducing kernel function . The characteristic formula of (6) is given by . Then the characteristic values are with multiplicity . So, let the expression form of the reproducing kernel function be as defined in (4). On the other hand, for (6), let satisfy the equation , . Integrating from to with respect to and letting , we have the jump degree of at given by . Through the last descriptions, the unknown coefficients of (4) can be obtained. However, by using MAPLE software package, the representation form of the reproducing kernel function is provided byThis completes the proof.
Definition 6 (see ). The inner product space is defined as is absolutely continuous real-valued function on and . On the other hand, the inner product and the norm in are defined, respectively, by , and , where .
Theorem 7 (see ). The Hilbert space is a complete reproducing kernel and its reproducing kernel function can be written as
Reproducing kernel functions possess some important properties such as being symmetric, unique, and nonnegative. The reader is asked to refer to [23–35] in order to know more details about reproducing kernel functions, including their mathematical and geometrical properties, their types and kinds, and their applications and method of calculations.
3. Formulation of Linear Operator
In this section, the formulation of a differential linear operator and the implementation method are presented in the reproducing kernel space . After that, we construct an orthogonal function system of the space based on the use of the Gram-Schmidt orthogonalization process in order to obtain the exact and approximate solutions of (1) and (2) using RKHS method.
First, as in [23–35], we transform the problem into a differential operator. To do this, we define a differential operator as such that . As a result, (1) and (2) can be converted into the equivalent form as follows: where and in which and for , , and . It is easy to show that is a bounded linear operator from the space into the space .
Initially, we construct an orthogonal function system of . To do so, put and , where is dense on and is the adjoint operator of . In terms of the properties of reproducing kernel function , one obtains , , .
For the orthonormal function system of the space , it can be derived from the Gram-Schmidt orthogonalization process of as follows: where are orthogonalization coefficients and are given as
Clearly, . Thus, can be written in the form , where indicates that the operator applies to the function of .
Theorem 8. If is dense on , then is a complete function system of the space .
Proof. For each fixed , let . In other words, one can write . Note that is dense on ; therefore . It follows that , , from the existence of . So, the proof of the theorem is complete.
Lemma 9. If , then there exist positive constants such that , , , where .
Proof. For any , we have . By the expression form of the kernel function , it follows that . Thus, . Hence, , , .
The internal structure of the following theorem is as follows: firstly, we will give the representation form of the exact solutions of (1) and (2) in the form of an infinite series in the space . After that, the convergence of approximate solutions to the exact solutions , , will be proved.
Theorem 10. For each , in the space , the series is convergent in the sense of the norm of . On the other hand, if is dense on , then the following hold:(i)the exact solutions of (9) could be represented by (ii)the approximate solutions of (9) and , , are converging uniformly to the exact solutions and their derivatives as , respectively.
Proof. For the first part, let be solutions of (9) in the space . Since , is the Fourier series expansion about normal orthogonal system , and is the Hilbert space, then the series is convergent in the sense of . On the other hand, using (10), it easy to see that Therefore, the form of (12) is the exact solutions of (9). For the second part, it is easy to see that by Lemma 9, for any , where and are positive constants. Hence, if as , the approximate solutions and , , , are converged uniformly to the exact solutions and their derivatives, respectively. So, the proof of the theorem is complete.
4. Construction of Iterative Method
In this section, an iterative method of obtaining the solutions of (1) and (2) is represented in the reproducing kernel space for linear and nonlinear cases. Initially, we will mention the following remark about the exact and approximate solutions of (1) and (2).
Case 2. If (1) is nonlinear, then in this case the exact and approximate solutions can be obtained by using the following iterative algorithm.
Algorithm 11. According to (12), the representation form of the solutions of (1) can be denoted by
where , . In fact, in (16) are unknown; one will approximate using known . For numerical computations, one defines the initial functions , put , and define the -term approximations to by
where the coefficients of , , , are given as
Here, we note that, in the iterative process of (17), we can guarantee that the approximations satisfy the periodic boundary conditions (2). Now, the approximate solutions can be obtained by taking finitely many terms in the series representation of and
Now, we will proof that in the iterative formula (17) are converged to the exact solutions of (1). In fact, this result is a fundamental in the RKHS theory and its applications. The next two lemmas are collected in order to prove the prerecent theorem.
Lemma 12. If , as , and is continuous in with respect to , for and , then , as .
Proof. Firstly, we will prove that in the sense of . Since By reproducing property of , we have and . Thus, . From the symmetry of , it follows that as . Hence, as soon as . On the other hand, by Theorem 10 part (ii), for any , it holds that as . Therefore, in the sense of as and . Thus, by means of the continuation of , it is obtained that , as .
Lemma 13. For , one has , .
Proof. The proof of will be obtained by induction as follows: if , then . Using the orthogonality of , it yields that Now, if , then . Again, if , then ,. Thus, . Indeed, it is easy to see by using mathematical induction that , . But on the other hand, from Theorem 10, converge uniformly to . It follows that, on taking limits in (17), . Therefore, , where is an orthogonal projector from the space to Span . Thus, as and .
Proof. The proof consists of the following three steps. Firstly, we will prove that the sequence in (17) is monotone increasing in the sense of . By Theorem 8, is the complete orthonormal system in the space . Hence, we have . Therefore, , , is monotone increasing. Secondly, we will prove the convergence of . From (17), we have . From the orthogonality of , it follows that . Since, the sequence is monotone increasing in the sense of . Due to the condition that is bounded, is convergent as . Then, there exist constants such that . It implies that , . On the other hand, since it follows for that Furthermore, . Consequently, as , we have . Considering the completeness of , there exists such that , as in the sense of . Thirdly, we will prove that are the solutions of (9). Since is dense on , for any , there exists subsequence , such that as . From Lemma 13, it is clear that . Hence, let ; by Lemma 12 and the continuity of , we have . That is, satisfies (1). Also, since , clearly, satisfies the periodic boundary conditions (2). In other words, are the solutions of (1) and (2), where and are given by (18). The proof is complete.
According to the internal structure of the present method, it is obvious that if we let denote the exact solutions of (9), denote the approximate solutions obtained by the RKHS method as given by (17), and denote the difference between and , where and , then and or . Consequently, this shows the following theorem.
Theorem 15. The difference , , is monotone decreasing in the sense of the norm of .
5. Numerical Examples
In this section, the theoretical results of the previous sections are illustrated by means of some numerical examples in order to illustrate the performance of the RKHS method for solving systems of first-order periodic BVPs and justify the accuracy and efficiency of the method. To do so, we consider the following three nonlinear examples. These examples have been solved by the presented method with different values of and . Results obtained by the method are compared with the exact solution of each example by computing the absolute and relative errors and are found to be in good agreement with each other. In the process of computation, all experiments were performed in MAPLE software package.
Example 1. Consider the following first-order nonlinear differential system: subject to the periodic boundary conditions The exact solutions are and .
The present method enables us to approximate the solutions and their derivatives at every point of the range of integration. Hence, it is possible to pick any point in and as well the approximate solutions and their derivatives will be applicable. Next, new numerical results for Example 1 which include the absolute error at some selected gird points in for approximating and , where , , , and , are given in Table 3.
Example 2. Consider the following first-order nonlinear differential system: subject to the periodic boundary conditions The exact solutions are and .
Example 3. Consider the following first-order nonlinear differential system: subject to the periodic boundary conditions The exact solutions are , , and .
From the previous tables, it can be seen that the RKHS method provides us with the accurate approximate solutions. On the other aspect as well, it is clear that the accuracy obtained using the mentioned method is advanced by using only a few tens of iterations.
Here, we use the RKHS method to solve systems of first-order periodic BVPs. The solutions were calculated in the form of a convergent series in the space with easily computable components. In the proposed method, the -term approximations are obtained and proved to converge to the exact solutions. Meanwhile, the error of the approximate solutions is monotone decreasing in the sense of the norm of . It is worthy to note that, in our work, the approximate solutions and their derivatives converge uniformly to the exact solutions and their derivatives, respectively. On the other aspect as well, the present method enables us to approximate the solutions and their derivatives at every point of the range of integration. The results show that the present method is an accurate and reliable analytical technique for solving systems of first-order periodic BVPs.
Conflict of Interests
The authors declare that there is no conflict of interests.
The authors would like to express their thanks to unknown referees for their careful reading and helpful comments.
- E. Coddington and N. Levinson, Theory of Ordinary Differential Equations, McGraw-Hill, New York, NY, USA, 1955.
- H. I. Freedman and J. H. Wu, “Periodic solutions of single-species models with periodic delay,” SIAM Journal on Mathematical Analysis, vol. 23, no. 3, pp. 689–701, 1992.
- J. Mawhin and J. R. Ward, “Nonuniform nonresonance conditions at the two first eigenvalues for periodic solutions of forced Lienard and Duffing equations,” Rocky Mountain Journal of Mathematics, vol. 12, no. 4, pp. 643–654, 1982.
- A. C. Lazer, “Application of a lemma on bilinear forms to a problem in nonlinear oscillations,” Proceedings of the American Mathematical Society, vol. 33, pp. 89–94, 1972.
- K. Abd-Ellateef, R. Ahmed, and Z. Drici, “Generalized quasilinearization for systems of nonlinear differential equations with periodic boundary conditions,” Dynamics of Continuous, Discrete & Impulsive Systems A, vol. 12, no. 1, pp. 77–85, 2005.
- C. C. Tisdell, “Existence of solutions to first-order periodic boundary value problems,” Journal of Mathematical Analysis and Applications, vol. 323, no. 2, pp. 1325–1332, 2006.
- R. P. Agarwal and J. Chen, “Periodic solutions for first order differential systems,” Applied Mathematics Letters, vol. 23, no. 3, pp. 337–341, 2010.
- R. Chen, R. Ma, and Z. He, “Positive periodic solutions of first-order singular systems,” Applied Mathematics and Computation, vol. 218, no. 23, pp. 11421–11428, 2012.
- C. P. Gupta, “Periodic solutions for coupled first order nonlinear differential systems of Hamiltonian type,” Nonlinear Analysis: Theory, Methods & Applications, vol. 8, no. 11, pp. 1271–1285, 1984.
- Q. Kong and M. Wang, “Positive solutions of even order system periodic boundary value problems,” Nonlinear Analysis: Theory, Methods & Applications, vol. 72, no. 3-4, pp. 1778–1791, 2010.
- A. Boucherif and N. Merabet, “Boundary value problems for first order multivalued differential systems,” Archivum Mathematicum, vol. 41, no. 2, pp. 187–195, 2005.
- J. M. Ortega and W. C. Rheinboldt, Iterative Solution of Nonlinear Equations in Several Variables, Academic Press, New York, NY, USA, 1970.
- W. G. Li, “Solving the periodic boundary value problem with the initial value problem method,” Journal of Mathematical Analysis and Applications, vol. 226, no. 1, pp. 259–270, 1998.
- J. Mawhin, Topological Degree Methods in Nonlinear Boundary Value Problems, vol. 40 of CBMS Regional Conference Series in Mathematics, American Mathematical Society, Providence, RI, USA, 1979.
- I. T. Kiguradze, “On periodic solutions of th order ordinary differential equations,” Nonlinear Analysis: Theory, Methods & Applications, vol. 40, no. 1–8, pp. 309–321, 2000.
- O. abu Arqub, A. el-Ajou, Z. al Zhour, and S. Momani, “Multiple solutions of nonlinear boundary value problems of fractional order: a new analytic iterative technique,” Entropy, vol. 16, no. 1, pp. 471–493, 2014.
- Z. abo-Hammour, O. abu Arqub, S. Momani, and N. Shawagfeh, “Optimization solution of Troesch's and Bratu's problems of ordinary type using novel continuous genetic algorithm,” Discrete Dynamics in Nature and Society, vol. 2014, Article ID 401696, 15 pages, 2014.
- O. abu Arqub, Z. abo-Hammour, and S. Momani, “Application of continuous genetic algorithm for nonlinear system of second-order boundary value problems,” Applied Mathematics and Information Sciences, vol. 8, pp. 235–248, 2014.
- O. abu Arqub, Z. abo-Hammour, S. Momani, and N. Shawagfeh, “Solving singular two-point boundary value problems using continuous genetic algorithm,” Abstract and Applied Analysis, vol. 2012, Article ID 205391, 25 pages, 2012.
- A. Berlinet and C. Thomas-Agnan, Reproducing Kernel Hilbert Spaces in Probability and Statistics, Kluwer Academic, Boston, Mass, USA, 2004.
- M. Cui and Y. Lin, Nonlinear Numerical Analysis in the Reproducing Kernel Space, Nova Science, New York, NY, USA, 2008.
- A. Daniel, Reproducing Kernel Spaces and Applications, Springer, Basel, Switzerland, 2003.
- F. Geng, “Solving singular second order three-point boundary value problems using reproducing kernel Hilbert space method,” Applied Mathematics and Computation, vol. 215, no. 6, pp. 2095–2102, 2009.
- L.-H. Yang and Y. Lin, “Reproducing kernel methods for solving linear initial-boundary-value problems,” Electronic Journal of Differential Equations, vol. 2008, pp. 1–11, 2008.
- Y. Z. Lin, M. G. Cui, and L. H. Yang, “Representation of the exact solution for a kind of nonlinear partial differential equation,” Applied Mathematics Letters, vol. 19, no. 8, pp. 808–813, 2006.
- W. Wang, M. Cui, and B. Han, “A new method for solving a class of singular two-point boundary value problems,” Applied Mathematics and Computation, vol. 206, no. 2, pp. 721–727, 2008.
- W. Jiang and Z. Chen, “Solving a system of linear Volterra integral equations using the new reproducing kernel method,” Applied Mathematics and Computation, vol. 219, no. 20, pp. 10225–10230, 2013.
- F. Geng and M. Cui, “A reproducing kernel method for solving nonlocal fractional boundary value problems,” Applied Mathematics Letters, vol. 25, no. 5, pp. 818–823, 2012.
- F. Z. Geng and S. P. Qian, “Reproducing kernel method for singularly perturbed turning point problems having twin boundary layers,” Applied Mathematics Letters, vol. 26, no. 10, pp. 998–1004, 2013.
- W. Jiang and Z. Chen, “A collocation method based on reproducing kernel for a modified anomalous subdiffusion equation,” Numerical Methods for Partial Differential Equations, vol. 30, no. 1, pp. 289–300, 2014.
- F. Z. Geng, S. P. Qian, and S. Li, “A numerical method for singularly perturbed turning point problems with an interior layer,” Journal of Computational and Applied Mathematics, vol. 255, pp. 97–105, 2014.
- N. Shawagfeh, O. abu Arqub, and S. Momani, “Analytical solution of nonlinear second-order periodic boundary value problem using reproducing kernel method,” Journal of Computational Analysis and Applications, vol. 16, pp. 750–762, 2014.
- M. al-Smadi, O. abu Arqub, and S. Momani, “A computational method for two-point boundary value problems of fourth-order mixed integrodifferential equations,” Mathematical Problems in Engineering, vol. 2013, Article ID 832074, 10 pages, 2013.
- O. abu Arqub, M. al-Smadi, and N. Shawagfeh, “Solving Fredholm integro-differential equations using reproducing kernel Hilbert space method,” Applied Mathematics and Computation, vol. 219, no. 17, pp. 8938–8948, 2013.
- O. abu Arqub, M. al-Smadi, and S. Momani, “Application of reproducing kernel method for solving nonlinear Fredholm-Volterra integrodifferential equations,” Abstract and Applied Analysis, vol. 2012, Article ID 839836, 16 pages, 2012.