On the solutions of three-point boundary value problems using variational-fixed point iteration method

Given a three-point fourth-order boundary value problems y(iv)+p(x)y″′+q(x)y″+r(x)y′+s(x)y=f(x),a≤x≤b\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} y^{(iv)}+p(x)y^{\prime \prime \prime }+q(x)y^{\prime \prime }+r(x)y^{\prime }+s(x)y=f(x),a \le x \le b \end{aligned}$$\end{document}such that y(a)=y(b)=y″(b)=y″(α)=0,a≤α≤b;\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} y(a)=y(b)=y^{\prime \prime }(b)=y^{\prime \prime }(\alpha )=0,a \le \alpha \le b; \end{aligned}$$\end{document}where p,q,r,s,f∈C[a,b]\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p,q,r,s,f \in C[a,b]$$\end{document}, we combine the application of variational iteration method and fixed point iteration process to construct an iterative scheme called variational-fixed point iteration method that approximates the solution of three-point boundary value problems. The success of the variational or weighted residual method of approximation from a practical point of view depends on the suitable selection of the basis function. The method is self correcting one and leads to fast convergence. Problems were experimented to show the effectiveness and accuracy of the proposed method.


Introduction
The numerical solution of boundary value problems is of great importance as a result of its wide application in scientific and technological research [1]. Many researchers have developed various numerical methods, especially iterative methods to approximate different types of differential equations, see [2][3][4][5]. In recent years, there has been a growing interest in the treatment of iterative approximation such as variational iteration method, fixed point iteration and so on [6], variational iteration method has been used over the years to obtain an approximate solution of some boundary value problems, see [7,8]. On the other hand, fixed point iteration is a method of computing fixed point of iterated function, it is a wellknown method of approximation whose version is the variation iteration method, see [9,10]. Obviously, these methods have been proved by many researchers to be powerful tools in solving boundary value problems [11]. However, there are noticeable shortcomings in implementations of these methods especially the use of arbitrary function as a starting value, an in-appropriate choice of starting function may affect the rate of convergence, see [11,12].
In this paper, the propose method is an elegant combination of variational iteration method and fixed point iteration method with the use of finite element method to determine the starting function.
To illustrate the basic technique of variational iteration method, we consider the following general differential equation where L is linear operator, N is nonlinear operator and g(x) is forcing term. According to variational iteration method, see [13,14], a correctional functional of (1) can be constructed as follows where k is the Lagrange multiplier [15,16] which can be determined by variational theory, thus kðLu n ðsÞ þ Nu n ðsÞ À gðsÞÞds its stationary conditions can be obtained using integration by parts. The second term in the right is called correction andũu n is considered as a restricted variation, i.e du n ¼ 0.
Thus, the Lagrange's multiplier [17,19] can be identified as kðxÞ ¼ ðÀ1Þ n ðs À xÞ nÀ1 =ðn À 1Þ!: The given solution is considered as the fixed point of the following functional [11] under the suitable choice of initial approximation [21] u 0 ðxÞ at n ¼ 0. We use finite element method to determine the starting function to avoid the arbitrary choice of starting function. After single iteration process we obtained u 1 ; repeating the process iteratively the term to be integrated become larger and cumbersome to operate and the result of each iteration step diverges from exact solution. Based on this fact we introduce the application of fixed point iterative process to overcome it. Thus we have the following theorem: Theorem 1 Let (E, d) be a complete metric space and T be a self map on E. Further, let y o 2 E and let y nþ1 ¼ f ðT; y n Þ denotes an iteration procedure which gives a sequence fy n g . Then T is an iteration process and defined for arbitrary y 0 by y nþ1 ¼ f ðT; y n Þ ¼ ð1 À a n Þy n þ a n Ty n ; n ! 0 ð3Þ where fa n g is a real sequence satisfying a n ¼ 1; 0 a n 1 for n ! 0 and P 1 n¼0 a n ¼ 1: Proof Let y ðivÞ þ pðxÞy 000 þ qðxÞy 00 þ rðxÞy 0 þ sðxÞy ¼ tðxÞ such that yðaÞ ¼ yðbÞ ¼ y 00 ðbÞ ¼ y 00 ðaÞ ¼ 0; a a b; where p; q; r; s; f 2 C½a; b, then the scheme y ðivÞ nþ1 ¼ ð1 À k n Þy ðivÞ n þ k n ð1 À aÞy ðivÞ n ð4Þ obtained by harnessing Mann and Banach fixed point iteration [6,20], to yield y ðivÞ nþ1 ¼ kðtðxÞ À pðxÞy 000 n À qðxÞy 00 n À rðxÞy 0 n À sðxÞÞy n þ ð1 À kÞy ðivÞ then T is an operator such that any y(x) is a solution of (6) at fixed point of T and can be referred as a fixed point operator. For convergence of (3) and (5) we let y nþ1 ¼ ð1 À a n Þy n þ a n Ty n y 0 nþ1 ¼ ð1 À a n Þy 0 n þ a n Ty 0 n y 00 nþ1 ¼ ð1 À a n Þy 00 n þ a n Ty 00 n y 000 nþ1 ¼ ð1 À a n Þy 000 n þ a n Ty 000 n y ðivÞ nþ1 ¼ ð1 À a n Þy ðivÞ n þ a n Ty ðivÞ n : Therefore Eqs. (9) and 10 become y ðivÞ nþ1 ¼ ð1 À a n Þy ðivÞ n þ a n Z x 2 x 1 d 4 =dx 4 ðGðx; sÞf ðt; yðtÞ; y 0 ðtÞ; Â y 00 ðtÞ; y 000 ðtÞÞdt þ vðxÞ ðivÞ : Also, combining Eqs. (7) and (11) yields y ðivÞ nþ1 ¼ ð1 À a n Þy ðivÞ n þ a n Ty ðivÞ n : ð12Þ Therefore scheme (3) and (5) are convergent. This scheme is use to approximate boundary value problems iteratively with the use of arbitrary initial approximation y 0 at n ¼ 0. Since the variational iteration method and fixed point iteration methods are similar [10]. We let y 0 ¼ u 1 to avoid the assumption of the arbitrary function y 0 where the process will be carried out iteratively until convergence is obtained or the iteration is terminated. However, the finite element methods are the Galerkin method, collocation method, Raleigh-Ritz method, etc. Galerkin method is an approximate solution of boundary value problems suggested by Galerkin [18], based on the requirement that the basis function / 0 ; / 1 ; / 2 ; . . .; / n be orthogonal to the residual Z wðx i ; a 0 ; a 1 ; . . .; a n Þ/ i dx ¼ 0 . . .; n. This gives rise to the following system of linear algebraic equations for the coefficients of the approximation solution y n ðxÞ ¼ / 0 ðxÞ þ a 1 / 1 ðxÞ þ a 2 / 2 ðxÞ þ Á Á Á þ a n / n ðxÞ of the boundary value problem Ly ¼ y ðivÞ þ pðxÞy 000 þ qðxÞy 00 þ rðxÞy 0 þ sðxÞy ¼ f ðxÞ; 0 x 1 such that Therefore, . . a 1 Lð/ 1 ; / n Þ þ a 2 L/ 2 ; / n þ Á Á Á þ a n L/ n ; / n ¼ ðf À Lð/ 0 ; / n ÞÞ: The weight functions are taken with the concept of inner product and orthogonality, it is obvious that the inner product of the two function in a certain domain is which is used to determine the starting function of the variational iteration method instead of an arbitrary choice. h

Numerical examples
In this section, two experiments are considered to demonstrate present methods: subject to the boundary conditions yð0Þ ¼ y 00 ð1=2Þ ¼ y 000 ð0Þ ¼ 0; yð1Þ ¼ 1: The following must be observed: Galerkin method of approximation is used to determine the initial approximation of variational iteration method whose trial function is called the basis function, the approximate solution we sought. Where U 0 ðxÞ ¼ x; We differentiate the Eq. (14) successively to obtain second and fourth derivatives and then substitute in Eq. (13) to get the residual A weight function is chosen within the bases function with the concepts of inner product and orthogonality These are sets of n-order linear equations which is solved to obtain all C i coefficients as follows: we obtain Also solving these Eqs. (16) and (17) simultaneously we obtain C ¼ 308=4331; C 2 ¼ 42=4331. We substitute these constants in (14), hence the approximation solution we sought for.
The second step is the use variational iteration method to determine the starting function of the fixed point iterative procedure; we construct the correct functional of (13) as follows: t nþ1 ðxÞ ¼ t n ðxÞ þ Z x 0 kðt ð4Þ n ðsÞ À t ð2Þ n ðsÞ À gðsÞÞds at n ¼ 0 We let we differentiate (20), i.e t 0 ðxÞ successively to obtain its second and fourth derivative and substitute same in Eq. (19) to have When the process is repeated for further iterations, the function to be integrated is getting larger and complex, where iterated values diverge from the analytical solution. Based on this fact, We let y 0 ¼ U 1 ðxÞ, where U 1 ðxÞ is the iterative function obtained after single iteration taken as an initial values for fixed point iterative technique. The scheme in (12) the fixed point iterative procedure can be used as follows y ðivÞ nþ1 ðxÞ ¼ y 00 n ðxÞ À 2 ð22Þ at n ¼ 0 y ðivÞ 1 ðxÞ ¼ y 00 0 ðxÞ À 2: But We differentiate the Eq. (24) twice to get y 00 0 and then substituted it into the given Eq. (23) to obtain y ðivÞ 1 ðxÞ ¼ À1:29324819x 2 þ 1:939505888x 3 À 3:057954283x 4 À 0:4040637266x 5 À 2: To obtain y 1 ðxÞ, we integrate (25) four times successively and imposing the boundary conditions to get The repeat the process at n ¼ 2; 3. . . until the iteration is terminated or it converges with the analytical solution, after few iterations we get y 5 ðxÞ ¼ 0:9632277022x þ 0:1131826518x 2 À 0:0739028367x 4 À 0:002462970407x 6 À 0:00004406199500x 8 À 4:8098132 Â 10 À7 x 10 À 4:175351396 Â 10 À9 x 12 À 1:223616824 Â 10 À12 x 16 À 1:145095958 Â 10 À13 x 17 þ 3:559612579 Â 10 À12 x 15 À 2:294149120 Â 10 À11 x 14 : Example 2 Consider y ðivÞ À y 00 ¼ À12x subject to the boundary conditions yð0Þ ¼ y 00 ð0Þ ¼ y 000 ð1=2Þ ¼ 0; yð1Þ ¼ 13: We let which is the basis function where U 0 ðxÞ ¼ 13x; We differentiate Eq. (27) successively to obtain second and fourth derivatives and then substitute same in (26) the given differential equation to get the residual A weight function are chosen within the basis function with the concepts of inner product and orthogonality These are sets of n-order linear equations which is solved to obtain all C i coefficients as follows: Also solving Eqs. (29) and (30) simultaneously to obtain C 1 ¼ À635=4331; C 2 ¼ 504=4331. Then we substitute these constants in (27) called the basis function and is the approximate solution we want Next we apply variational iteration method by constructing a correctional functional of (26) t nþ1 ðxÞ ¼ t n ðxÞ þ Z x 0 kðt ð4Þ n ðsÞ À t ð2Þ n ðsÞ À gðsÞÞds at n ¼ 0 We let we differentiate (33), i.e t 0 ðxÞ successively to obtain it second and fourth derivative and substitute in the given Eq. (32) to have we simplify to get t 1 ðxÞ ¼ 12:96975294x À 0:2932348187x 2 þ 1:163703533x 3 þ 0:7597475208x 4 À 6:130182374x 5 þ 4:267898383x 6 þ 0:9699769053x 7 : As the number of iterations increases, the function to be integrated is getting larger and cumbersome making the functions diverging from the analytical solutions at each iteration steps. Based on this fact, we apply the fixed point iterative technique where we use U 1 ðxÞ ¼ y 0 as a starting value after single iteration to avoid the use of arbitrary function (Fig. 1).

Discussion
The accuracy and convergence of the method are of great significance in a numerical experiment of such type. Accuracy measures the degree of closeness of the numerical solution to the theoretical solution while convergence measures even the closer approach of successive iteration to the exact solution as the number of iteration increases. To asses the success of our method, the scheme was tested with some numerical examples whose results presented as Tables 1 and 2. These tables show the comparison with exact method, Galerkin method, variational iteration method and the approximate solution obtained by variational-fixed point iteration method. It is observed that with few iterations, the order of the error is quite encouraging, which indicates fast rate of convergence. It is clearly seen that as the iterations proceeds, the error decreases and convergence is assured (Fig. 2).

Conclusion
In this paper we have shown the performance of the variational-fixed point iterative scheme for the solution threepoint boundary value problems with help of some experiments which indicates to be very powerful and an efficient technique with good convergence property when compared with the existing method. We can conclude the method have an advantage over the existing methods.

Compliance with ethical standards
Conflict of interest The authors declare that they have no competing interests.
Author's contributions Both authors jointly worked on deriving the results and approved the final manuscript.