In this paper, the following NTBVPs are considered
where f ∈ (C[a, b] × R2, R) is a nonlinear function, la, ma, lb, and mb are given positive constants and ka, kb are given constants.
We consider the above NTBVPs which occur in applied mathematics, theoretical physics, engineering, control and optimization theory. By various researchers, different numerical methods have been proposed. Among these are the Adomian decomposition method (ADM) , He-Laplace method , Homotopy perturbation method (HPM) . There are many methods for solving of two-point boundary value problems ODEs.
Additionally, many papers stated the BCM. In Ref. , authors utilized the BCM for solving delay differential equation (DDE) and optimal control of switched systems numerically. In Ref. , authors stated the utilization of BCM on some linear optimal control systems with pantograph delays. For solving the quadratic Riccati differential equation BCM is utilized (see Ref. ). In Ref. , authors presented BCM for solving singular nonlinear equations of Lane-Emden type. In this work, one may utilize BCM and orthonormal Bernstein polynomial (OBP) method for solving NTBVPs.
The Bernstein polynomials are not orthonormal therefore Gram-Schmidt technique could be utilized to construct OBP.
The organization of this paper is classified as follows: Problem statement is introduced in Section 2. Orthonormal Bernstein polynomials of arbitrary order are presented in Section 3. Numerical examples are solved in Section 4. Finally, Section 5 will give a conclusion briefly.
2 Problem statement
Our strategy is utilizing BCM to approximate the solutions y(x) where y(x) is given below. Define the Bezier polynomials of degree n over the interval [x0, xf] as follows:
where h = xf – x0, and
is the Bernstein polynomial of degree n over the interval [x0, xf], ar represents the control points. Suppose that MB is the coefficient matrix of Br,n(x), r = 0, 1, …, n, where MB(i, j) is the coefficient of the Bi,n(x) with respect to the monomial xj–1, then by Eq. (3) one may have
By substituting y(x) in Eq. (1), one may define R1(x, a0, a1, …, an) for x ∈ [x0, xf] as follows:
The convergence of this method is proven where n → ∞ (see ).
Now, we define the residual function over the interval [x0, xf] as follows
3 Orthonormal Bernstein polynomials of arbitrary order
In this section, OBP is constructed. For arbitrary n, Eq. (6) can be obtained by induction. Now, let n = 3, for simplicity Bi = Bi(x) (i = 0, 1, 2, 3) and OBi = OBi(x), hence
where OBi (i ≥ 0) can be obtained from Gram-Schmidt orthonormalization technique
then the following OBP is considered as
4 Numerical examples
Now, some numerical examples of NTBVPs are stated to illustrate BCM. All results are obtained by utilizing Maple 14.
First, the following NTBVP is considered (see )
Applying the proposed technique with n = 4, one may have
where the obtained error is zero. The graphs of approximated and exact solution y are plotted in Figure 1.
Third, the following NTBVP is considered (see )
Applying the proposed technique with n = 4, one may have
This paper deals with the approximate solution of NTBVPs via BCM and OBP. The unity of the function of orthogonality for Bernstein polynomials (BPs) and the simplicity of applying BPs are great merits that make the technique very attractive and easy to use. Illustrative examples are solved to express the validity and applicability of the proposed methods. The approximate solution is coming very close to the exact solution.
The authors would like to thank the anonymous reviewers of this paper for their careful reading, constructive comments and nice suggestions, which have improved the paper very much. The Second author wishes to acknowledge financial support from the University of Venda.
Wazwaz A.M., Adomian decomposition method for a reliable 5 treatment of the Bratu type equations, Appl. Math. Comp., 2005, 166, 652-663. Google Scholar
Mishra H.K., He-Laplace method for the solution of two-point boundary value problems, Amer. J. Math. Anal., 2014, 2(3), 45-49. Google Scholar
Xu L., He’s homotopy perturbation method for a boundary layer equation in unbounded domain, Comp. Math. Appl., 2007, 54, 1067-1070. Google Scholar
Cuomo S., Marasco A., A numerical approach to nonlinear two-point boundary value problems for ODEs, Comp.Math. Appl., 2008, 55, 2476-2489. Google Scholar
Ghomanjani F., Farahi M.H., Optimal control of switched systems based on bezier control points, Int. J. Intell. Syst. Appl., 2012, 4(7), 16-22. Google Scholar
Ghomanjani F., Farahi M.H., Kamyad A.V., Numerical solution of some linear optimal control systems with pantograph delays, IMA J. Math. Control Inf., 2015, 32(1), 225-243. Google Scholar
Ghomanjani F., Khorram E., Approximate solution for quadratic Riccati differential equation, J. Taibah Univ. Sci., 2017, 11 246-250. Google Scholar
Adiyaman M., Oger V., A residual method using Bezier curves for singular nonlinear equations of Lane-Emden type, Kuwait J. Sci., 2017, 44(4), 9-18. Google Scholar
Ghomanjani F., Farahi M.H., Gachpazan M., Bezier control points method to solve constrained quadratic optimal control of time varying linear systems, Comp. Appl. Math., 2012, 31(3), 1-24. Google Scholar
About the article
Published Online: 2018-06-26
Conflict of InterestConflict of Interests: The author declares that there is no conflict of interest in the use of the above mentioned software.
Citation Information: Open Physics, Volume 16, Issue 1, Pages 371–374, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0050.
© 2018 Fateme Ghomanjani and Stanford Shateyi, published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0