- Research
- Open access
- Published:
Existence and approximation of solution for a nonlinear second-order three-point boundary value problem
Boundary Value Problems volume 2016, Article number: 169 (2016)
Abstract
A nonlinear second-order ordinary differential equation with four cases of three-point boundary value conditions is studied by investigating the existence and approximation of solutions. First, the integration method is proposed to transform the considered boundary value problems into Hammerstein integral equations. Second, the existence of solutions for the obtained Hammerstein integral equations is analyzed by using the Schauder fixed point theorem. The contraction mapping theorem in Banach spaces is further used to address the uniqueness of solutions. Third, the approximate solution of Hammerstein integral equations is constructed by using a new numerical method, and its convergence and error estimate are analyzed. Finally, some numerical examples are addressed to verify the given theorems and methods.
1 Introduction
Nonlocal boundary value problems for linear and nonlinear ordinary differential equations are arising in the theory of mathematical physics and some engineering applications [1–3]. They have attracted much attention and lots of interesting observations have been given [4–7]. The existence and approximation of solutions are very important in order to understand various phenomena in physics, engineering, and so on [8]. Here we generally consider the following nonlinear second-order ordinary differential equation:
where \(\psi(x,\varphi(x))\) and \(g(x)\) are known, and \(\varphi(x)\) is the unknown to be computed. The nonlinear ordinary differential equation (1) can be understood as the case that a linear problem is imposed upon a nonlinear loading \(\psi(x,\varphi(x))\). Moreover, the three-point boundary value condition is equipped and it is assumed to be one of the following four cases:
where α, β, ξ, and k are the known constants. It is found that II and IV can be reformulated by using the variable transformation \(u(x)=a+b-x\) from I and III, respectively. That is, under the assumption of \(\bar{\varphi}(u)=\varphi(a+b-u)=\varphi(x)\), one has
by using case I, and
from case III. Since the constants α, β, ξ, and k are arbitrary, the cases of I and III will be considered mainly in the following sections. When \(k=0\), the three-point boundary value problems degenerate to two-point boundary value problems.
It is noted that when \(g(x)=0\), the existence of solutions for equation (1) with various boundary value conditions has been studied widely. For example, the method of lower and upper solutions is developed by Ma [9] and the multiplicity solutions for a three-point boundary value problem at resonance were given. Xu [10, 11] considered the singular three-point and m-point boundary value problems, respectively. The multiplicity results and existence of positive solutions were analyzed by using a fixed point index theory. Yao [12] investigated the existence of positive solutions for a second-order three-point boundary value problem and a successive iteration method was given for computing the solutions. Nieto [13] studied the existence of solution for a second-order nonlinear ordinary differential equation with three-point boundary value conditions at resonance. In addition, the problems of nonlinear second-order ordinary differential equations with m-point and integral boundary value conditions were further investigated in [14–19] and so on. On the other hand, when \(\psi(x,\varphi(x))=a(x)F(x,\varphi(x))\) and \(\psi (x,\varphi(x))\) degenerates to \(a(x)\psi(\varphi(x))\), the multipoint boundary value problems of nonlinear second-order ordinary differential equations were dealt with in [20–27]. Recently, the special linear case of \(\psi(x,\varphi(x))=q(x)\varphi(x)\) has been studied in [28] and an approximate solution has been given. The case of \(g(x)=0\) in equation (1) with impulse three-point boundary value conditions have been studied in [29], where the existence conditions for obtaining a nontrivial solution have been given.
As shown in the above-mentioned work, the existence of the solutions for second-order multipoint boundary value problems is always focused on. Moreover, the approximate solutions of boundary value problems are very important in engineering applications. A monotone iterative technique was developed for the approximate solution of a second-order three-point boundary value problem in [30]. This paper generally focuses on the nonlinear second-order ordinary differential equation with four cases of three-point boundary value conditions in (1)-(5). A general method is proposed to transform the nonlinear three-point boundary value problems into nonlinear Hammerstein integral equations. The existence and uniqueness of solutions for the obtained Hammerstein integral equations are considered by using the Schauder fixed point theorem and the contraction mapping theorem, respectively. A new numerical method is further proposed to construct the approximate solutions of Hammerstein integral equations. Some numerical examples are computed to show the effectiveness of the proposed methods.
2 Hammerstein integral equations
In this section, we will transform the nonlinear second-order ordinary differential equation (1) with three-point boundary value conditions in (2)-(5) into Hammerstein integral equations. Then the existence and uniqueness of solutions for the obtained Hammerstein integral equations will be investigated.
2.1 Transformations
In the following, we apply the integration method to get the following four theorems.
Theorem 1
When \((b-a)+k(\xi-a)\neq0\), the three-point boundary value problem
can be transformed into the Hammerstein integral equation as follows:
where
Proof
By integrating twice, from a to x, both sides of the differential equation in (8) with respect to x, one gets
Then the unknowns \(\varphi'(a)\) and \(\varphi(a)\) in (10) will be determined by using the boundary value conditions in (8), respectively. By setting \(x=b\) in (10) one obtains
Applying (11) to (10), it follows that
Now we let \(x=\xi\) in (10) and obtain
In virtue of (8), (12), and (13), it gives
Since the right side of equation (14) is a function with respect to x, we let
After some computations, one verifies that (14) and (15) can be rewritten as indicated in (9). □
Furthermore, one can obtain the following theorem and the proof has been omitted for saving space.
Theorem 2
Under the condition of \((b-a)+k(\xi-a)\neq0\), the three-point boundary value problem
is equivalent to the following Hammerstein integral equation:
where
In the following, by considering the boundary value conditions in the cases of III and IV, we give Theorems 3 and 4, respectively.
Theorem 3
In virtue of \(k+(b-a)\neq0\), the ordinary differential equation with three-point boundary value conditions
can be changed to the Hammerstein integral equation. That is, one has
where
Proof
Performing the integration procedures similar to those in Theorems 1 and 2, we obtain
It is assumed that \(x=a\) in (20) and one arrives at
Insertion of (21) into (20) yields
With the knowledge of \(\varphi(b)\) in (12), and using the boundary value conditions in (18), one further has
The proof is completed by rewriting (23) as (19). □
Theorem 4
If \((b-a)-k\neq0\), the following boundary value problem:
is equivalent to the Hammerstein integral equation as follows:
where
The proof can be completed similar to that of Theorem 3.
It is seen from Theorems 1-4 that the nonlinear second-order three-point boundary value problems have been transformed into Hammerstein integral equations. We remark that the integration method is uniform and enough to transform any nonlocal boundary value problem of ordinary differential equations into an integral equation [19, 31]. In the end, as a check, we consider the special case of \(g(x)=\alpha=\beta=0\) and obtain the boundary value problems as those in [32]. Based on Theorems 1-4, the solutions and Green’s functions of the nonlinear three-point boundary value problems in [32] can be determined easily.
2.2 Existence and uniqueness of solution
For a nonlinear equation, the fixed point theorems are always used to address the existence and uniqueness of solutions [33–35]. Here since the considered three-point boundary value problems have been transformed into the Hammerstein integral equations, it is natural to study the existence and uniqueness of solutions for the obtained Hammerstein integral equations. Moreover, it is seen that the existence and uniqueness of solutions for Hammerstein integral equations have been investigated widely such as those in the book [36] and the recent results on \(\mathbf{L}_{1}\) spaces [37]. In the present paper, for the obtained Hammerstein integral equations, we will use the Schauder fixed point theorem to address the existence of the solutions, and apply the Banach fixed point theorem to investigate the uniqueness of the solutions. The obtained results are related to the considered nonlinear boundary value problems of ordinary differential equations.
Now it is convenient to rewrite equations (9), (17), (19), and (25) in the following general form, namely:
where \(\rho_{i}(x)\) (\(i=1,2,3,4\)) are the linear functions with respect to x. It is easy to see that \(K_{i}(x,t)\) and \(\rho_{i}(x)\) are continuous on \([a,b;a,b]\) and \([a,b]\), respectively. Furthermore, we assume that
Based on the Schauder fixed point theorem, we first give the following theorem to address the existence of the solutions in (26).
Theorem 5
It is assumed that \(S=\{\varphi|\varphi\in\mathbf{L}_{2}[a,b],\|\varphi \|\le M\}\). \(\forall\varphi\in S\), one has
\(\forall\varepsilon>0\), \(\varphi_{1},\varphi_{2}\in S\), \(\exists\delta (\varepsilon)>0\), when \(\|\varphi_{1}-\varphi_{2}\|<\delta(\varepsilon)\), it gives
If \(A_{i}B(b-a)+D_{i}\sqrt{b-a}< M\), the nonlinear integral equation (26) at least has a solution in S.
Proof
Obviously S is a closed convex set. For convenience, we define the integral operator T from \(\mathbf{L}_{2}[a,b]\) to \(\mathbf {L}_{2}[a,b]\) such that
First, one can see that T is a mapping from S to S. In fact, for \(\|\varphi\|\le M\), we have
then
Second, we prove that T is continuous. \(\forall\varepsilon >0\), \(\varphi_{1},\varphi_{2}\in S\), \(\exists\delta(\varepsilon)>0\), such that when \(\|\varphi_{1}-\varphi_{2}\|<\delta(\varepsilon)\), it follows that
and
In the end, let us prove that \(\mathbf{T}(S)\) is relatively compact. Since \(K_{i}(x,t)\) and \(\rho_{i}(x)\) are continuous on \([a,b;a,b]\) and \([a,b]\), respectively, they are uniformly continuous. Consequently, \(\forall\varepsilon>0\), \(\varphi\in S\), \(\exists\eta(\varepsilon)>0\), for \(|x_{1}-x_{2}|<\eta(\varepsilon)\), one has
It is seen that \(\mathbf{T}(S)\) is uniformly bounded and equicontinuous. According to Ascoli-Arzela theorem [38], \(\{ \mathbf{T}\varphi(x)\}\) has a subsequence with uniform convergence, so \(\mathbf{T}(S)\) is relatively compact.
By using the Schauder fixed point theorem, there exists at least a point \(\varphi\in S\) such that \(\mathbf{T}\varphi=\varphi\). □
In addition, we can change some conditions in Theorem 5 to obtain the following corollary.
Corollary 1
Let \(S=\{\varphi|\varphi\in\mathbf{L}_{2}[a,b],\|\varphi\|\le M\}\), and
\(\forall\varepsilon>0\), \(\varphi_{1},\varphi_{2}\in S\), \(\exists\delta (\varepsilon)>0\), when \(\|\varphi_{1}-\varphi_{2}\|<\delta(\varepsilon)\), one has
If \(A_{i}(E+F)(b-a)+D_{i}<\frac{M}{\sqrt{b-a}}\), the nonlinear integral equation (26) at least has a solution in S.
Furthermore, strengthening the conditions of the nonlinear term \(\psi (x,\varphi(x))\), we have the uniqueness theorem of solution in Banach spaces.
Theorem 6
Suppose that \(g(x)\in\mathbf{L}_{2}[a,b]\), and \(\psi(x,\varphi(x))\) satisfies the following Lipschitz condition:
with the Lipschitz constant \(L>0\). Moreover, one has the constraint conditions as
where N and \(C_{i}\) are positive constants. When \(LC_{i}<1\), the nonlinear integral equation (26) has a unique solution in \(\mathbf{L}_{2}[a,b]\).
Proof
Since \(g(x)\in\mathbf{L}_{2}[a,b]\), we can easily get \(f_{i}(x)\in\mathbf {L}_{2}[a,b]\). The kernel \(K_{i}(x,t)\) is a polynomial function with respect to x and t, so \(K_{i}(x,t)\in\mathbf{L}_{2}[a,b;a,b]\). Assume that T is an operator form \(\mathbf{L}_{2}[a,b]\) to \(\mathbf{L}_{2}[a,b]\), and \(\mathbf{T}\varphi=f_{i}-\boldsymbol{\kappa}_{i}\varphi \), where
By using
it is seen that \({\boldsymbol{\kappa}_{i}}\) is a bounded operator form \(\mathbf{L}_{2}[a,b]\) to \(\mathbf{L}_{2}[a,b]\). \(\forall\varphi_{1},\varphi_{2}\in\mathbf{L}_{2}[a,b]\), we further have
Then it gives
and
When \(LC_{i}<1\), T is a contraction operator. According to the fixed point theorem in Banach spaces, one can see that \(\mathbf{T}\varphi=\varphi\) has a unique solution in \(\mathbf{L}_{2}[a,b]\). □
3 Approximation of the solution
In practical applications, of much interest is how to obtain the solutions except for the existence of the solutions. However, the closed-form solutions of the Hammerstein integral equations in (9), (17), (19), and (25) cannot be determined easily due to the complexity of the kernels. Thus it is interesting to obtain numerical solutions of Hammerstein integral equations and many methods have been proposed [8, 39–46]. Moreover, it is noted that a simple Taylor-series expansion method has been proposed in [47] and modified in [48, 49] for numerically solving linear Fredholm integral equations of the second kind. Recently, by using the idea of piecewise approximation, the simple Taylor-series expansion method has been further modified in [28]. Here the proposed method in [28] is further extended and applied to solve the nonlinear integral equation of Hammerstein type. The convergence and error estimate of the approximate solution will be made. Moreover, it is seen from Theorems 5 and 6 that a solution in \(\mathbf{L}^{2}[a,b]\) is only determined by using the given conditions. Based on the proposed numerical method, the solution \(\varphi(x)\) should have more smoothing property and here it is assumed \(\varphi(x)\in\mathbf{C}^{n+1}[a,b]\) (\(n\geq0\)). Indeed, the case of \(\varphi(x)\in\mathbf{C}^{n+1}[a,b]\) (\(n\geq0\)) is important in practical applications. Two examples will be given in Section 4 to verify the extended numerical method by comparing a difference format.
3.1 Constructing the approximate solution
Generally, we write the Hammerstein integral equation as
where \(K(x,t)\) and \(f(x)\) are known functions. It is convenient to define the integral operator κ as
and it is compact from \(C^{(n+1)}[a,b]\) into \(C^{(n+1)}[a,b]\) for \(K(x,t)\in C^{(n+1)}[a,b; a,b]\).
Similar to those in [28], we choose a series of quadrature points as \(a=x_{0}< x_{1}<\cdots< x_{m}=b\) for \(m\geq1\). The integral operator κ can be further expressed as the following sum:
For the simplified case, the equidistant quadrature points are always chosen as
By letting \(t=x_{q}+hs\), equation (29) is reexpressed as
Now it is assumed that \(\psi(x_{q}+hs,\varphi(x_{q}+hs))\) can be expanded as the following Taylor series:
where \(R_{n}(\theta_{q},h,s)\) denotes the Lagrange remainder,
By eliminating the Lagrange remainder, the operator \((\boldsymbol{\kappa}\varphi)(x)\) can be approximated by using
Moreover, we suppose that
where the superscript \((i)\) denotes the ith-order differentiation with respect to x and \((\boldsymbol{\kappa}\varphi)^{(0)}(x)=(\boldsymbol{\kappa}\varphi)(x)\). Making use of (32), \((\boldsymbol{\kappa}\varphi)^{(i)}(x)\) can be approximated by
Now we further have the following theorem.
Theorem 7
Assume that one has the following conditions:
where \(i=0,1,2,\ldots,n\). The sequence \((\boldsymbol{\kappa}_{n}\varphi )^{(i)}(x)\) is convergent, namely
Proof
Applying equations (35) and (36), we get
From (38) it follows that \(\|(\boldsymbol{\kappa}_{n}\varphi)^{(i)}(x)-(\boldsymbol{\kappa}\varphi )^{(i)}(x)\|_{\infty}\to0\) with \(n\rightarrow+\infty\) and the proof is completed. □
In the end, let us give the approximate solution of Hammerstein integral equations. From equations (34) and (36), the discretization format of the derivatives of the Hammerstein integral equation (27) can be expressed as
with \(i=0,1,\ldots,n\) and \(l=0,1,\ldots,(m-1)\). Hereafter \(\frac{d^{j} \psi}{d y^{j}}|_{(x_{q},\varphi^{(j)}_{q})}\) means that the variable y equals \(x_{q}\) and the exact value \(\varphi ^{(j)}(x_{q})\) is replaced by the approximate one \(\varphi^{(j)}_{q}\). Once the solution of the nonlinear system (39) is given, the approximate solution of \(\varphi(x)\) can be further constructed as
where \(a\le x \le b\).
As shown in [28], one can see from equation (40) that the approximate solution has two parameters. The proposed method is based on the discretization points \(x_{q}\) (\(q=0,1,\ldots, m-1\)), which is different from the simple Taylor-series expansion method in [47–49]. The effectiveness and advantage of the new method will be shown in the given numerical examples of Section 4. Furthermore, in order to give the approximate solution in (40), the convergence of the nonlinear system (39) is requisite. Under some conditions, the nonlinear system (39) is convergent and it will be proved in the next subsection about the error estimate of the approximate solution.
On the other hand, it should be pointed out that when the proposed method is applied to solve the Hammerstein integral equations in (9), (17), (19), and (25), the derivatives \(K_{x}^{(i)}(x,t)\) (\(i=1,2,\ldots,n\)) for \(x=t\) in (39) must be dealt with again. The reason is based on the fact that the derivatives \(\partial ^{i}K_{j}(x,t)/ \partial x^{i}\) (\(i=1,2,\ldots,n\); \(j=1,2,3,4\)) for \(x=t\) are not existing. As shown in [50], for the practical computations, it is reasonable to adopt the following method:
In addition, we can rewrite equation (26) as
where \(\rho_{i}(x)\) (\(i=1,2,3,4\)) are the linear functions with respect to x. For generic nonlinear functions \(\psi(t,\varphi(t))\) and \(g(t)\), the proposed numerical methods can be used similarly for equation (42).
3.2 Convergence and error estimate
From the viewpoint of mathematical theory and practical applications, the convergence and error estimate of the approximate solution are all important. For the approximation method, we have the following theorem.
Theorem 8
It is assumed that \(d^{j}\psi(y,\varphi(y))/dy^{j}\) satisfies the Lipschitz conditions as follows:
with the Lipschitz constants \(L_{\nu}>0\) and \(j=0,1,\ldots,n\). One further has the following conditions:
where \(\bar{L}=\max_{\nu=0,1,\ldots,j}L_{v}\) and \(i=0,1,2,\ldots,n\). The approximate solution \(\varphi_{m,n}(x)\) in (40) is convergent to the exact solution \(\varphi(x)\). That is, we get
and
Moreover, the following error estimate can be obtained:
where \(\bar{M}=\max_{i=0,1,\ldots,n}M_{i}\).
Proof
Equation (39) can be further rewritten as
where
with
for \(i=0,1,\ldots,n\) and \(l=0,1,\ldots,(m-1)\).
On the other hand, application of equations (27) and (32) leads to
where
with
Then it is further found that
where \(\bar{M}=\max_{i=0,1,\ldots,n}M_{i}\) and \(\bar{L}=\max_{\nu =0,1,\ldots,j}L_{v}\).
Now let us assume that
Based on (43) and (44), it follows that
Furthermore, one has
One can see from (48) that the nonlinear system (39) is convergent with \(n\rightarrow+\infty\) and \(h\rightarrow0\), respectively.
In the following, we define a sequence of numerical integration operators \(\bar{\boldsymbol{\kappa}}_{n}\) such as
From Theorem 7, it is seen that
It is found from (50) that when \(n\rightarrow+\infty\) or \(m\rightarrow+\infty\) (i.e. \(h\rightarrow0\)), one always has \(\|\varphi(x)-\varphi_{m,n}(x)\|_{\infty}\rightarrow0\). This completes the proof. □
As shown in Theorem 8, one can choose a pair of feasible values for m (i.e. h) and n to obtain a good approximation of the exact solution. The above observations will be further verified by using the numerical examples in the next section.
4 Numerical results
In order to show the effectiveness of the proposed methods, we give two numerical examples corresponding to cases I and IV, respectively. The existence and uniqueness of the solution will be considered, and the approximate solution will be calculated numerically. All the computations are made by using the programming language of MATLAB (R2014).
Example 1
Assume that a second-order three-point boundary value problem is given as
where \(x\in[0,1]\), \(\varphi(x)\in\mathbf{C}^{\infty}[0,1]\) and \(|\varphi (x)|<3\). The exact solution is \(\varphi(x)=xe^{x}\).
According to Theorem 1, the nonlinear ordinary differential equation with three-point boundary value conditions in (51) can be transformed into the following Hammerstein integral equation:
where
We first prove the uniqueness of the solution and choose the Lipschitz constant as
Then one obtains
Since \(LC_{1}<6\sqrt{0.0050347}=0.42573<1\), the nonlinear three-point boundary value problem has a unique solution according to Theorem 6.
Now the proposed numerical method is applied to solve the obtained Hammerstein integral equation (52). It is seen that because \(f_{1}(x)\) in (52) is explicit, the reformulation of the problem in (42) is unnecessary and the proposed numerical method is carried out directly. Similar to those in [28], it is convenient to write the parameters m and n in vector style \((m,n)\). Figure 1 shows the exact solution and the approximate one with \((m,n)=(1,0)\). It is seen that the approximate solution is approached to the exact one. Moreover, the cases of \((m,n)=(2,2)\), \((m,n)=(2,4)\), \((m,n)=(4,2)\), and \((m,n)=(4,4)\) are chosen to compute, respectively. The absolute errors between the approximate and exact solutions are listed in Table 1. One can see from Table 1 that we have given a good approximation of the exact solution and the proposed numerical method is effective. On the other hand, the central difference format is used to compute the nonlinear boundary value problem and the obtained results are given in Table 1 for comparisons. A system of nonlinear algebraic equations is constructed and it is solved by using the Broyden iterative method [51]. The initial vectors are \({ \mathbf{x}_{0}}=[1,1,\ldots ,1]_{10\times1}^{T}\) and \({\mathbf{x}_{0}}=[1,1,\ldots,1]_{20\times 1}^{T}\) for the step \(h=0.1\) and \(h=0.05\), respectively. The results are obtained by iterating five times. The comparisons reveal that the proposed method is effective.
Example 2
Consider a second-order three-point boundary value problem as follows:
where \(x\in[0,1]\), \(\varphi(x)\in\mathbf{C}^{\infty}[0,1]\) and \(|\varphi (x)|\le1\). The exact solution is \(\varphi(x)=x^{2}\).
Making use of Theorem 4, the boundary value problem (53) can be transformed into the following Hammerstein integral equation:
where
If we choose the Lipschitz constant as
it follows that
Under the consideration of \(LC_{4}\le e\sqrt{0.016840}=0.35275<1\), one can see that the three-point boundary value problem has a unique solution by using Theorem 6.
For numerical computations, as shown in (42), equation (54) should be rewritten as
In addition, the exact solution and the approximate one for \((m,n)=(1,0)\) are depicted in Figure 2. The absolute errors of the exact and approximate solutions by using \((m,n)=(2,1)\), \((m,n)=(2,2)\), \((m,n)=(4,1)\), and \((m,n)=(4,2)\), and the central difference format are shown in Table 2. It is found from Figure 2 and Table 2 that a good approximation solution is determined. When m or n is increasing, the absolute error \(|\varphi(x)-\varphi_{m,n}(x)|\) is decreasing. The observation is in accordance with the theoretical analysis in Theorem 8 and that in [28].
5 Conclusions
Four cases of nonlinear second-order three-point boundary value problems have been investigated and they are transformed into the Hammerstein integral equations by using the integration method. Based on the Schauder fixed point theorem, the sufficient conditions for the existence of the solutions have been given. The uniqueness of the solutions has been considered by using the Banach fixed point theorem. Furthermore, we have constructed the approximate solution of Hammerstein integral equations by applying a novel numerical method, which depends on the values of two parameters. The convergence and error estimate of the approximate solution have been made, and they show that one can get a good approximation of the exact solution by choosing a pair of the parameters. Two examples have been carried out numerically and the obtained results have revealed that the proposed methods are effective. In the future, the proposed method will be extended to solve nonlinear second-order differential equations with various nonlocal boundary value conditions.
References
Picone, M: Su un problema al contorno nelle equazioni differenziali lineari ordinarie del secondo ordine. Ann. Sc. Norm. Super. Pisa, Cl. Sci. 10, 1-95 (1908)
Boucherif, A: Nonlinear three-point boundary value problems. J. Math. Anal. Appl. 77, 577-600 (1980)
Štikonas, A: A survey on stationary problems, Green’s functions and spectrum of Sturm-Liouville problem with nonlocal boundary conditions. Nonlinear Anal., Model. Control 19, 301-334 (2014)
Conti, R: Recent trends in the theory of boundary value problems for ordinary differential equations. Boll. Unione Mat. Ital. 22, 135-178 (1967)
Ma, R: A survey on nonlocal boundary value problems. Appl. Math. E-Notes 7, 257-279 (2001)
Ntouyas, SK: Nonlocal initial and boundary value problems: a survey. In: Handbook of Differential Equations: Ordinary Differential Equations, vol. II, pp. 461-557. Elsevier, Amsterdam (2005)
Whyburn, WM: Differential equations with general boundary conditions. Bull. Am. Math. Soc. 48, 692-704 (1942)
Ascher, UM, Mattheij, RMM, Rusell, RD: Numerical Solution of Boundary Value Problems for Ordinary Differential Equations. SIAM, Philadelphia (1995)
Ma, R: Multiplicity results for a three-point boundary value problem at resonance. Nonlinear Anal. 53, 777-789 (2003)
Xu, X: Positive solutions for singular m-point boundary value problems with positive parameter. J. Math. Anal. Appl. 291, 352-367 (2004)
Xu, X: Multiplicity results for positive solutions of some semi-positone three-point boundary value problems. J. Math. Anal. Appl. 291, 673-689 (2004)
Yao, Q: Successive iteration and positive solution for nonlinear second-order three-point boundary value problems. Comput. Math. Appl. 50, 433-444 (2005)
Nieto, JJ: Existence of a solution for a three-point boundary value problem for a second-order differential equation at resonance. Bound. Value Probl. 2013, 130 (2013)
Liu, X, Qiu, J, Guo, Y: Three positive solutions for second-order m-point boundary value problems. Appl. Math. Comput. 156, 733-742 (2004)
Jiang, W, Guo, Y: Multiple positive solutions for second-order m-point boundary value problems. J. Math. Anal. Appl. 327, 415-424 (2007)
Liu, H, Ouyang, Z: Existence of solutions for second-order three-point integral boundary value problems at resonance. Bound. Value Probl. 2013, 197 (2013)
Karakostas, GL, Tsamatos, PC: Existence of multiple positive solutions for a nonlocal boundary value problem. Topol. Methods Nonlinear Anal. 19, 109-121 (2002)
Karakostas, GL, Tsamatos, PC: Multiple positive solutions of some Fredholm integral equations arisen from nonlocal boundary-value problems. Electron. J. Differ. Equ. 2002, 30 (2002)
Webb, JRL, Infante, G: Positive solutions of nonlocal boundary value problems: a unified approach. J. Lond. Math. Soc. 74, 673-693 (2006)
Ma, R: Multiplicity of positive solutions for second-order three-point boundary value problems. Comput. Math. Appl. 40, 193-204 (2000)
He, X, Ge, W: Triple solutions for second-order three-point boundary value problems. J. Math. Anal. Appl. 268, 256-265 (2002)
Guo, Y, Shan, W, Ge, W: Positive solutions for second-order m-point boundary value problems. J. Comput. Appl. Math. 151, 415-424 (2003)
Han, X: Positive solutions for a three-point boundary value problem at resonance. J. Math. Anal. Appl. 336, 556-568 (2007)
Sun, Y, Liu, L, Zhang, J, Agarwal, RP: Positive solutions of singular three-point boundary value problems for second-order differential equations. J. Comput. Appl. Math. 230, 738-750 (2009)
Webb, JRL: Positive solutions of some three point boundary value problems via fixed point index theory. Nonlinear Anal. 47, 4319-4332 (2001)
Bai, D, Feng, H: Eigenvalue for a singular second order three-point boundary value problem. J. Appl. Math. Comput. 38, 443-452 (2012)
Webb, JRL: Existence of positive solutions for a thermostat model. Nonlinear Anal., Real World Appl. 13, 923-938 (2012)
Zhong, XC, Huang, QA: Approximate solution of three-point boundary value problems for second-order ordinary differential equations with variable coefficients. Appl. Math. Comput. 247, 18-29 (2014)
Yang, X, Wang, Z, Shen, J: Existence of solution for a three-point boundary value problem for a second-order impulsive differential equation. J. Appl. Math. Comput. 47, 49-59 (2015)
Verma, AK, Singh, M: A note on existence results for a class of three-point nonlinear BVPs. Math. Model. Anal. 20(4), 457-470 (2015)
Corduneanu, C: Integral Equations and Applications. Cambridge University Press, Cambridge (1991)
Zhao, ZQ: Solutions and Green’s functions for some linear second-order three-point boundary value problems. Comput. Math. Appl. 56, 104-113 (2008)
Granas, A, Dugundji, J: Fixed Point Theory. Springer, New York (2003)
Wang, F, Zhang, F: Existence of positive solutions of Neumann boundary value problem via a cone compression-expansion fixed point theorem of functional type. J. Appl. Math. Comput. 35, 341-349 (2011)
Zhang, F, Lu, H, Wang, F: Existence and multiplicity of positive solutions for second-order self-adjoint boundary value problem with integral boundary conditions at resonance. Fixed Point Theory 13(2), 669-680 (2012)
Guo, DJ, Sun, JX: Nonlinear Integral Equation. Shandong Science and Technology Press, Jinan (1987) (in Chinese)
Latrach, K, Taoudi, MA: Existence results for a generalized nonlinear Hammerstein equation on \(\mathbf{L}_{1}\) spaces. Nonlinear Anal. 66, 2325-2333 (2007)
Trench, WF: Introduction to real analysis. Library of Congress Cataloging-in-Publication Data (2010)
Abdou, MA, El-Borai, MM, El-Kojok, MM: Toeplitz matrix method and nonlinear integral equation of Hammerstein type. J. Comput. Appl. Math. 223, 765-776 (2009)
Atkinson, KE: A survey of numerical methods for solving nonlinear integral equations. J. Integral Equ. Appl. 4, 15-46 (1992)
Golberg, MA: Solution Methods for Integral Equations: Theory and Applications. Plenum, New York (1979)
Chidume, CE, Djitté, N: Iterative approximation of solutions of nonlinear equations of Hammerstein type. Nonlinear Anal. 70, 4086-4092 (2009)
Chidume, CE, Djitté, N: An iterative method for solving nonlinear integral equations of Hammerstein type. Appl. Math. Comput. 219, 5613-5621 (2013)
Chidume, CE, Osilike, MO: Iterative solution of nonlinear integral equations of the Hammerstein type. J. Niger. Math. Soc. 11, 9-19 (1992)
Kaneko, H, Noren, RD, Novaprateep, B: Wavelet applications to the Petrov-Galerkin method for Hammerstein equations. Appl. Numer. Math. 45, 255-273 (2003)
Moore, C: The solution by iteration of nonlinear equations of Hammerstein type. Nonlinear Anal. 49, 631-642 (2002)
Ren, Y, Zhang, B, Qiao, H: A simple Taylor-series expansion method for a class of second kind integral equations. J. Comput. Appl. Math. 110, 15-24 (1999)
Huabsomboona, P, Novaprateep, B, Kaneko, H: On Taylor-series expansion methods for the second kind integral equations. J. Comput. Appl. Math. 234, 1466-1472 (2010)
Li, XF: Approximate solution of linear ordinary differential equations with variable coefficients. Math. Comput. Simul. 75, 113-125 (2007)
Lanczos, C: Linear Differential Operator. Van Nostrand, New York (1961)
Broyden, CG: A class of methods for solving nonlinear simultaneous equations. Math. Comput. 19, 577-593 (1965)
Acknowledgements
The authors would like to thank the Associate Editor and the reviewers for the valuable suggestions in improving the paper. The work was supported by the National Natural Science Foundation of China (No. 11362002), the Natural Science Foundation of Guangxi (No. 2016GXNSFAA380012), the Innovation Project of Guangxi Graduate Education (No. YCSZ2015030), the project of outstanding young teachers’ training in higher education institutions of Guangxi, and the project of Guangxi Colleges and Universities Key Laboratory of Mathematics and Its Applications.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
XCZ and QA made equal contributions to theory, computations, and writing of the article. HMW made main contributions to the numerical computations by using a difference format. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Zhong, XC., Huang, QA. & Wei, HM. Existence and approximation of solution for a nonlinear second-order three-point boundary value problem. Bound Value Probl 2016, 169 (2016). https://doi.org/10.1186/s13661-016-0678-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-016-0678-4