JAMJournal of Applied Mathematics1687-00421110-757XHindawi Publishing Corporation35925110.1155/2012/359251359251Research ArticleOn Six Solutions for m-Point Differential Equations System with Two Coupled Parallel Sub-Super SolutionsLiuJian1,2SuHua2WangShuli3ChenRudong1School of Mathematical SciencesQufu Normal UniversityQufu, Shandong 273165Chinaqfnu.edu.cn2School of Mathematics and Quantitative EconomicsShandong University of Finance and EconomicsJinan, Shandong 250000Chinasufe.edu.cn3College of ScienceYanshan University at QinhuangdaoQinhuangdao 066004Chinaysu.edu.cn201222122011201215082011051020112012Copyright © 2012 Jian Liu 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.

Under the assumption of two coupled parallel subsuper solutions, the existence of at least six solutions for a kind of second-order m-point differential equations system is obtained using the fixed point index theory. As an application, an example to demonstrate our result is given.

1. Introduction

In this paper, we consider the following second-order m-point boundary value problems of nonlinear equations system -φ′′(t)=f1(φ(t))+f2(ψ(t)),  t[0,1],-ψ′′(t)=g1(ψ(t))+g2(φ(t)),t[0,1],φ(0)=0,φ(1)=i=1m-2αiφ(ξi),ψ(0)=0,ψ(1)=i=1m-2αiψ(ξi), where fi,gi:11(i=1,2) are continuous and αi,ξi satisfying

i=1m-2αi(0,1) with αi(0,+) for i=1,2,,m-2 and 0<ξ1<ξ2<<ξm-2<1.

Multipoint boundary value problems arise in many applied sciences for example, the vibrations of a guy wire composed of N parts with a uniform cross-section throughout, but different densities in different parts can be set up as a multipoint boundary value problems (see ). Many problems in the theory of elastic stability can be modelled by multipoint boundary value problems (see ). The study of multipoint boundary value problems for linear second-order ordinary differential equations was initiated by Il'in and Moiseev . Subsequently, Gupta  studied certain three-point boundary value problems for nonlinear second-order ordinary differential equations. Since then, the solvability of more general nonlinear multipoint boundary value problems has been discussed by several authors using various methods. We refer the readers to  and the references therein.

In the recent years, many authors have studied existence and multiplicity results for solutions of multipoint boundary value problems via the well-ordered upper and lower solutions method, see [8, 13, 14] and the references therein. However, only in very recent years, some authors considered the multiplicity of solutions under conditions of non-well-ordered upper and lower solutions. For some abstract results concerning conditions of non-well-ordered upper and lower solutions, the readers are referred to recent works .

In , Xu et al. considered the following second-order three-point boundary value problem y′′(t)+f(t,y)=0,t[0,1],y(0)=0,y(1)-αy(η)=0, where 0<η<1,  0<α<1,  fC([0,1]×1,1). He obtained the following result. First, let us give the following condition (H0) to be used later.

(H0) There exists M>0 such that f(t,x2)-f(t,x1)-M(x2-x1),t[0,1],x2x1. Let the function e be e=e(t)=t for t[0,1].

Theorem 1.1.

Suppose that (H0) holds, u1 and u2 are two strict lower solutions of (1.2), v1 and v2 are two strict upper solutions of (1.2), and u1<v1,  u2<v2,  u2v1,  u1v2. Moreover, assume -ξ0eu2-u1ξ0e,-ξ0ev2-v1ξ0e, for some ξ0>0. Then, the three-point boundary value problem (1.2) has at least six solutions.

We would also like to mention the result of Yang , in . Yang studied the following integral boundary value problem -(au)+bu=g(t)f(t,u),(cosγ0)u(0)-(sinγ0)u(0)=01u(τ)dα(τ),(cosγ1)u(1)+(sinγ1)u(1)=01u(τ)dβ(τ), where γ0[0,π/2] and γ1[0,π/2], 01u(τ)dα(τ) and 01u(τ)dβ(τ) denote the Riemann-Stieltjes integrals of u with respect to α and β, respectively. Some sufficient conditions for the existence of either none, or one, or more positive solutions of the problem (1.5) were established. The main tool used in the proofs of existence results is a fixed point theorem in a cone, due to Krasnoselskii and Zabreiko.

At the same time, we note that Webb and Lan  have considered the first eigenvalue of the following linear problem u′′(t)+λu(t)=0,0<t<1,u(0)=0,  u(1)=i=1m-2αiu(ηi), they also investigated the existence and multiplicity of positive solutions of several related nonlinear multipoint boundary value problems. Furthermore, Ma and O'Regan  studied the spectrum structure of the problem (1.6), and the authors obtained the concrete computational method and the corresponding properties of real eigenvalue of (1.6) by constructing an auxiliary function. Their work is very fundamental to further study for multipoint boundary value problems. By extending and improving the work in , Rynne  showed that the associated Sturm-Liouville problem consisting of (1.6) has a strictly increasing sequence of simple eigenvalues {λn}n=0 with eigenfunctions ϕn(t)=sin(λnt).

Very recently, Kong et al.  were concerned with the general boundary value problem with a variable wu′′(t)+w(t)f(u)=0,t(a,b),cosαu(a)-sinαu(a)=0,α[0,π),u(b)=i=1m-2kiu(ηi) By relating (1.7) to the eigenvalues of a linear Sturm-Liouville problem with a two-point separated boundary condition, the existence and nonexistence of nodal solutions of (1.7) were obtained. We also point out that Webb  made the excellent remark on some existence results of symmetric positive solutions obtained in some recent papers and the author also corrected the values of the principle eigenvalue previously given in some examples.

In this paper, by means of two coupled parallel subsuper solutions, we obtain some sufficient conditions for the existence of six solutions for (1.1) and our main tool is based on the fixed point index theory. At the end of this paper, we will give an example which illustrates that our work is true. Our method stems from the paper .

2. Preliminaries and a Lemma

In the section, we shall give some preliminaries and a lemma which are fundamental to prove our main result.

Let E be an ordered Banach space in which the partial ordering ≤ is induced by a cone PE. A cone P is said to be normal if there exists a constant N>0, such that θxy implies xNy, the smallest N is called the normal constant of P. P is called solid, if intP, that is, P has nonempty interior. Every cone P in E defines a partial ordering in E given by xy if and only if y-xP. If xy and xy, we write x<y; if cone P is solid and y-x intP, we write xy. P is called total if E=P-P¯. Let B:EE be a bounded linear operator. B is said to be positive if B(P)P. An operator A is strongly increasing, that is, x<y implies AxAy. If A is a linear operator, A is strongly increasing implying A is strongly positive.

Let E be an ordered Banach space, P a total cone in E, the partial ordering ≤ induced by P. B:EE is a positive completely continuous linear operator. Let r(B)>0 the spectral radius of B, B* the conjugated operator of B, and P* the conjugated cone of P. Since PE is a total cone (i.e., E=P-P¯), according to the famous Krein-Rutman theorem (see ), we infer that if r(B)0, then there exist φ¯P{θ} and g*P*{θ}, such that Bφ¯=r(B)φ¯,B*g*=r(B)g*. Fixed φ¯P{θ},g*P*{θ} such that (2.1) holds. For δ>0, letP(g*,δ)={xP,  g*(x)δx}, then P(g*,δ) is also a cone in E. One can refer  for definition and properties about the cones.

Definition 2.1 (see [<xref ref-type="bibr" rid="B21">29</xref>]).

Let B be a positive linear operator. The operator B is said to satisfy condition H, if there exist φ¯P{θ},g*P*{θ}, and δ>0 such that (2.1) holds, and B maps P into P(g*,δ).

Lemma 2.2 (see [<xref ref-type="bibr" rid="B14">10</xref>]).

Suppose that d=1-i=1m-2αi0. Then, the BVP -u′′(t)=0,t(0,1),u(0)=0,u(1)=i=1m-2αiu(ξi), has Green's function G(t,s)=G̃(t,s)+i=1m-2αiG̃(ξi,s)1-i=1m-2αi, where G̃(t,s)={1-t,0st1,1-s,0ts1.

For convenience, we list the following hypotheses which will be used in our main result.

fi,  gi(i=1,2) are strictly increasing;

there exist constants k>0,  l>0 and D>0 such that

|f1(±k)±l|<N-1k,

|g1(±k)±D|<N-1k,

|f2(±k)|l, and (iv) |g2(±k)|D, where N=maxt[0,1]01G(t,s)ds;

there exist constants 0<c1<k,  -k<c2<0,  0<c3<k,  -k<c4<0, such that, for all t[0,1], we have

c1<01G(t,s)f1(c1)ds-Nl,

c2>01G(t,s)g1(c2)ds+ND,

c3<01G(t,s)g1(c3)ds-ND, and

c4>01G(t,s)f1(c4)ds+Nl;

lim|φ|+(f1(φ)+f2(ψ))/φ2λ1 uniformly for ψ,

lim|Ψ|+(g1(ψ)+g2(φ))/ψ2λ1 uniformly for φ, where λ1 is the first eigenvalue of the following boundary value problem:-u′′(t)=λu,  t(0,1),u(0)=0,u(1)=i=1m-2αiu(ξi).

It is well known that λ1=r-1(H), where linear operator H:C[0,1]C[0,1] is defined as Hu(t)=01G(t,s)u(s)ds.

3. Main ResultsTheorem 3.1.

Assume (H0), (H1)(H4) hold, then BVP (1.1) has at least six distinct continuous solutions.

Proof.

It is easy to check that BVP (1.1) is equivalent to the following integral equation systems: φ(t)=01G(t,s)[f1(φ(s))+f2(ψ(s))]ds,ψ(t)=01G(t,s)[g1(ψ(s))+g2(φ(s))]ds, where G(t,s) is defined as in Lemma 2.2. By (H0), we know that G(t,s)0,  for  all  t,s[0,1].

Let E=C[0,1]×C[0,1], define the norm in E as (φ,ψ)=φ+ψ. Then, E is a Banach space with this norm. Let P={(φ,ψ)Eφ(t)0,ψ(t)0, for  all  t[0,1]},Q={φC[0,1]φ(t)0,  for  all    t[0,1]}. Then, P=Q×Q is a normal and solid cone. Set T:EE, such that T(φ,ψ)=(01G(t,s)[f1(φ(s))+f2(ψ(s))]ds,01G(t,s)[g1(ψ(s))+g2(φ(s))]ds), it is clear that the solutions of (1.1) are equivalent to the fixed points of T.

Set 0<ξ1λ1, let f(φ,ψ)=(f1(φ)+f2(ψ),g1(ψ)+g2(φ))λ1+ξ1,K(φ,ψ)=((λ1+ξ1)Hφ,(λ1+ξ1)Hψ)=(H1φ,H1ψ), where H1=(λ1+ξ1)H, then T=Kf. It is easy to see that H is a strongly positive completely continuous operator, and it follows from G(t,s)0 and the continuity of G(t,s) that K is a strongly positive completely continuous operator. Since fi,  gi:R1R1(i=1,2) are strictly increasing continuous functions, we know that f is a strictly increasing continuous bounded operator. By T=Kf, we can prove that T is completely continuous. We infer from the increasing properties of K and f that T is increasing.

Let φ1c1,  ψ1-k,    φ2k,    ψ2c2,φ3-k,    ψ3c3,    φ4c4,    ψ4k, then (φi,ψi)(i=1,2,3,4) satisfy (φ1,ψ1)<(φ2,ψ2),(φ3,ψ3)<(φ4,ψ4),(φ1,ψ1)(φ4,ψ4),(φ3,ψ3)(φ2,ψ2). By (H2)(iii) and (H3)(i), we have 01G(t,s)[f1(φ1)+f2(ψ1)]ds01G(t,s)f1(c1)ds-Nl>c1=φ1. It follows from (H2)(ii),(iv), and the increasing property of g2 that 01G(t,s)[g1(ψ1)+g2(φ1)]ds>(-k)N-101G(t,s)ds(-k)N-1N=-k=ψ1. Equations (3.2), (3.5), and (3.6) imply that (φ1,ψ1)<T(φ1,ψ1). Similarly, by (H1)(H3), we obtain T(φ2,ψ2)<(φ2,ψ2),(φ3,ψ3)<T(φ3,ψ3),T(φ4,ψ4)<(φ4,ψ4).

By [20, Lemma 3], we get that H1 satisfies condition H. Therefore, there exist j0*Q*{θ},  δ>0, such that H1*j0*=r(H1)j0*,j0*(H1φ)δH1φ,φQ. By the definition of spectral radius of completely continuous operator, we have r(K)=r(H1), and combining (3.9), we infer that H1*j0*=r(K)j0*. Let j*((φ,ψ))=j0*(φ)+j0*(ψ),  for  all  (φ,ψ)E, then j*P*{θ}. According to the proof in , we can get that K satisfies condition H.

By condition (H4), we obtain that there exists C>0, such that f1(φ)+f2(ψ)(λ1+ξ1)φ,  φC,ψR,  f1(φ)+f2(ψ)(λ1+ξ1)φ,φ-C,ψR,  g1(ψ)+g2(φ)(λ1+ξ1)ψ,ψC,φR,  g1(ψ)+g2(φ)(λ1+ξ1)ψ,ψ-C,φR. Equations (3.12)–(3.15) imply (f1(φ)+f2(ψ),g1(ψ)+g2(φ))(λ1+ξ1)(φ,ψ),  φ,ψC,(f1(φ)+f2(ψ),g1(ψ)+g2(φ))(λ1+ξ1)(φ,ψ),φ,ψ-C. Since f1(φ),g1(ψ) are continuous in [-(c1+k),C], so they are bounded, then there exists h>0 such that f1(φ)-h,g1(ψ)-h,-(φ1,ψ1)=-(c1+k)φ,ψC. By virtue of (3.18) and the increasing properties of f2 and g2, one shows (f1(φ)+f2(ψ),g1(ψ)+g2(φ))(-h-|f2(-(φ1,ψ1))|,-h-|g2(-(φ1,ψ1))|),-(φ1,ψ1)φ,ψC.

In addition, if φ,ψ satisfy φC,  -(φ1,ψ1)ψC, then it follows from (3.12), (3.18), and the increasing property of g2 that (f1(φ)+f2(ψ),g1(ψ)+g2(φ))((λ1+ξ1)φ,-h+g2(-(φ1,ψ1)))(λ1+ξ1)(φ,ψ)-((λ1+ξ1)C+h+|g2(-(φ1,ψ1))|,(λ1+ξ1)C+h+|g2(-(φ1,ψ1))|)=(λ1+ξ1)(φ,ψ)-(d1,d1), where d1=(λ1+ξ1)C+h+|g2(-(φ1,ψ1))|. Similarly, if φ,  ψ satisfy ψC,  -(φ1,ψ1)φC, then combining the increasing property of f2 with (3.14) and (3.18), we know that (f1(φ)+f2(ψ),g1(ψ)+g2(φ))(λ1+ξ1)(φ,ψ)-(d2,d2), where d2=(λ1+ξ1)C+h+|f2(-(φ1,ψ1))|. Let d=max{d1,d2}. By (3.21) and (3.22), we get that if φC,  -(φ1,ψ1)ψC or ψC,  -(φ1,ψ1)φC, it is obvious that (f1(φ)+f2(ψ),g1(ψ)+g2(φ))(λ1+ξ1)(φ,ψ)-(d,d). It follows from (3.16), (3.19), and (3.23) that (f1(φ)+f2(ψ),g1(ψ)+g2(φ))(λ1+ξ1)(φ,ψ)-(d,d),φ,ψ-(φ1,ψ1). In a similar way, from (3.12) and (3.14), we can show that there exists e>0 such that (f1(φ)+f2(ψ),g1(ψ)+g2(φ))(λ1+ξ1)(φ,ψ)-(e,e),  φ,ψ-(φ3,ψ3). Let a=max{d,e}. It follows from (3.24) and (3.25) that if φ,ψ-(φ1,ψ1) or φ,ψ-(φ3,ψ3), then (f1(φ)+f2(ψ),g1(ψ)+g2(φ))(λ1+ξ1)(φ,ψ)-(a,a). In a similar way, from (3.13) and (3.15), we can prove that there exists constant ã such that if φ,ψ(φ2,ψ2) or φ,ψ(φ4,ψ4), then (f1(φ)+f2(ψ),g1(ψ)+g2(φ))(λ1+ξ1)(φ,ψ)-(ã,ã). Let P((φ1,ψ1))={(u,v)E(u,v)(φ1,ψ1)}, P((φ3,ψ3))={(u,v)E(u,v)(φ3,ψ3)}, P((φ2,ψ2))={(u,v)E(u,v)(φ2,ψ2)}, P((φ4,ψ4))={(u,v)E(u,v)(φ4,ψ4)}, for  all  (φ,ψ)P((φ1,ψ1))P((φ3,ψ3)), then φ(t),ψ(t)-(φ1,ψ1) or φ(t),ψ(t)-(φ3,ψ3), for  all  t[0,1]; therefore, in virtue of expression of B,f, and (3.26), we have Kf(φ(t),ψ(t))K(φ(t),ψ(t))-(λ1+ξ1)a(λ1+ξ1)(01G(t,s)ds,01G(t,s)ds)=(r-1(K)+ϵ)K(φ(t),ψ(t))-(u1(t),v1(t)), where ϵ=1-r-1(K),  (u1(t),v1(t))=a(01G(t,s)ds,01G(t,s)ds). Since 0<ξ1λ1, one can show 0<ϵ=1-λ1λ1+ξ1=ξ1λ1+ξ1λ1λ1+ξ1=r-1(K).

This implies that there exist (u1,v1)E and 0<ξ2r-1(K) such that Kf(φ,ψ)(r-1(K)+ξ2)K(φ,ψ)-(u1,v1),(φ,ψ)P((φ1,ψ1))P((φ3,ψ3)). Similarly, we get by (3.27) that there exist (u2,v2)E and 0<ξ3r-1(K) such that Kf(φ,ψ)(r-1(K)+ξ3)K(φ,ψ)-(u2,v2),(φ,ψ)P((φ2,ψ2))P((φ4,ψ4)). We get by (3.7) that T:P((φ1,ψ1))P((φ1,ψ1)). Let E1={(u,v)E(φ1,ψ1)(u,v)(φ2,ψ2)}. Since P is normal, then E1 is bounded (see ). Choose M1>0 such that M1>max{sup(u,v)E1sup(u,v)E11ξ2δ(ξ2δ(φ1,ψ1)+ξ2δT(φ1,ψ1)+r-1(K)g*((u1,v1))-ξ2g*(T(φ1,ψ1))),sup(u,v)E1(u,v)+(φ1,ψ1)}. Let Ω1={(u,v)P((φ1,ψ1))(u,v)-(φ1,ψ1)<M1,(u,v)(φ3,ψ3)}, then E1Ω1 and Ω1 is a bounded open set. By the proof of Theorem 2.1 in , we can show that (u,v)-T(u,v)λK(u¯,v¯),λ0,  (u,v)Ω1, where (u¯,v¯) satisfies K(u¯,v¯)=r(K)(u¯,v¯).

Equation (3.34) implies that T has no fixed point on Ω1. It is easy to prove that P((φ1,ψ1)) is a retract of E, which together with (3.32) implies that the fixed point index i(T,Ω1,P((φ1,ψ1))) over Ω1 with respect to P((φ1,ψ1)) is well defined, and a standard proof yields i(T+sK(u¯,v¯),Ω1,P((φ1,ψ1)))=0. Set H¯(t,(u,v))=(1-t)T(u,v)+t(T(u,v)+sK(u¯,v¯)),  (t,(u,v))[0,1]×Ω1¯, then for any (u,v)Ω1¯,    t[0,1], we have H¯(t,(u,v))P((φ1,ψ1)). It follows from (3.34) that H(t,(u,v))(u,v),  for  all  (t,(u,v))[0,1]×Ω1, and by (3.35) and the homotopy invariance of the fixed point index, we get i(T,Ω1,P((φ1,ψ1)))=i(T+sK(u¯,v¯),Ω1,P((φ1,ψ1)))=0. Let W1={(u,v)P((φ1,ψ1))(u,v)(φ2,ψ2)}. By means of usual method (see ), we get that i(T,W1,P((φ1,ψ1)))=1. It is evident that A has no fixed point on W1, by (3.36), (3.37), and the additivity of the fixed point index, we have i(T,Ω1W1¯,P((φ1,ψ1)))=i(T,Ω1,P((φ1,ψ1)))-i(T,W1,P((φ1,ψ1)))=0-1=-1. Set E2={(u,v)E(φ3,ψ3)(u,v)(φ4,ψ4)}, and choose M2>0 such that M2>max{sup(u,v)E1sup(u,v)E11ξ2δ(ξ2δ(φ3,ψ3)+ξ2δT(φ3,ψ3)+r-1(K)g*((u1,v1))-ξ2g*(T(φ3,ψ3))),sup(u,v)E2(u,v)+(φ3,ψ3)}. Let Ω2={(u,v)P((φ3,ψ3))(u,v)-(φ3,ψ3)<M2,(u,v)(φ1,ψ1)},W2={(u,v)P((φ3,ψ3))(u,v)(φ4,ψ4)}. Similarly to the proof of (3.37) and (3.38), we get that i(T,W2,P((φ3,ψ3)))=1,i(T,Ω2W2¯,P((φ3,ψ3)))=-1. Choose M3,M4>0 such that M3>max{sup(u,v)E11ξ3δ(ξ3δ(φ2,ψ2)+ξ3δT(φ2,ψ2)+r-1(K)g*((u2,v2))-ξ3g*(T(φ2,ψ2)))sup(u,v)E1(u,v)+(φ2,ψ2)},M4>max{sup(u,v)E1  1ξ3δ(ξ3δ(φ4,ψ4)+ξ3δT(φ4,ψ4)+r-1(K)g*((u2,v2))-ξ3g*(T(φ4,ψ4))),sup(u,v)E2(u,v)+(φ4,ψ4)}. Set Ω3={(u,v)P((φ2,ψ2))(u,v)-(φ2,ψ2)<M3,(u,v)(φ4,ψ4)},Ω4={(u,v)P((φ4,ψ4))(u,v)-(φ4,ψ4)<M4,(u,v)(φ2,ψ2)},W3={(u,v)P((φ2,ψ2))(φ1,ψ1)(u,v)},W4={(u,v)P((φ4,ψ4))(φ2,ψ2)(u,v)}. By virtue of (3.31) and the same method as that for (3.34), we have (u,v)-T(u,v)-μK(u¯,v¯),μ0,  (u,v)Ω3Ω4. By (3.44), similarly to the proof of (3.38), we can prove that i(T,Ω3W3¯,P((φ2,ψ2)))=i(T,Ω4W4¯,P((φ4,ψ4)))=-1. Equations (3.37)–(3.41), (3.45) imply that T has at least six distinct fixed points, that is, the system of differential equations (1.1) has at least six solution in C[0,1]×C[0,1].

4. An Example

In this section, we present a simple example to explain our results.

Consider the following second-order three-point BVP for nonlinear equations system: -φ′′(t)=f1(φ(t))+f2(ψ(t)),t[0,1],-ψ′′(t)=g1(ψ(t))+g2(φ(t)),  t[0,1],φ(0)=0,φ(1)=12φ(14),ψ(0)=0,  ψ(1)=12ψ(14), where f1(φ)=(φ3/30000)+5φ1/3, g1(ψ)=(ψ3/30000)+4ψ1/5, f2(ψ)=(1/6)·1005+1/12arctan6ψ, g2(φ)=(1/16)·1003+(1/30)arctan7φ, α1=1/2,    ξ1=1/4,  N=maxt[0,1]01G(t,s)ds=31/32, fi,gi:11 are strictly increasing continuous functions, and condition (H1) is satisfied. Choose k=100,  l=5/12+π/24,D=5/16+π/60. Some direct calculations show |f1(±100)±l|1063104+51003+512+π24<3231100=N-1k,|g1(±100)±D|1003+41005+516+π60<34+452<3231100=N-1k,|f2(±100)|161005+112arctan6100<1652+112π2=512+π24=l,|g2(±100)|1161003+130arctan7100<1165+130π2=516+π60=D, Therefore, condition (H2) is satisfied.

Choosing c1=1/8,  c2=-1,  c3=1/4,  c4=-1, it is easy to check that 01G(t,s)f1(c1)ds-Nl1532(1833104+52)-3132(512+π24)153252-(512+16)=7564-712>18=c1,01G(t,s)  g1(c2)ds+ND1532(-13104-4)+3132(516+π60)1532(-4)+616=-32<-1=c2,01G(t,s)g1(c3)ds-ND1532[1433104+4(14)1/5]  -3132(516+π60)15324145-381532434-38=2332>14=c3,01G(t,s)f1(c4)ds+Nl1532(-13104-5)+3132(512+π24)1532(-5)+712=-7532+712<-1=c4. Therefore, condition (H3) is satisfied. At last, we will check condition (H4), by the method of [9, 31], and we consider the linear eigenvalue problem u′′+λu=0,0<t<1,u(0)=0,u(1)=12u(14). Let Γ(s)=coss-(1/2)cos(s/4). By the paper , we know that the sequence of positive eigenvalue of (4.4) is exactly given by λn=sn2,    n=1,2,, where sn is the sequence of positive solutions of Γ(s)=0. In , Han obtained that s1=1.0675, moreover λ1=s12=1.1396. It is easy to know that lim|x|f1(x)+f2(y)x=lim|x|(x3/30000)+5x1/3+(1/6)1005+(1/12)arctan6yx=+21.1396=2λ1, uniformly for y, lim|y|g1(y)+g2(x)y=lim|y|(y3/30000)+4y1/5+(1/16)1003+(1/30)arctan7xy=+21.1396=2λ1, uniformly for x. Therefore, condition (H4) is also satisfied. Consequently, all conditions of Theorem 3.1 are satisfied, and we get the system of differential equations (4.1) has at least six solutions in C[0,1]×C[0,1].

Acknowledgments

The authors are very grateful to the anonymous referees for their valuable suggestions. This project is supported by the National Natural Science Foundation of China (10971046), the University Science and Technology Foundation of Shandong Provincial Education Department (J10LA62), the Natural Science Foundation of Shandong Province (ZR2009AM004, ZR2010AL014), and the Doctor of Scientific Startup Foundation for Shandong University of Finance (08BSJJ32).

MoshinskyM.Sobre los problems de condiciones a la frontiera en una dimension de caracteristicas discontinuasBoletín de la Sociedad Matemática Mexicana195071250043985TimoshenkoS. P.Theory of Elastic Stability1961New York, NY, USAMcGraw-Hill0134026Il’inV. A.MoiseevE. I.onlocal boundary value problem of the first kind for a Sturm-Liouville operator in its differential and finite difference aspectsJournal of Difference Equations198723803810GuptaC. P.Solvability of a three-point nonlinear boundary value problem for a second order ordinary differential equationJournal of Mathematical Analysis and Applications19921682540551117601010.1016/0022-247X(92)90179-HZBL0763.34009GuptaC. P.A generalized multi-point boundary value problem for second order ordinary differential equationsApplied Mathematics and Computation1998891–3133146149169910.1016/S0096-3003(97)81653-0ZBL0910.34032Aykut HamalN.YorukF.Positive solutions of nonlinear m-point boundary value problems on time scalesJournal of Computational and Applied Mathematics2009231192105253265210.1016/j.cam.2009.02.003InfanteG.PietramalaP.Existence and multiplicity of non-negative solutions for systems of perturbed Hammerstein integral equationsNonlinear Analysis2009713-41301131010.1016/j.na.2008.11.0952527550ZBL1169.45001KhanR. A.WebbJ. R. L.Existence of at least three solutions of a second-order three-point boundary value problemNonlinear Analysis20066461356136610.1016/j.na.2005.06.0402200497ZBL1101.34005MaR.Existence of solutions of nonlinear m-point boundary-value problemsJournal of Mathematical Analysis and Applications20012562556567182175710.1006/jmaa.2000.7320MaR.O'ReganD.Solvability of singular second order m-point boundary value problemsJournal of Mathematical Analysis and Applications20053011124134210592510.1016/j.jmaa.2004.07.009XuF.ChenZ.XuF.Multiple positive solutions for nonlinear second-order m-point boundary-value problems with sign changing nonlinearitiesElectronic Journal of Differential Equations20081122392949YangZ.Positive solutions to a system of second-order nonlocal boundary value problemsNonlinear Analysis20056271251126510.1016/j.na.2005.04.0302154107ZBL1089.34022XianX.Three solutions for three-point boundary value problemsNonlinear Analysis20056261053106610.1016/j.na.2005.04.0172152997ZBL1076.34011ZhangZ.WangJ.The upper and lower solution method for a class of singular nonlinear second order three-point boundary value problemsJournal of Computational and Applied Mathematics200214714152193176810.1016/S0377-0427(02)00390-4ZBL1019.34021SunJ.ZhangK.On the number of fixed points of nonlinear operators and applicationsJournal of Systems Science and Complexity20031622292351979577ZBL1131.47308SunJ.CuiY.Multiple solutions for nonlinear operators and applicationsNonlinear Analysis20076691999201510.1016/j.na.2006.02.0382304975ZBL1120.47050XianX.Solutions for an operator equation under the conditions of pairs of paralleled lower and upper solutionsNonlinear Analysis20086972251226610.1016/j.na.2007.08.0052437623ZBL1165.47044ZhangK. M.SunJ. X.Multiple solutions for superlinear operator equations in Banach spaces and applicationsActa Mathematica Sinica2005481991082128657ZBL1117.47304XuX.O'ReganD.SunJ.Multiplicity results for three-point boundary value problems with a non-well-ordered upper and lower solution conditionMathematical and Computer Modelling2007451-2189200227344910.1016/j.mcm.2006.05.003ZBL1140.34009YangZ.Existence and nonexistence results for positive solutions of an integral boundary value problemNonlinear Analysis20066581489151110.1016/j.na.2005.10.0252246350ZBL1104.34017WebbJ. R. L.LanK. Q.Eigenvalue criteria for existence of multiple positive solutions of nonlinear boundary value problems of local and nonlocal typeTopological Methods in Nonlinear Analysis2006271911152236412ZBL1146.34020MaR.O'ReganD.Nodal solutions for second-order m-point boundary value problems with nonlinearities across several eigenvaluesNonlinear Analysis20066471562157710.1016/j.na.2005.07.0072200159RynneB. P.Spectral properties and nodal solutions for second-order, m-point, boundary value problemsNonlinear Analysis200767123318332710.1016/j.na.2006.10.0142350888KongL.KongQ.WongJ. S. W.Nodal solutions of multi-point boundary value problemsNonlinear Analysis201072138238910.1016/j.na.2009.06.0722574948ZBL1195.34028WebbJ. R. L.Remarks on a non-local boundary value problemNonlinear Analysis20107221075107710.1016/j.na.2009.07.0472579370ZBL1186.34029DeimlingK.Nonlinear Functional Analysis1985Berlin, GermanySpringer787404GuoD. J.LakshmikanthamV.Nonlinear Problems in Abstract Cones1988Boston, Mass, USAAcademic Press959889GuoD.Nonlinear Functional Analysis20012ndJinan, ChinaShandong Science and Technology PressSunJ. X.Nontrivial solutions of superlinear Hammerstein integral equations and their applicationsChinese Annals of Mathematics. Series A198675528535886319AmannH.Fixed point equations and nonlinear eigenvalue problems in ordered Banach spacesSIAM Review1976184620709041543210.1137/1018114ZBL0345.47044HanX.-l.Positive solutions of a three-point boundary value problemJournal of Mathematical Research and Exposition20072734975042349495ZBL1144.34311