\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2009(2009), No. 103, pp. 1--13.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2009 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2009/103\hfil Singular $\phi$-Laplacian BVPs]
{Existence and multiplicity results for Singular
$\phi$-Laplacian BVPs on the positive half-line}

\author[S. Djebali, K. Mebarkig\hfil EJDE-2009/103\hfilneg]
{Sma\"il Djebali, Karima Mebarki}  % in alphabetical order

\address{Sma\"il Djebali \newline
Department of Mathematics, Ecole Normale Superieure \\
P.O. Box 92, 16050 Kouba, Algiers, Algeria}
\email{djebali@ens-kouba.dz, djebali@hotmail.com}

\address{Karima Mebarki \newline
Faculty of Fundamental Sciences \\
U.M.B.B., 35000. Boumerd\`es, Algeria}
\email{mebarqi@hotmail.fr}

\thanks{Submitted May 14, 2009. Published August 27, 2009.}
\subjclass[2000]{34B10, 34B15, 34B16, 34B18, 34B40}
\keywords{$\phi$-Laplacian; three point BVPs; positive solutions;
singularity; \hfill\break\indent fixed point in a cone}

\begin{abstract}
 This work  proves the existence and multiplicity of
 positive solutions for a second-order nonlinear three-point
 $\phi$-Laplacian boundary-value problem posed on the positive
 half-line. The nonlinearity depends on the solution and its
 derivative and may exhibit a time singularity at the origin.
 Existence of single and multiple nontrivial positive
 solutions is proved using fixed point index theory.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}{Lemma}[section]
\newtheorem{remark}{Remark}[section]

\section{Introduction}

 This paper concerns the existence of positive
solutions to the following three-point $\phi$-Laplacian boundary
value problem posed on the positive half-line:
\begin{equation}\label{GP}
\begin{gathered}
-(\phi(y'))'(t)=m(t)f(t,y(t),y'(t)),\quad t\in I\\
y(0)=\alpha y'(\eta),\quad\lim_{t\to+\infty}y'(t)=0,
\end{gathered}
\end{equation}
where $\alpha\geq0$ and $\eta\in(0,\infty)$ are given real
numbers. The interval $I:=(0,+\infty)$ denotes the set of positive
real numbers and $\mathbb{R}^+:=[0,+\infty)$. The function
$f:I\times\mathbb{R}^+\times\mathbb{R}\to\mathbb{R}^+$ is
continuous; the function $m:I\to \mathbb{R}^+$ is continuous but
is allowed to have a singularity at $t=0$. $\phi$ is a nonlinear
operator of derivation generalizing the $p$-Laplacian operator.

Boundary value problems on the half-line arise in many
applications of physical phenomena and in chemistry and biology. A
general survey of existence theory is well developed in the books
by Agarwal  et al  \cite{AO, AOW}. In case of second-order
differential equations corresponding to $\phi=I_d$ , Problem
\eqref{GP} has been extensively studied in the literature. Using
the theory of fixed point index on cones of Banach spaces, the
authors obtained in \cite{DMe2, DMe3} some existence results for
the generalized Fisher equation
$-y''+cy'+\lambda y=f(t,y(t),y'(t))$ ($c,\lambda>0$)
subject to Dirichlet or Neumann
limit condition at positive infinity; see also \cite{TianGe1} for
the case of a multi-point condition at the origin. Indeed, since
the pioneer works of Gupta \cite{Gup, GupTro, TianGeShan}, much
attention has been devoted to three-point and more generally to
multi-point boundary value problems (see \cite{GuoGe, LianGe}).
When the derivative operator is generalized to $(q(t)x'(t))'$ and
$f$ depends on the first derivative, the Mawhin coinc\"{\i}dence
theory is applied in \cite{LiuLiFang} to get the existence of at
least one solution.

However, some interesting recent works have been developed for the
case of the so-called $p$-Laplacian operator $\phi_p(s)=|
s|^{p-1}s\;(p>1)$. While the theory is well developed for
$p$-Laplacian problems on bounded intervals (see e.g., \cite{BDM,
BM} and the references therein), less results are known for BVPs
posed on infinite intervals. By means of the three-functional
fixed point theorem, existence of three positive solutions are
obtained in \cite{GuoYuWa} when the nonlinearity depends on the
first derivative and in \cite{LiangZhang} when it does not; some
local growth conditions are assumed on the nonlinearity. Using the
upper and lower solution technique, existence results of single
and double solutions are obtained in \cite{DSa}. The nonlinear
alternative of Leray and Schauder has been recently employed in
\cite{TianGe2} to prove existence of positive solutions for a
multi-point boundary-value problem associated with the equation
$$
\big(\rho(t)| x'(t)^{p-2}| x'(t)\big)'+f(t,x(t),x'(t))=0
$$
on $(0,+\infty)$ when $\rho\in C[0,+\infty)\cap C^1(0,+\infty)$
is positive and satisfies
$$
\int_0^\infty\phi_p^{-1}(1/\rho(t))dt<\infty.
$$

In this work, we assume that the nonlinear map
$\phi:\mathbb{R}\to\mathbb{R}$ is an increasing homeomorphism such
that $\phi(0)=0$ and
\begin{equation}\label{PHI1}
|\phi^{-1}(x)|\le\phi^{-1}(|x|),\quad \forall x\in \mathbb{R}.
\end{equation}
In the second part of this paper, we further assume that $\phi$ is
sub-multiplicative; i.e.,
\begin{equation}\label{PHI2}
\forall\,\alpha,\beta\in \mathbb{R}^+,\quad
\phi(\alpha\beta)\le\phi(\alpha)\phi(\beta).
\end{equation}

Note that if $\phi$ is sub-multiplicative, then the converse
$\phi^{-1}$ is super-multiplicative, that is
\begin{equation}\label{PHI3}
\forall\,\alpha,\beta\in \mathbb{R}^+,\quad
\phi^{-1}(\alpha\beta)\ge\phi^{-1}(\alpha)\phi^{-1}(\beta).
\end{equation}
Clearly, $\phi$ is an extension of the usual p-Laplacian nonlinear
operator which is sub-multiplicative and super-multiplicative,
hence multiplicative.

Finally, recall that various physiological processes are modelled by
singular differential equations. For instance,  the electrical
potential in an isolated neutral atom is governed by the following
problem derived in 1927 by L.H. Thomas \cite{Thomas} and Fermi
\cite{Fermi},
\begin{gather*}
y''=\sqrt{y^3/t}\\
y(0)=1,\quad y(+\infty)=0.
\end{gather*}
A more general survey on singular boundary value problems can be
found in \cite{ORegan}. This is our main motivation of considering
the case when the factor $m$ is time-singular. Our objective is
then to prove some existence results of nontrivial positive
solutions for Problem \eqref{GP} under suitable conditions on the
functions $f$ and $m$. Throughout, by a solution we mean a
positive solution $y\in C^1[0,+\infty)$ such that $\phi(y')\in
C^1(0,+\infty)$ with $y(t)\ge 0$ on $[0,+\infty)$ and the equation
in \eqref{GP} is satisfied. Some preliminaries including the main
assumptions, the problem transformation and a compactness
criterion are gathered together in Section 2. Section 3 is devoted
to proving two existence theorems, one of a single positive
solution and the other one of three positive solutions. An example
of application with a nonlinear operator of derivation ends this
paper in Section 4.

\section{Preliminaries}

\subsection{Functional framework}

In this section, we present some definitions and lemmas we need in
the proofs of the main results. For some real parameter $\theta>0$,
consider the space
$$
X=\big\{y\in C^1([0,\infty),\mathbb{R}):\;
\lim_{t\to+\infty}\frac{y(t)}{e^{\theta t}}\, \text{ exists and
}\, \lim_{t\to+\infty}y'(t)=0\big\}
$$
with the norm
$$
\|y\|_{\theta}=\max\big\{\Vert y\Vert_1,\Vert y\Vert_2\big\},
$$
where
$$
\Vert y\Vert_1=\sup_{t\in[0,\infty)}\frac{|y(t)|}{e^{\theta
t}},\quad \Vert y\Vert_2=\sup_{t\in [0,\infty)}|y'(t)|\,.
$$

\begin{remark}\label{rem1} \rm
Clearly $X$ is a Banach space. Moreover, if $y\in X$ is such that $y(0)=\alpha
y'(\eta)$, then
\begin{align*}
\frac{y(t)}{e^{\theta t}}
&= e^{-\theta t}\big\{\int_0^ty'(s)\,ds+y(0)\big\}\\
&= e^{-\theta t}\big\{\int_0^ty'(s)\,ds+\alpha y'(\eta)\big\}\\
&\leq e^{-\theta t}\big\{t\|y\|_2+\alpha \|y\|_2\big\}\\
&= \frac{t+\alpha}{e^{\theta t}}\|y\|_2,\quad\forall\,
t\in\mathbb{R}^+.
\end{align*}
Hence $\|y\|_1\le K\|y\|_2$ where
$K=\sup_{t\in\mathbb{R}^+}\gamma(t)$ with
$\gamma(t)=\frac{t+\alpha}{e^{\theta t}}$; more
precisely,
$$
K=\begin{cases}
\alpha, & \text{if } \theta\alpha \ge1\\
\gamma\big(\frac{1-\theta\alpha}{\theta}\big), & \text{if
} 0\le\theta\alpha<1.
\end{cases}
$$
As a consequence
$$
\|y\|_2\le\|y\|_{\theta}\le \max\{1,K\}\|y\|_2
$$
and
$\lim_{t\to+\infty}\frac{y(t)}{e^{\theta t}}=0$.
\end{remark}

\subsection{Integral formulation}

In order to transform \eqref{GP} into a fixed point problem,
we need the following auxiliary lemma the proof of which is omitted.

\begin{lemma}\label{lemm1}
Let $v\in L^1(I)$. Then $y\in C^1(I)$ is a solution of
\begin{equation}\label{GP1}
\begin{gathered}
-(\phi(y'))'(t)=v(t),\quad t\in I\\
y(0)=\alpha y'(\eta),\quad\lim_{t\to+\infty}y'(t)=0
\end{gathered}
\end{equation}
if and only if
\begin{equation}\label{eq1}
y(t)=C+\int_0^t\phi^{-1}\Big(
\int_s^{+\infty}v(\tau)\,d\tau\Big)ds,\quad t\in I,
\end{equation}
where $C=\alpha \phi^{-1}\Big(\int_\eta^{+\infty}v(\tau)\,d\tau\Big)$.
\end{lemma}

\subsection{General assumptions and a fixed point operator}

Assume first that the nonlinearity satisfies the following
hypotheses:
\begin{itemize}
\item[(H1)] The function
$f:I\times\mathbb{R}^+\times\mathbb{R}\to\mathbb{R}^+$ is
continuous and when $y,z$ are bounded, $f(t,e^{\theta t}y,z)$ is
bounded on $[0,+\infty)$.

\item[(H2)] The function $m:I\to\mathbb\mathbb{R}^+$ is continuous
and does not vanish identically on any subinterval of $I$. It may
be singular at $t=0$ but satisfies
\begin{equation}\label{integralm}
A:=\int_0^{+\infty}m(s)ds<\infty.
\end{equation}
\end{itemize}
The first hypothesis means that $f$ is bounded in term of the
variable $t$ and is justified by the fact the unboundedness of the
nonlinearity is carried by the singular factor $m$. The second
hypothesis means that the singularity is integrably bounded.

For a bounded subset $\Omega\subset X$, define the integral
operator
\begin{equation}\label{operator}
Fy(t)=C+\int_0^t\phi^{-1}\Big(
\int_s^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau\Big)ds,\quad
t\in I,
\end{equation}
where $C=\alpha \phi^{-1}\big(\int_\eta^{+\infty}m(\tau)
f(\tau,y(\tau),y'(\tau))d\tau\big)$.
By Lemma \ref{lemm1}, all solutions of  \eqref{GP} are
fixed points of $F$ on $X$ and conversely. Define the
cone
$$
P:=\big\{y\in X: y(t)\ge 0 \text{ on } \mathbb{R}^+\text{ with
}\;y(0)=\alpha y'(\eta)\big\}.
$$
We have

\begin{lemma}\label{lemma1}
Under Assumptions {\rm (H1), (H2)}, $F$ maps the set
$\overline{\Omega}\cap P$ into $P$.
\end{lemma}

\begin{proof}
First we show that $F:\overline{\Omega}\cap P\to X$ is well
defined. Let $y\in\overline{\Omega}\cap P$. Then, there exists
$M>0$ such that $\|y\|_{\theta}\le M$. By Assumption (H1), let
$$
S_M=\sup\{f(t,e^{\theta t}y,z),\; t\in I, (y,| z|)\in [0,M]^2\}.
$$
Since, for any $t\ge0$, $0\le y(t)e^{-\theta t}\le M$ and
$| y'(t)|\le M$,
Assumption (H2) implies
$$
\int_0^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau=
\int_0^{+\infty}m(\tau)f\big(\tau,e^{\theta
\tau}y(\tau)e^{-\theta\tau},y'(\tau)\big)d\tau\le AS_M.
$$
Hence for any fixed $t\in(0,+\infty)$,
$$
\int_0^t\phi^{-1}\Big(
\int_s^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau\Big)ds\le
\int_0^t\phi^{-1}(AS_M)ds<\infty.
$$
In addition, we can easily prove that for each
$y\in\overline{\Omega}\cap P$,
$$
Fy\in C^1([0,\infty),\mathbb{R}),\; Fy(t)\ge0,\;t\in I,$$
$$
Fy(0)=C=\alpha \phi^{-1}\Big(
\int_\eta^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau\Big)=\alpha
(Fy)'(\eta),
$$
and
$$
\lim_{t\to+\infty}(Fy)'(t)=\lim_{t\to +\infty}\phi^{-1}\Big(
\int_t^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau\Big)=\phi^{-1}(0)=0.
$$
By Remark \ref{rem1}, we obtain
\begin{align*}
0&\le\lim_{t\to+\infty}\frac{Fy(t)}{e^{\theta
t}} \le \lim_{t\to+\infty}\gamma(t)\sup_{t\in
[0,\infty)}|(Fy)'(t)|\\
&\le \lim_{t\to+\infty}\gamma(t)\phi^{-1}(AS_M)=0.
\end{align*}
\end{proof}

\subsection{A compactness criterion}

To investigate the compactness of the operator $F$, we
appeal to the following result.

\begin{lemma}[\cite{Cord}] \label{lem2}
Let $M\subseteq C^1_{\infty}(\mathbb{R}^+,\mathbb{R})$. Then the set
$M$ is relatively compact in $C^1_{\infty}(\mathbb{R}^+,\mathbb{R})$
if the following conditions hold:
\begin{itemize}
\item[(a)] $M$ is uniformly bounded in
$C^1_{\infty}(\mathbb{R}^+,\mathbb{R})$.
\item[(b)] The functions belonging to the sets
$$
A=\{y: y(t)=\frac{x(t)}{e^{\theta t}},\; x \in M\} \text{ and }
B=\{z: z(t)=x'(t),\; x\in M\}
$$
are almost equi-continuous on $\mathbb{R}^+$.
\item[(c)] The functions from $A$ and $B$
are equi-convergent at $+\infty$.
\end{itemize}
\end{lemma}

Next, we state two compactness results depending on whether or not
the function $m$ is singular at the origin.

\begin{lemma}[The regular case]\label{lemma2}
Assume $m:[0,\infty)\to [0,\infty)$ is continuous. Then, the
mapping $F:\overline{\Omega}\cap P\to P$ is completely continuous.
\end{lemma}

\begin{proof}
\textbf{Claim 1.} $F$ is continuous on $P$. Assume
$\lim_{n\to+\infty}y_n\to y$ in $P;$ then there exists $N>0$
independent of $n$ such that
$\max\{\|y\|_{\theta},\sup_{n\ge1}\|y_n\|_{\theta}\}\le N$.
Letting
$$S_N=\sup\{f(t,e^{\theta t}y,z),\; t\in [0,+\infty),
(y,| z|)\in [0,N]^2\},
$$
we get
$$
\int_0^{+\infty}m(s)(f(s,y_n(s),y_n'(s))-f(s,y(s),y'(s)))ds\le
2AS_N.
$$
Then, the Lebesgue's dominated convergence theorem both with the
continuity of $f$ and $\phi^{-1}$ imply
\begin{align*}
&|(Fy_n)'(t)-(Fy)'(t))|\\
&=\big|\phi^{-1}\Big(
\int_t^{+\infty}m(\tau)f(\tau,y_n(\tau),y_n'(\tau))d\tau\Big)
-\phi^{-1}\Big(\int_t^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau
\Big)\big|\\
&\to 0, \quad \text{as } n\to +\infty.
\end{align*}
Consequently,
$$
\|Fy_n-Fy\|_{\theta}\le
\max\{1,K\}\|Fy_n-Fy\|_2\to 0,\quad \text{ as } n\to +\infty,
$$
yielding our claim.

\textbf{Claim 2.} $F$ is compact provided it maps bounded sets
into relatively compact sets. Let $\Omega$ be any bounded subset
of $X;$ then there exists $M>0$ such that $\|y\|_{\theta}\le M$
for all $y\in\overline{\Omega}\cap P$. On one hand
$$
\|Fy\|_2=\phi^{-1}\Big(
\int_0^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))\,d\tau\Big)\le
\phi^{-1}(AS_M),\; \forall\,y\in\overline{\Omega}\cap P.
$$
Hence $\|Fy\|_{\theta}\le \max\{1,K\}\|Fy\|_2\le
\max\{1,K\}\phi^{-1}(AS_M)$ which proves that
$F(\overline{\Omega}\cap P)$ is uniformly bounded. On the other
hand, for any $y\in\overline{\Omega}\cap P$, any
$T\in(0,+\infty)$ and $t_1,t_2\in[0,T]$, we have
\begin{align*}
& |\frac{Fy(t_2)}{e^{\theta t_2}}-\frac{Fy(t_1)}{e^{\theta
t_1}}|\\
&= | C(e^{-\theta t_2}-e^{-\theta t_1})
+e^{-\theta t_2}\int_{0}^{t_2}\phi^{-1}
\Big(\int_s^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau\Big)ds
\\
&\quad -e^{-\theta t_1}\int_{0}^{t_1}\phi^{-1}
\Big(\int_s^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau\Big)ds|\\
&\le \Big(C+\int_{0}^{t_2}\phi^{-1}(AS_M)ds\Big)
| e^{-\theta t_2}-e^{-\theta t_1}|+e^{-\theta t_1}
|\int_{t_2}^{t_1}\phi^{-1}(AS_M)ds|
\end{align*}
which tends to $0$ as $| t_1-t_2|\to 0$. Also, by the
continuity of $\phi^{-1}$,
\begin{align*}
&|(Fy)'(t_2)-(Fy)'(t_1)|\\
&= \Big|\phi^{-1}\Big(\int_{t_2}^{+\infty}m(s)f(s,y(s),y'(s))ds\Big)
-\phi^{-1}\Big(\int_{t_1}^{+\infty}m(s)f(s,y(s),y'(s))ds\Big)\Big|\\
&= \Big|\phi^{-1}\Big(\int_{t_2}^{t_1}m(s)f(s,y(s),y'(s))ds
+\int_{t_1}^{+\infty}m(s)f(s,y(s),y'(s))ds\Big)\\
&\quad -\phi^{-1}\Big(\int_{t_1}^{+\infty}m(s)f(s,y(s),y'(s))ds\Big)\Big|
\end{align*}
tends to $0$ as $| t_1-t_2|\to 0$. This proves that
$F(\overline{\Omega}\cap P)$ is equi-continuous. Since (H2)
yields
$$
\lim_{t\to+\infty}\int_t^{+\infty}m(s)\,ds=0,
$$
the Lebesgue dominated convergence theorem implies
\begin{align*}
\lim_{t\to+\infty}|\frac{Fy(t)}{e^{\theta
t}}-\lim_{s\to+\infty}\frac{Fy(s)}{e^{\theta s}}|
&\leq \lim_{t\to \infty}\gamma(t) \sup_{t\in
[0,\infty)}|(Fy)'(t)|\\
&\leq \lim_{t\to \infty}
\gamma(t)\phi^{-1}\Big(S_M\int_0^{+\infty}m(s)\,ds\Big)=0
\end{align*}
and
\begin{align*}
\lim_{t\to+\infty}|(Fy)'(t)-\lim_{s\to+\infty}(Fy)'(s)|
&=\lim_{t\to+\infty}|
\phi^{-1}\Big(\int_t^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))
\,d\tau\Big)|\\
&\leq \phi^{-1}\Big( S_M
\lim_{t\to+\infty}\int_t^{+\infty}m(s)\,ds\Big)=0.
\end{align*}
This means that $F(\overline{\Omega}\cap P)$ is equi-convergent at
$\infty$. By Lemma \ref{lem2}, $F(\overline{\Omega}\cap P)$ is
relatively compact.
\end{proof}

\begin{lemma}[The singular case]\label{lemma3}
Let $m$ be singular at $t=0$. Then, the mapping $F$ given by
\eqref{operator} is completely continuous.
\end{lemma}

\begin{proof} For each $n\geq1$, define the approximating operator
$F_n$ on $\overline{\Omega}\cap P$ by
$$
F_ny(t)=C+\int_{\frac{1}{n}}^t\phi^{-1}\Big(
\int_s^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau\Big)ds,\quad
t\in I.
$$
Thus it suffices to prove that $F_n$ converges uniformly to $F$ on
$\overline{\Omega}\cap P$. For any $ t\in I$ and
$y\in\overline{\Omega}\cap P$ satisfying $\|y\|_{\theta}\le M$, by
(H1) the following estimates hold
\begin{align*}
|F_ny(t)-Fy(t))|e^{-\theta t}
&= \big|\int_0^{\frac{1}{n}}e^{-\theta t}\phi^{-1}\Big(
\int_s^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))d\tau\Big)ds\big|\\
&\leq \frac{1}{n}e^{-\theta t}\phi^{-1}\big(AS_{M}\big),
\end{align*}
and
$$
|(F_ny)'(t)-(Fy)'(t))|=0,\quad\forall\,n\in\mathbb{N}.
$$
Consequently, Assumption (H2) both with the Cauchy criterion for
convergent integrals imply that
$$
\|F_ny-Fy\|_{\theta}=
\max\{\|F_ny-Fy\|_1,\|F_ny-Fy\|_2\}\to 0,\quad \text{as }n\to +\infty.
$$
Since from Lemma \ref{lemma2}, the operator
$F_n:\overline{\Omega}\cap P\to P$ is completely continuous for
each $n\geq1$ and $F_n$ converges to $F$ uniformly on any closed,
bounded subset of $\overline{\Omega}\cap P$, the uniform limit
operator $F$ is completely continuous, proving the lemma.
\end{proof}

\section{Existence results}

\subsection{One positive solution}

The following Lemma is needed in this section.
Detailed properties of the fixed point index on cones of
Banach spaces may be found in \cite{Deim, Zeid}.

\begin{lemma}[\cite{AMO, Deim, GuoLak, Zeid}] \label{lemA}
Let $\Omega$ be a bounded open subset of a real Banach space
$E,P$ a cone of $E,\,\theta\in \Omega$ and
$A: \overline{\Omega}\cap P\to P$ a completely continuous operator.
Suppose that
$$
Ax\neq\lambda x,\quad\forall\,x\in\partial\Omega\cap
P,\;\lambda\ge1.
$$
Then the index $i(A,\Omega\cap P,P)=1$.
\end{lemma}

The main existence result in this section is

\begin{theorem}\label{th1}
Assume  {\rm (H1), (H2)} and
\begin{itemize}
\item[(H3)] for all $(t,y,z)\in I\times\mathbb{R}^+\times\mathbb{R}$,
\[
0\le f(t,y,z)\le a(t)\phi(e^{-\theta
t}y)+b(t)\phi(| z|)+c(t),
\]
where
$a, b, c\in C^0(\mathbb{R}^+)$,  $mb,mc\in L^1(I)$,  and
there exists $R>0$, such that
\begin{equation}\label{Rhypothesis}
\phi^{-1}\left((|ma|_{L_1}+|mb|_{L_1})\phi(R)+|mc|_{L_1}\right)<
\frac{R}{\max\{K,1\}}\,.
\end{equation}
\end{itemize}
Then  \eqref{GP} has at least one nonnegative, concave, and
nondecreasing solution $y\in P\cap{B}_\theta(0,R)$ where
${B}_\theta(0,R)$ is the open ball centered at the origin with
radius $R$ in the $\theta-$weighted space $X$.

If further $\min_{I\times[0,R]\times[0,R]}f(t,e^{\theta t}y,z)\ge1$,
then
$y(t)\ge w(t)$ for all $t\in I$,
where
\begin{equation}\label{omega}
\omega(t):=\alpha
\phi^{-1}\Big(\int_\eta^{+\infty}m(\tau)d\tau\Big)
+\int_0^t\phi^{-1}\Big(\int_s^{+\infty}m(\tau)d\tau\Big)\,ds,\quad
t\in I.
\end{equation}
\end{theorem}

\begin{remark} \label{rm3.1} \rm
The properties of $\phi$ and $m$ imply that $\omega\ge0$ and
$\omega$ is well defined. Moreover $\omega$ is the unique solution
of Problem (\ref{GP1}) for $v\equiv m$.
\end{remark}

\begin{proof}
Consider the open ball
$\Omega:=\{y\in X:\|y\|_\theta<R\}$.
 From Lemmas \ref{lemma2} and \ref{lemma3}, the mapping
$ F: \overline{\Omega}\cap P\to P$ is completely
continuous.

\textbf{Claim 1.} $Fy\neq\lambda y$, for any
$y\in\partial\Omega\cap P$  and $\lambda\geq1$.

Let $y\in\partial\Omega\cap P$. By Assumption (H3), the
following estimates hold
\begin{align*}
|(Fy)'(t)|
&=  \phi^{-1}\Big(
\int_t^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))\,d\tau\Big)
\\
&\leq \phi^{-1}\Big(
\int_0^{+\infty}m(\tau)\left(a(\tau)\phi(e^{-\theta
\tau}y(\tau))+b(\tau)\phi(|y'(\tau)|)+c(\tau)\right)\,d\tau\Big)
\\
&\leq \phi^{-1}\Big(
\left(|ma|_{L_1}\phi(\|y\|_1)+|mb|_{L_1}\phi(\|y\|_2)\right)+|mc|_{L_1}\Big)
\\
&\leq \phi^{-1}\left(
\left(|ma|_{L_1}+|mb|_{L_1}\right)\phi(\|y\|_{\theta})+|mc|_{L_1}\right)
\\
&\leq \phi^{-1}\left(
\left(|ma|_{L_1}+|mb|_{L_1}\right)\phi(R)+|mc|_{L_1}\right)
\\
&< \frac{R}{\max\{1,K\}}=\frac{\Vert y\Vert_\theta}{\max\{1,K\}}\,.
\end{align*}
Passing to the supremum over $t$, we infer that
$$
\|Fy\|_2<\frac{1}{\max\{1,K\}}\|y\|_{\theta},\quad
\forall\,y\in\partial\Omega\cap P.
$$
Hence,
\begin{equation}\label{EQ2é}
\|Fy\|_\theta\le\max\{1,K\}\|Fy\|_2<\|y\|_\theta,\quad
\forall\,y\in\partial\Omega\cap P.
\end{equation}
As a consequence
\begin{equation}\label{eq12}
Fy\neq\lambda y,\quad \forall\,y\in\partial\Omega\cap
P,\;\forall\,\lambda\geq1.
\end{equation}
Indeed, on the contrary there would exist some
$y_0\in \partial\Omega\cap P$ and $\lambda_0\geq1$ such that
$Fy_0=\lambda_0y_0$. Thus
$$
\|Fy_0\|_\theta=\lambda_0\|y_0\|_\theta\geq\|y_0\|_\theta=R,
$$
contradicting (\ref{EQ2é}). This implies that (\ref{eq12}) holds.
Therefore, Lemma \ref{lemA} yields
\begin{equation}\label{eq13}
i\,(F,\Omega\cap P,P)=1.
\end{equation}
Hence (\ref{eq13}) and the solution property of the fixed point
index imply that the operator $F$ has a fixed point $y$ which
belongs to ${\Omega}\cap P$. Moreover, we have that
$(\phi(y'))'(t)=(\phi((Fy)'))'(t)=-m(t)f(t,y(t),y'(t))\le0,\, t\in
I$, which implies that $y$ is concave on $I$. Next, we show that
$y$ is a nontrivial solution.

\textbf{Claim 2.} For this fixed point $y$, we claim
that $Fy(t)\ge\omega(t)$ on $I$, where $\omega$ is as given by
(\ref{omega}).
Otherwise, we have
$$
\sup_{t\in I}\{\omega(t)-Fy(t)\}>0.
$$
Now, we distinguish between two cases.

\textbf{Case 1.}
$\lim_{t\to+\infty}\{\omega(t)-Fy(t)\}
=\sup_{t\in I}\{\omega(t)-Fy(t)\}>0$.
Under the assumption that
$\min_{I\times[0,R]\times[0,R]}f(t,e^{\theta t}y,z)\ge1$ and using
the fact that $\phi^{-1}$ is nondecreasing, we get
\begin{align*}
\lim_{t\to+\infty}\{\omega(t)-Fy(t)\}
&=  \alpha \phi^{-1}\Big(\int_\eta^{+\infty}m(\tau)\,d\tau\Big)
+\int_0^{+\infty}\phi^{-1}\Big(\int_s^{+\infty}m(\tau)\,d\tau\Big)ds\\
&\quad -\alpha \phi^{-1}\Big(
\int_\eta^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))\,d\tau\Big)\\
&\quad -\int_0^{+\infty}\phi^{-1}\Big(
\int_s^{+\infty}m(\tau)f(\tau,y(\tau),y'(\tau))\,d\tau\Big)ds
\leq  0.
\end{align*}
This is a contradiction to the assumption in Case 1.

\textbf{Case 2.} There exists a real number $t_1\ge 0$ such that
$$
\omega(t_1)-Fy(t_1) =\sup_{t\in\mathbb{R}^+}\{\omega(t)-Fy(t)\}>0.
$$
Arguing as in case 1, we can easily check that
$\omega(t_1)-Fy(t_1)\le0$ and a contradiction is reached.
Therefore $y$ is a nonnegative, concave, and nondecreasing
solution to Problem \eqref{GP} and satisfies
$$
0\le\|y\|_{\theta}< R,\quad
y(t)\ge \omega(t),\quad \forall\,t\in I.
$$
\end{proof}

\subsection{A multiplicity result}

The following Lemma is needed in this section.

\begin{lemma}[\cite{AMO, Deim, GuoLak}] \label{lemB}
Let $\Omega$ be a bounded open set in a real Banach space $E$, $P$
a cone of $E,\theta\in \Omega$ and
$A: \overline{\Omega}\cap P\to P$ a completely continuous mapping.
Assume that
$$
Ax\not\leq x, \quad \forall\,x\in\partial\Omega\cap P.
$$
Then the index $i(A,\Omega\cap P,P)=0$.
\end{lemma}

\begin{theorem}\label{th2}
Assume {\rm (H1)--(H3)}, and, instead of (\ref{Rhypothesis}),
there exist two constants $0<R_1<R_2$ such that
\begin{equation}\label{Rhypothesis2}
\phi^{-1}\left((|ma|_{L_1}+|mb|_{L_1})\phi(R_i)+|mc|_{L_1}\right)<
\frac{R_i}{\max\{1,K\}},\quad \text{for }i=1,2.
\end{equation}
\begin{itemize}
\item[(H4)]  Assume that for some $0<\gamma<\delta$,
$$
f(t,e^{\theta t}y,z)\ge g(t,y),\quad
\forall\,(t,y,z)\in[\gamma,\delta]\times\mathbb{R}^+\times\mathbb{R},
$$
where $g\in C([\gamma,\delta]\times\mathbb{R}^+)$ and
there exists  $\eta>R_1$ such that
\begin{equation}\label{eq14}
g(t,y)>\frac{\phi(1/\gamma)}{\int_\gamma^\delta
m(\tau)\,d\tau}\,\phi(e^{\theta t}y),\quad\text{for }
y\in[0,\eta],\; t\in[\gamma,\delta].
\end{equation}
\end{itemize}
Then for each constant $R_1<r<\min(R_2,\eta)$, Problem \eqref{GP}
has at least three nonnegative, concave, and nondecreasing
solutions $y_1,y_2,y_3\in P$ satisfying
$0\le\|y_1\|_{\theta}<R_1<\|y_2\|_{\theta}<r<\|y_3\|_{\theta}<R_2.
$ If further $\min_{I\times[0,R_1]\times[0,R_1]}f(t,e^{\theta t}y,z)\ge1$,
then
$$
y_1(t)\ge w(t),\quad \forall\, t\in I.
$$
\end{theorem}

\begin{proof}
\textbf{Claim 1.} Consider the open balls
$\Omega_{R_i}:=\{y\in X:\|y\|_\theta<R_i\}$, $i=1,2$.
Arguing as in Claim 1 of Theorem \ref{th1} and using
(\ref{Rhypothesis2}),
we can check that
\begin{equation}\label{eq2S}
i(F,\Omega_{R_i}\cap P,P)=1,\quad  i=1,2.
\end{equation}

\textbf{Claim 2.} Let
$\ell:=\phi(1/\gamma)/\int_\gamma^\delta m(\tau)\,d\tau$,
$R_1<r<\min(R_2,\eta)$ and consider the open
ball $ \Omega_r:=\{y\in X:\|y\|_\theta<r\}$.
We claim that $Fy\not\leq  y$,  for any $y\in\partial\Omega_r\cap P$.
Otherwise, let $y_0\in\partial\Omega_r\cap P$ be such that
\begin{equation}\label{eq15}
Fy_0\leq y_0.
\end{equation}
Then $0\le e^{-\theta t}y_0(t)\le r<\eta,\;\forall\,t\in[\gamma,\delta]$.
Moreover, by virtue of (H4), (\ref{eq14}) and (\ref{eq15}) together with the
property (\ref{PHI3}) of $\phi^{-1}$ and the definition of $\ell$,
we obtain successively the following estimates:
for every $t\in[\gamma,\delta]$,
\begin{align*}
y_0(t)&\ge
C+\int_0^t\phi^{-1}\Big(
\int_s^{+\infty}m(\tau)f(\tau,y_0(\tau),y_0'(\tau))d\tau\Big)ds\\
&\geq \int_0^\gamma\phi^{-1}\Big(
\int_s^{+\infty}m(\tau)f(\tau,y_0(\tau),y_0'(\tau))d\tau\Big)ds\\
&\geq \int_0^\gamma\phi^{-1}\Big(
\int_\gamma^{\delta}m(\tau)f(\tau,y_0(\tau),y_0'(\tau))d\tau\Big)ds\\
&\geq \int_0^\gamma\phi^{-1}\Big(
\int_\gamma^\delta m(\tau)g(\tau,e^{-\theta\tau}y_0(\tau))d\tau\Big)ds\\
&> \gamma\,\phi^{-1}\Big(
\int_\gamma^\delta m(\tau)\ell\phi(y_0(\tau))\,d\tau\Big)\\
&\geq \gamma\,\phi^{-1}\Big(
\phi(\min_{t\in[\gamma,\delta]}y_0(t))\Big)\phi^{-1}
\Big(\ell\int_\gamma^\delta m(\tau)\,d\tau\Big)\\
&\geq \gamma\,\phi^{-1}\Big(\ell\int_\gamma^\delta m(\tau)\,d\tau\Big)
\min_{t\in[\gamma,\delta]}y_0(t)\\
&\geq \min_{t\in[\gamma,\delta]}y_0(t).
\end{align*}
Hence for any $t\in[\gamma,\delta]$,
$y_0(t)>\min_{t\in[\gamma,\delta]}y_0(t)$, contradicting the
continuity of the function $y_0$ on the compact interval
$[\gamma,\delta]$. This implies that (\ref{eq12}) holds. As a
consequence, Lemma \ref{lemB} yields
\begin{equation}\label{eq17}
i\,(F,\Omega_{r}\cap P,P)=0.
\end{equation}
To sum up, from (\ref{eq2S}), (\ref{eq17}) and the fact that
$\overline{\Omega}_{R_1}\subset
\Omega_{r},\overline\Omega_{r}\subset{\Omega}_{R_2}$, we deduce that
$i(F,(\Omega_{r}\setminus\overline{\Omega}_{R_1})\cap
P,P)=-1$  and
$i(F,(\Omega_{R_2}\setminus\overline{\Omega}_{r})\cap P,P)=1$.
Therefore, there exist three fixed points $y_1,y_2,y_3\in P$
satisfying
$0\le\|y_1\|_{\theta}<R_1<\|y_2\|_{\theta}<r<\|y_3\|_{\theta}<R_2$.
In addition, if
$$
\min_{I\times[0,R_1]\times[0,R_1]}f(t,e^{\theta t}y,z)\ge1,
$$
then we can check as in Theorem \ref{th1} that
$y_1(t)\ge w(t)$ for all $ t\in I$.
\end{proof}

\begin{remark} \label{rmk3.2} \rm
Notice that at least the two solutions $y_2$ and $y_3$ are positive,
whence nontrivial and that, due to the range of values the constant
$r$ may take,
we can obtain as much pairs of solutions $y_2, y_3$ as we need.
\end{remark}

\section{Example}

Consider the increasing homeomorphism defined by
$$
\phi(x)=\begin{cases}
(\frac{1}{8}\times10^{-2}x^2)+(\frac{1}{4}\times10^{-2}),
&x\ge 1; \\
\frac{3}{8}\times10^{-2}x, &  x\le 1.
\end{cases}
$$
Let $a(t)=b(t)=e^{-kt}$ ($k>0$), $c(t)=2$, and
$$
m(t)=\begin{cases}
\frac{1}{t^{10^{-2}}}, &0<t\le 1; \\
\frac{1}{t^{10^{2}}}, & t\ge 1.
\end{cases}
$$
To check the inequality (\ref{Rhypothesis}) in Assumption
(H3), we take $\alpha=1/2$, $\eta=1$, $\gamma=1/5$ and
$\delta=1/3;$ thus we can choose $\theta=1$, $k=30$, and $R=50$.
Moreover
$$
\phi^{-1}(x)=\begin{cases}
\sqrt{8\times10^{2}x-2}, & x\ge\frac{3}{8}\times10^{-2} \\
\frac{8}{3}\times10^2x, & x\le\frac{3}{8}\times10^{-2},
\end{cases}
$$
$\max\{1,K\}= \max\{1,1/\sqrt{e}\}=1$, and
$$
\phi^{-1}\left((|ma|_{L_1}+|mb|_{L_1})\phi(R)+|mc|_{L_1}\right)
=42.4724<50.
$$
Therefore, Assumptions (H1)-(H3) are satisfied. As a consequence,
the singular boundary value problem
\begin{equation}\label{example}
\begin{gathered}
-(\phi(y'))'(t)=m(t)f(t,y(t),y'(t)),\quad t\in I \\
y(0)=\frac{1}{2}y'(1),\quad\lim_{t\to+\infty}y'(t)=0,
\end{gathered}
\end{equation}
where $f(t,y,z)=a(t)\phi(e^{-\theta t}y)+b(t)\phi(| z|)+c(t)$,
$(t,y,z)\in I\times\mathbb{R}^+\times\mathbb{R}$ has at least one
nonnegative, concave,
and nondecreasing solution $y$. Moreover
$$
f(t,e^{\theta t}y,z)\ge1,\quad
\forall\,(t,y,z)\in I\times\mathbb{R}^+\times\mathbb{R}.
$$
Hence
$y(t)\ge w(t)$,  for $t\in I$ where
$w(t):=1.2330+\int_0^t\phi^{-1}(\int_s^{+\infty}m(\tau)\,d\tau)ds$.

To check the inequality (\ref{Rhypothesis2}) in Theorem \ref{th2},
we take $R_1=42$ and $R_2=50$ to get the numerical values
\begin{gather*}
\phi^{-1}\left((|ma|_{L_1}+|mb|_{L_1})\phi(R_1)+|mc|_{L_1}\right)
=41.8670<42.
\\
\phi^{-1}\left((|ma|_{L_1}+|mb|_{L_1})\phi(R_2)+|mc|_{L_1}\right)
=42.4724<50.
\end{gather*}
Moreover,
$$
f(t,y(t),y'(t))\ge g(t,e^{-\theta t}y),\quad \forall\,(t,y,z)\in
I\times\mathbb{R}^+\times\mathbb{R}
$$
with $g(t,e^{-\theta t}y)=a(t)\phi(e^{-\theta t}y(t))+c(t)$. If we
take $\eta=45>R_1$, then for any $y\le\eta$ and
$t\in[\gamma,\delta]$, we have the estimates
$$
\frac{g(t,e^{-\theta t}y)}{\phi(y)}\ge\frac{c(t)}{\phi(y)}\ge\frac{2}{\phi(\eta)}
=0.7893>\frac{\phi(\frac{1}{\gamma})}{\int_\gamma^\delta
m(\tau)\,d\tau}=0.2498.
$$
Therefore,  (H4) in Theorem \ref{th2} is satisfied.
All the computations have been done using Matlab 7. As a
consequence,  for any $r\in(42,45)$, the singular boundary-value
problem (\ref{example})
has at least three nonnegative, concave, and nondecreasing solutions
$y_1,y_2,y_3\in P$ satisfying
$$
0\le\|y_1\|_{\theta}<42<\|y_2\|_{\theta}<r<\|y_3\|_{\theta}<50.
$$

\subsection*{Acknowledgments}
The authors are thankful to the anonymous referees for their careful
reading of the original manuscript, which led to substantial
improvements.

\begin{thebibliography}{00}

\bibitem{AMO} {R. P. Agarwal, M. Meehan and D. O'Regan,} \textit{Fixed Point
Theory and Applications,} Cambridge Tracts in Mathematics, Vol.
{\bf141}, Cambridge University Press, 2001.

\bibitem{AO} {R. P. Agarwal, D. O'Regan,}
\textit{Infinite Interval Problems for Differential, Difference and
Integral Equations,} Kluwer Academic Publisher Dordrecht, 2001.

\bibitem{AOW} {R. P. Agarwal, D. O'Regan, and P.J.Y. Wong,}
\textit{Positive Solutions of Differential, Difference and Integral
Equations,} Kluwer Academic Publisher, Dordrecht, 1999.

\bibitem{BDM} {A. Benmezai, S. Djebali and T. Moussaoui,} \textit{Existence results for one-dimensional
Dirichlet $\phi-$Laplacian BVPs: a fixed point approach,} Dynamic
Syst. Appl. {\bf17} (2008) 149--166.

\bibitem{BM} {C. Bereanu, J. Mawhin,} \textit{Existence and multiplicity results for some nonlinear problems with
singular $\phi-$laplacian,} J. Diff. Equations {\bf243(2)} (2007)
536--557.

\bibitem{Cord} {C. Corduneanu,} \textit{Integral Equations and Stability of
Feedback Systems,} Academic Press, New York, 1973.

\bibitem{Deim} {K. Deimling,} \textit{Nonlinear Functional Analysis,}
Springer-Verlag, Berlin, Heidelberg, 1985.

\bibitem{DMe2} {S. Djebali, K. Mebarki,}
\textit{Multiple positive solutions for singular BVPs on the
positive half-line,} Comp. Math. Appl., {\bf55(12)} (2008)
2940--2952

\bibitem{DMe3} {S. Djebali, K. Mebarki,} \textit{
On the singular generalized Fisher-like equation with derivative
depending nonlinearity,} Appl. Math Comput., {\bf205(1)} (2008)
336--351

\bibitem{DSa} {S. Djebali, O. Saifi,} \textit{Positive solutions for singular $\phi-$Laplacian BVPs
on the positive half-line,} submitted.

\bibitem{Fermi} {E. Fermi,}
\textit{Un methodo statistico par la determinazione di alcune
propriet\'a dell'atome,} Rend. Accad. Naz. del Lincei. CL. sci.
fis., mat. e nat., {\bf6}(1927) 602--607.

\bibitem{GuoGe} {Y. Guo, W. Ge,}
\textit{Positive solutions for three-point boundary value problems
with dependence on the first order derivative,} J. Math. Anal. Appl.
{\bf290} (2004) 291--301.

\bibitem{GuoYuWa} {Y. Guo, C. Yu, and J. Wang,}
\textit{Existence of three positive solutions for $m-$point boundary
value problems on infinite intervals,} Nonlin. Anal., {\bf71} (2009)
717--722.

\bibitem{GuoLak} {D. Guo, V. Lakshmikantham,}
\textit{Nonlinear Problems in Abstract Cones}, Academic Press, San
Diego, 1988.

\bibitem{Gup} {C. P. Gupta,} \textit{A  note on a second order three-point boundary value problem,}
J. Math. Anal. Appl. {\bf186} (1994) 277--281.

\bibitem{GupTro} {C. P. Gupta, S. I. Trofimchuk,} \textit{A sharper condition for the
solvability of a three-point second-order boundary value problem,}
J. Math. Anal. Appl {\bf205} (1997) 586--584.

\bibitem{Kras1} {M. A. Krasnosel'skii,} \textit{Positive
Solutions of Operator Equations}, Noordhoff, Groningen, The
Netherlands, 1964.

\bibitem{LianGe} {H. Lian, W. Ge,}
\textit{Solvability for second-order three-point boundary value
problems on a half-line,} Appl. Math. Lett. 19 (2006) 1000--1006.

\bibitem{LiangZhang} {S. Liang, J. Zhang,}
\textit{The existence of countably many positive solutions for
one-dimensional $p$-Laplacian with infinitely many singularities
on the half-line,}  Appl. Math. Comput. {\bf201} (2008) 210--220.

\bibitem{LiuLiFang} {Y. Liu, D. Li, and M. Fang,}
\textit{Solvability for second-order $m-$point boundary value
problems at resonance on the half-line,}  E.J.D.E. {\bf13} (2009)
1--11.

\bibitem{ORegan} {D. O'Regan,}
\textit{Theory of Singular Boundary Value Problems}, World
Scientific, Singapore, 1994.

\bibitem{Thomas} {L. H. Thomas,}
\textit{The calculation of atomic fields,} Proc. Camb. Phil. Soc.,
{\bf23} (1927) 542--548.

\bibitem{TianGe1} {Y. Tian, W. Ge,}
\textit{Positive solutions for multi-point boundary value problem on
the half-line,} J. Math. Anal. Appl. {\bf325} (2007) 1339--1349.

\bibitem{TianGe2} {Y. Tian, W. Ge,}
\textit{Existence and uniqueness of positive solutions for a bvp
with a $p$-Laplacian on the half-line,} E.J.D.E {\bf11} (2009)
1--11.

\bibitem{TianGeShan} {Y. Tian, W. Ge, and W. Shan,}
\textit{Positive solutions for three-point boundary value problem on
the half-line,} Comput. Math. with Appl. {\bf53} (2007) 1029--1039.

\bibitem{Zeid} {E. Zeidler,} \textit{Nonlinear Functional Analysis and
its Applications. Vol. I: Fixed Point Theorems,} Springer-Verlag,
New York, 1986.

\end{thebibliography}

\end{document}
