\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2010(2010), No. 58, pp. 1--12.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu  (login: ftp)}
\thanks{\copyright 2010 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2010/58\hfil A curve of stable solutions]
{A global curve of stable, positive
solutions for a $p$-Laplacian problem}

\author[B. P. Rynne\hfil EJDE-2010/58\hfilneg]
{Bryan P. Rynne} 

\address{Department of Mathematics and the Maxwell Institute for
Mathematical Sciences, Heriot-Watt University,
Edinburgh EH14 4AS, Scotland}
\email{bryan@ma.hw.ac.uk}

\thanks{Submitted August 13, 2009. Published April 28, 2010.}
\subjclass[2000]{34B15}
\keywords{Ordinary differential equations; $p$-Laplacian; \hfill\break\indent
nonlinear boundary value problems;  positive solutions; stability}

\begin{abstract}
We consider the boundary-value problem
\begin{gather*}
 - \phi_p (u'(x))' = \lambda f(x,u(x))  , \quad  x \in (0,1),\\
u(0) = u(1) = 0,
\end{gather*}
where
$p>1$ ($p \ne 2$), $\phi_p(s)  :=  |s|^{p-1} \mathop{\rm sign} s$, $s \in \mathbb{R}$,
$\lambda \ge 0$,
and the function
$f : [0,1] \times \mathbb{R} \to \mathbb{R}$ is $C^1$ and
satisfies
\begin{gather*} 
f(x,\xi) > 0, \quad (x,\xi) \in [0,1] \times \mathbb{R} ,\\
(p-1)f(x,\xi) \ge f_\xi(x,\xi) \xi ,
\quad  (x,\xi) \in [0,1] \times (0,\infty) .
\end{gather*}
These assumptions on $f$ imply that the trivial solution
$(\lambda,u)=(0,0)$ is
the only solution
with $\lambda=0$ or $u=0$,
and if $\lambda > 0$ then any solution $u$ is {\em positive},
that is, $u > 0$ on $(0,1)$.

We prove that  the set of nontrivial solutions
consists of a $C^1$ curve of positive solutions in
$(0,\lambda_{\rm max}) \times C^0[0,1]$,
with a parametrisation of the form
$\lambda \to (\lambda,u(\lambda))$,
where $u$ is a $C^1$ function defined on $(0,\lambda_{\rm max})$,
and $\lambda_{\rm max}$ is a suitable weighted eigenvalue of the
$p$-Laplacian
($\lambda_{\rm max}$ may be finite or $\infty$),
and $u$ satisfies
$$
\lim_{\lambda\to 0} u(\lambda) = 0,
\quad
\lim_{\lambda \to \lambda_{\rm max}} |u(\lambda)|_0 = \infty .
$$
We also show that for each $\lambda \in (0,\lambda_{\rm max})$
the solution $u(\lambda)$ is globally asymptotically stable,
with respect to positive solutions
(in a suitable sense).
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
%\newtheorem{cor}[thm]{Corollary}
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{remark}[theorem]{Remark}


\section{Introduction}  \label{intro.sec}

We consider the boundary-value problem
\begin{gather}
 - \phi_p (u'(x))' = \lambda f(x,u(x))  , \quad  x \in (0,1),
\label{main.eq}\\
u(0) = u(1) = 0,
\label{bc.eq}
\end{gather}
where
$p>1$ ($p \ne 2$), $\phi_p(s)  :=  |s|^{p-1} \mathop{\rm sign} s$,
$s \in \mathbb{R}$, $\lambda \ge 0$,
and the function
$f : [0,1] \times \mathbb{R} \to \mathbb{R}$ is $C^1$ and
satisfies
\begin{gather}  \label{posf.eq}
f(x,\xi) > 0, \quad (x,\xi) \in [0,1] \times \mathbb{R} ,\\  \label{fcond.eq}
(p-1)f(x,\xi) \ge f_\xi(x,\xi) \xi ,
\quad  (x,\xi) \in [0,1] \times (0,\infty) .
\end{gather}
The condition \eqref{posf.eq} ensures that the trivial solution
$(\lambda,u)=(0,0)$ is the only solution
with $\lambda=0$ or $u=0$,
and if $\lambda > 0$ then any solution $u$ is {\em positive},
that is, $u > 0$ on $(0,1)$.

In the semilinear case ($p=2$) the problem
\eqref{main.eq}--\eqref{bc.eq}
has been considered in many papers, under various hypotheses on $f$,
see for example
\cite{BIS,CR,DANSTR,HER,K1,K2,K3,LAE,SHI,SS}.
When $f$ is independent of $x$, detailed results for this case
are obtained in \cite{BIS} and \cite{LAE}.
These papers use quadrature to derive explicit formulae for a $C^1$
curve of solutions in $[0,\infty) \times C^0[0,1]$, passing through
$(\lambda,u) = (0,0)$, with a parametrisation of the form
$s \to (\lambda(s),u(s))$,
where the parameter $s = |u(s)|_0$.
The results on the shape of the solution curve are then obtained by
investigating the function $s \to \lambda(s)$.
However, when $f$ depends on $x$, such a formula for the solutions is
not available.
Despite this, curves of solutions, with similar properties to those
in the $x$ independent case, have been constructed in, for example,
\cite[Section~4]{CR}, and  \cite{HER,K1,K2,K3,LAE,SHI,SS}
(again under a variety of hypotheses on $f$).
In these papers the strategy is to use the implicit
function theorem to construct a solution curve in
$[0,\infty) \times C^0[0,1]$ by continuation away from the solution
$(\lambda,u) = (0,0)$,
and then investigate the structure of this curve directly.


The case of general $p > 1$ ($p \ne 2$) with $f$ independent of $x$ has
been considered in many recent papers using the quadrature method,
see for example,
\cite{AW1,AW2,BA,CS,KS,WY}.
In this paper we consider the general $p$ case, with $f$ dependent on
$x$, and we use the continuation approach to prove the following
results.
There exists $\lambda_{\rm max} > 0$
($\lambda_{\rm max}$ may be finite or $\infty$)
such that the set of nontrivial solutions of
\eqref{main.eq}--\eqref{bc.eq} consists of a $C^1$ curve of
globally stable, positive solutions in
$(0,\lambda_{\rm max}) \times C^0[0,1]$, with a parametrisation of the
form $\lambda \to (\lambda,u(\lambda))$,
where $u$ is a $C^1$ function defined on the interval
$(0,\lambda_{\rm max})$.
Furthermore,
$$
\lim_{\lambda\to 0} \lambda^{-p^*} u(\lambda) = - \Delta_p^{-1}(f(0)),
\quad
\lim_{\lambda \to \lambda_{\rm max}} |u(\lambda)|_0 = \infty ,
$$
(where $\Delta_p^{-1}$ is the inverse of the $p$-Laplacian operator, and
$p^* := (p-1)^{-1}$),
so the curve meets  the point $(0,0)$.
We also characterise the value of $\lambda_{\rm max}$ as a
weighted eigenvalue of the $p$-Laplacian

Under the hypotheses
\eqref{posf.eq}, \eqref{fcond.eq},
similar results have been obtained in the semilinear case $p=2$,
and in the general $p$ case with $f$ independent of $x$.
Other hypotheses on $f$ yield so called `S-shaped' curves of solutions
(the form of the parametrisation described above shows that S-shaped
curves are precluded by \eqref{posf.eq}, \eqref{fcond.eq}),
see for example \cite{BIS,K2,WY}, and the references therein.

The continuation approach relies on the use of the implicit function
theorem.
In order to apply the implicit function theorem to the above problem we
require some recent results on differentiability of the inverse
of the $p$-Laplacian. These results will be described in
Section~\ref{K_diff.sec}, and the solution curve will then be
constructed in Section~\ref{main_res.sec}.
In Section~\ref{stability.sec} it will be shown that the solutions on
this curve are globally asymptotically stable with respect to positive
solutions
(in a sense to be made precise below).
Finally, in Section~\ref{bif_curve.sec}, we briefly consider the
situation when we change the condition \eqref{posf.eq} to
allow $f(\cdot,0) = 0$
(while retaining \eqref{fcond.eq}),
and we show that a similar $C^1$ curve of solutions exists.

\section{Preliminaries}  \label{prelim.sec}


For any integer $r \ge 0$, $C^r[0,1]$ will denote the standard
Banach
space of real valued, $r$-times continuously differentiable functions
defined on $[0,1]$, with the norm $|u|_r = \sum_{i=0}^r
|u^{(i)}|_0$,
where $|\cdot|_0$ denotes the usual sup-norm on $C^0[0,1]$
(throughout, all function spaces will be real).
For any $q \ge 1$,  $L^q(0,1)$ will denote the standard Banach
space of real valued functions on $[0,1]$ whose $q$th power is
integrable, with norm $\|\cdot\|_q$.
We let $W^{1,q}(0,1)$, with norm $\|\cdot\|_{1,q}$,
denote the usual Sobolev space of absolutely continuous
functions $u$ on $[0,1]$, with derivative $u' \in L^q(0,1)$,
while $W_0^{1,q}(0,1)$ denotes the set of functions in $W^{1,q}(0,1)$
satisfying \eqref{bc.eq}.

If $F : X \to Z$ is a function between Banach spaces  $X$ and $Z$,
then $Df(x) : X \to Z$ will denote the
Fr\'echet derivative of $F$ at $x$;
partial Fr\'echet derivatives will be indicated by subscripts,
for example, $D_x G(x,y),\, D_y G(x,y)$ will denote the partial
derivatives of a function $G$ depending on $x$ and $y$.


\subsection{The $p$-Laplacian and its inverse}  \label{K_diff.sec}

Letting
$$
\mathcal{D}_p : \{ u \in C^1[0,1] :
\text{$u$ satisfies \eqref{bc.eq} and $\phi_p(u') \in W^{1,1}(0,1)$\} ,}
$$
we define the $p$-Laplacian operator $\Delta_p : \mathcal{D}_p \to L^1(0,1)$ by
$$
\Delta_p(u) = \phi_p(u')' , \quad u \in \mathcal{D}_p .
$$
This operator is {\em $(p-1)$-homogeneous}, that is,
$\Delta_p(t u) = t^{p-1} \Delta_p(u)$,
for any $t \in \mathbb{R}$ and $u \in \mathcal{D}_p$.
The following invertibility result
is well known --- see, for example,
\cite[Theorem~3.1]{BR}, \cite[Theorem~20]{HM}
(these references prove the result for periodic boundary conditions, but
the proof can readily be modified to deal with Dirichlet boundary
conditions).

\begin{theorem} \label{De_inv.thm}
For any $h \in L^1(0,1)$, the problem
\begin{equation} \label{De_h.eq}
\Delta_p(u) = h,
\quad h \in L^1(0,1) ,
\end{equation}
has a unique solution $u = \Delta_p^{-1}(h) \in \mathcal{D}_p$.
The operator
$\Delta_p^{-1} : L^1(0,1) \to C^1[0,1]$
is continuous and $p^*$-homogeneous
$($recall that $p^* := (p-1)^{-1})$.
The operator
$\Delta_p^{-1} : L^1(0,1) \to C^0[0,1]$
is compact.
\end{theorem}


Next we discuss the differentiability of the operator $\Delta_p^{-1}$.
The following result is proved in \cite[Theorem~3.4]{BR}
for the periodic case; the proof in the Dirichlet case is similar
(but simpler).
A similar result is described in Theorem~5 and Corollary~6 of
\cite{GS}, however, the arguments in the proofs in \cite{GS}
seem to be incomplete.


\begin{theorem}  \label{diff_K.thm}
For $h \in L^1(0,1)$, let $u = u(h) := \Delta_p^{-1}(h)$.
\begin{itemize}
\item[(A)] Suppose that $p>2$
and $h \in C^0[0,1]$ is such that
$u'(x) = 0 \implies h(x) \ne 0$, for  $x \in [0,1]$.
Then there exists a neighbourhood $V$ of $h$ in $C^0[0,1]$
such that:
\begin{itemize}
\item[(a)]
for $h \in V$,
$|u(h)'|^{2-p} \in L^1(0,1);$
\item[(b)]
the mapping $h \to |u(h)'|^{2-p} : V \to L^1(0,\pi_p)$
is continuous;
\item[(c)]
the mapping $\Delta_p^{-1} : V \to W_0^{1,1}(0,1)$ is $C^1$
and
\begin{equation}  \label{DK_per_diff_form.eq}
w = D\Delta_p^{-1}(h) \bar h \implies  ( |u'|^{p-2} w' )' = p^* \bar h ,
\quad \bar h \in C^0[0,1] .
\end{equation}
\end{itemize}
\item[(B)]  Suppose that $1<p<2$.
Then the mapping $\Delta_p^{-1} : L^1(0,1) \to C^1[0,1]$ is $C^1$,
and \eqref{DK_per_diff_form.eq} holds for $\bar h \in L^1(0,1)$.
\end{itemize}
\end{theorem}

\begin{remark}  \label{diff_K.rem} \rm
In either case (A) or case (B) of
Theorem~\ref{diff_K.thm}, the implication in
\eqref{DK_per_diff_form.eq} is that the function
$|u'|^{p-2} w' \in W^{1,1}(0,1)$, so that the derivative
$(|u'|^{p-2} w')'$ is defined in the $L^1$ sense (at least).
This remark also applies to the function $w$ in Lemma~\ref{F_diff.lem}
below.
\end{remark}


\subsection{Principal eigenvalues of the $p$-Laplacian}
\label{princ_evals.sec}

We  briefly consider the weighted, nonlinear eigenvalue problem
\begin{equation}  \label{princ_eval.eq}
 -\Delta_p(w) = \mu \rho \phi_p(w) ,  \quad w \in \mathcal{D}_p ,
\end{equation}
where $\mu \in \mathbb{R}$ and the weight function $\rho \in L^1(0,1)$ satisfies
$\rho \ge 0$ on $[0,1]$,
and $\rho > 0$ on a set of positive measure.
A {\em principal eigenvalue} of  \eqref{princ_eval.eq} is an eigenvalue
$\mu$  for which there is a corresponding eigenfunction $w_\mu \ge 0$.
The following result is well known --- see, for example,
\cite[Sections 3-4]{CUE}.

\begin{lemma}  \label{princ_eval.lem}
Under the above hypotheses on the weight function $\rho$ the eigenvalue
problem
\eqref{princ_eval.eq}
has a unique principal eigenvalue $\mu_0(\rho)$.
In addition, $\mu_0(\rho) > 0$,  $w_{\mu_0(\rho)}$ is positive
and
$$
\int_0^1 |w'|^p  \ge \mu_0(\rho) \int_0^1 \rho |w|^p  ,
\quad  w \in W_0^{1,p}(0,1) .
$$
If $\rho_1,\,\rho_2$ are two such weight functions,
with
\[
\text{$\rho_1 \le \rho_2$ on $[0,1]$ and
$\rho_1 < \rho_2$ on a set of positive measure,}
\]
then
$\mu_0(\rho_1) > \mu_0(\rho_2).$
\end{lemma}

We can now define the number $\lambda_{\rm max}$ which will be shown
to characterise the right-hand end point of the curve of solutions.
Define
$g : [0,1] \times (0,\infty) \to [0,\infty)$ by
\begin{equation}  \label{g_defn.eq}
g(x,\xi) := f(x,\xi)/\xi^{p-1} ,
\quad (x,\xi) \in [0,1] \times (0,\infty).
\end{equation}
It follows from \eqref{fcond.eq} that
\begin{equation}  \label{foverxideriv.eq}
g_\xi(x,\xi) =
( f_\xi(x,\xi)\xi - (p-1) f(x, \xi) )  / \xi^p
\le 0 ,
\quad (x,\xi) \in [0,1] \times (0,\infty) ,
\end{equation}
which implies that the following limits exist
\begin{equation}  \label{gainfty_def.eq}
\gamma_\infty(x) := \lim_{\xi\to\infty}  g(x,\xi)  \ge 0,
 \quad x \in [0,1],
\end{equation}
and  $\gamma_\infty \in L^\infty(0,1)$.
Now let
$$
\lambda_{\rm max} :=
\begin{cases}
\mu_0(\gamma_\infty) < \infty , & \text{if $\gamma_\infty \ne 0$,}
\\
\infty ,                    & \text{if $\gamma_\infty = 0$.}
\end{cases}
$$

\section{Main results}  \label{main_res.sec}


For any $u \in C^0[0,1]$, we define $f(u) \in C^0[0,1]$ by
$f(u)(x) = f(x,u(x))$, $x \in [0,1]$
(that is, $f$ will denote both a function and its corresponding
Nemitskii operator).
Then \eqref{main.eq}--\eqref{bc.eq}  can be rewritten as
\begin{equation}  \label{reform.eq}
 - \Delta_p(u) = \lambda f(u),
 \quad (\lambda,u) \in [0,\infty) \times \mathcal{D}_p.
\end{equation}
Clearly, $(\lambda,u) = (0,0)$ is a solution of \eqref{reform.eq}
and, by Theorem~\ref{De_inv.thm} and the fact that $f(0)>0$, the only
solution of \eqref{reform.eq}
with $\lambda=0$ or $u=0$ is $(0,0)$.
Let
\[
\mathcal{S} := \{ \text{$(\lambda,u) \in (0,\infty) \times \mathcal{D}_p$ satisfying
\eqref{reform.eq}\}.}
\]
We first prove some basic positivity properties of solutions of
\eqref{reform.eq}.


\begin{lemma} \label{pos.lem}
Every solution $(\lambda,u) \in \mathcal{S}$ satisfies:
\begin{itemize}
\item[(a)]  $u > 0$ on $(0,1)$;
\item[(b)] $u'(0) > 0,$ $u'(1) < 0$.
\end{itemize}
\end{lemma}

\begin{proof}
For any $(\lambda,u) \in \mathcal{S}$,
it follows from the differential equation \eqref{main.eq} and the
positivity condition \eqref{posf.eq} that $u$ cannot have a local
minimum in the interval $(0,1)$,
so $u \ge 0$ on $[0,1]$,
and $u'(0) \ge 0,$ $u'(1) \le 0$.
If $u'(0) = 0$ then by \eqref{main.eq} there exists $\delta > 0$ such
that
$u' < 0$ on $(0,\delta)$, which contradicts the fact that $u \ge
0$.
Hence, $u'(0) > 0$ and similarly $u'(1) < 0$,
and it then follows, by the preceding argument, that $u > 0$ on
$(0,1)$.
\end{proof}

We now prove our main result on the structure of $\mathcal{S}$.

\begin{theorem} \label{un_sol_ex.thm}
There exists a $C^1$ function
$u : (0,\lambda_{\rm max}) \to C^0[0,1]$
such that:
\begin{itemize}
\item[(a)]
$\lim_{\lambda\to 0} \lambda^{-p^*} u(\lambda) = - \Delta_p^{-1}(f(0))$
and
$\lim_{\lambda \to \lambda_{\rm max}} |u(\lambda)|_0 = \infty$
\\$($if $p^* > 1$ then $u$ extends to a $C^1$ function through $0);$
\item[(b)]
if  $\lambda \in (0,\lambda_{\rm max})$ then $u(\lambda)$ is the
unique solution of \eqref{reform.eq};
\item[(c)]
if $\lambda \ge \lambda_{\rm max}$ then \eqref{reform.eq} has no
solution.
\end{itemize}
Hence,  $\mathcal{S}$ consists precisely of the curve of solutions
$\{(\lambda,u(\lambda)) : \lambda \in (0,\lambda_{\rm max})\}$.
\end{theorem}


Theorem~\ref{un_sol_ex.thm} will be proved by a series of intermediate
steps.
We first note that, by Theorem~\ref{De_inv.thm}, equation
\eqref{reform.eq} is equivalent to the equation
\begin{equation}  \label{Flau.eq}
F(\lambda,u) := u  + \lambda^{p^*} \Delta_p^{-1}(f(u)) = 0 ,
\quad (\lambda,u) \in [0,\infty) \times C^0[0,1] ,
\end{equation}
and the function
$F: \mathbb{R} \times C^0[0,1] \to C^0[0,1]$ is continuous.
Since $f(u) > 0$, for any $u \in C^0[0,1]$,
Theorem~\ref{diff_K.thm}  yields additional properties of
$F$ which, for reference, we state in the following lemma.
This lemma also unifies the cases $1<p<2$ and $p>2$ in
Theorem~\ref{diff_K.thm}, and is tailored to our use of the implicit
function theorem below
(the lemma would be trivial in the linear case $p=2$).


\begin{lemma} \label{F_diff.lem}
If $(\lambda,u) \in \mathcal{S}$ then $F$ is $C^1$ on a
neighbourhood of $(\lambda,u)$ in $(0,\infty) \times C^0[0,1]$.
The derivative $D_u F(\lambda,u) = I + \lambda^{p^*} K_u$,
where $I$ denotes the identity on $C^0[0,1]$ and
$K_u : C^0[0,1] \to C^0[0,1]$ is defined by
$$
K_u w := D \Delta_p^{-1}(f(u)) (f_\xi(u) w) ,
\quad w \in C^0[0,1] .
$$
The operator $K_u$ is compact, so $D_u F(\lambda, u)$ is singular if and
only
the null space $N(D_u F(\lambda, u)) \ne \{0\}$.
Also, $w \in N(D_u F(\lambda, u))$ if and only if $w \in W_0^{1,1}(0,1)$
and
satisfies the linear, weighted Sturm-Liouville problem
\begin{equation}  \label{nonsing.eq}
\begin{gathered}
 - ( |u'|^{p-2} w' )' = p^* \lambda f_\xi(u) w , \\
w(0) = w(1) = 0
\end{gathered}
\end{equation}
\end{lemma}

\begin{remark} \label{F_diff.rem} \rm
By Theorem~\ref{diff_K.thm}, the coefficient function
$|u'|^{p-2}$ in \eqref{nonsing.eq} satisfies the standard
linear Sturm-Liouville hypothesis: $1/|u'|^{p-2} \in L^1(0,1)$.
Hence, in particular, if $w$ is a non-trivial solution of
\eqref{nonsing.eq} then $|u'|^{p-2} w' \in W^{1,1}(0,1)$ and
$|u'|^{p-2} w'|_i \ne 0$, $i=0,\,1$
(where $|u'|^{p-2} w'|_i$ denotes the value of the continuous function
$|u'|^{p-2} w'$ at $x=i$).
\end{remark}


\begin{proposition} \label{surjective.prop}
If $(\lambda,u) \in \mathcal{S}$ then $D_u F(\lambda, u)$
is non-singular.
\end{proposition}

\begin{proof}
Suppose that $D_u F(\lambda, u)$ is singular.
By Lemma~\ref{F_diff.lem}, there exists a non-trivial $w \in W^{1,1}(0,1)$
satisfying \eqref{nonsing.eq}, with
$|u'|^{p-2} w'|_0 > 0$.
Let
\begin{align*}
x_0 &:= \inf \{x \in (0,1] : w(x) = 0 \}
\end{align*}
($w(1)=0$, so the set on the right is non-empty).
Clearly,  $x_0 > 0$,
$w(x_0) = 0$, $w(x) > 0$ for $x \in (0,x_0)$ and
$|u'|^{p-2} w'|_{x_0} < 0$.
Now let
$$
W := - \phi_p(u') w + u |u'|^{p-2} w' .
$$
By the properties of $u$ and $w$, we have
$W \in W^{1,1}(0,1)$,  $W(0)=W(1)=0$,
and, by \eqref{bc.eq}, \eqref{fcond.eq}, and \eqref{nonsing.eq}
\begin{equation}  \label{Wd.eq}
W' = \lambda( f(u) - p^* u f_\xi(u) ) w  \ge 0 ,
\quad  \text{on $(0,x_0)$,}
\end{equation}
with strict inequality when $x > 0$ is sufficiently close to $0$
(by \eqref{posf.eq}).
Combining these results yields
$$
 0 < W(x_0) = u |u'|^{p-2} w'|_{x_0} \le 0 ,
$$
and this contradiction proves Proposition~\ref{surjective.prop}.
\end{proof}

\begin{proposition} \label{lan_lamax.prop}
Suppose that $(\lambda_n,u_n) \in \mathcal{S}$, $n = 1,2,\dots,$ satisfies
$|u_n|_0 \to \infty$.
Then $\lambda_n \to \lambda_{\rm max}$.
\end{proposition}


\begin{proof}
Without loss of generality we may suppose that $\lambda_n \to
\lambda_\infty$,
for some $\lambda_\infty \ge 0$
(we allow $\lambda_\infty = \infty$),
and $|u_n|_0 > 1$ for each $n \ge 1$.

We first suppose that $\lambda_\infty < \infty$.
Let
$$
\widetilde u_n := u_n/|u_n|_0, \quad \widetilde f_n := f(
u_n)/|u_n|_0^{p-1},
\quad  n \ge 1 .
$$
Dividing \eqref{Flau.eq} by $|u_n|_0$ and using the $p^*$-homogeneity
of $\Delta_p^{-1}$ shows that
\begin{equation} \label{w_n.eq}
 - \widetilde u_n = \Delta_p^{-1}(\lambda_n \widetilde f_n) =
\Delta_p^{-1}(\lambda_n
g(u_n) \phi_p(\widetilde u_n)) ,
 \quad  n \ge 1 .
\end{equation}
By \eqref{foverxideriv.eq},
\begin{equation} \label{tf_bound.eq}
0 \le \widetilde f_n \le
|g(\cdot,1)|_0
+
\max \{ f(x,\xi) : (x,\xi) \in [0,1]^2 \},
 \quad  n \ge 1 ,
\end{equation}
so by \eqref{w_n.eq}, Theorem~\ref{De_inv.thm},
 and the compactness of the embedding
$C^1[0,1] \to C^0[0,1]$,
we may suppose that
$\widetilde u_n \to \widetilde u_\infty$ in $C^0[0,1]$,
for some $\widetilde u_\infty \in C^0[0,1]$ with $|\widetilde
u_\infty|_0 = 1$.

We now claim that
\begin{equation}  \label{fdivmodcvges.eq}
\text{$\widetilde f_n \to \gamma_\infty \phi_p(\widetilde u_\infty)$,
\quad in $L^1(0,1)$.}
\end{equation}
Results of this form are well known, so we simply sketch a proof.
By \eqref{gainfty_def.eq}, for any $x \in [0,1]$ and  $\epsilon>0$,
there exists $C(x,\epsilon) > 0$ such that, for  $n \ge 1$,
$$
\widetilde f_n(x)
\le
\frac{C(x,\epsilon) + (\gamma_\infty(x) + \epsilon)
u_n(x)^{p-1}}{|u_n|_0^{p-1}}
\to (\gamma_\infty(x) + \epsilon) \widetilde u_\infty(x)^{p-1} .
$$
Together with  a similar lower bound, and dominated convergence
(recall \eqref{tf_bound.eq}), this proves \eqref{fdivmodcvges.eq}.


Now, letting $n \to \infty$ in \eqref{w_n.eq}
(using Theorem~\ref{De_inv.thm} and \eqref{fdivmodcvges.eq}),
yields
$$
 - \widetilde u_\infty = \Delta_p^{-1}(\lambda_\infty \gamma_\infty
\phi_p(\widetilde u_\infty)) .
$$
Since $\widetilde u_\infty \ne 0$, this implies that $\gamma_\infty
\ne 0$,
and so we have shown that
$\lambda_\infty < \infty \implies \gamma_\infty \ne 0$
(and $\lambda_\infty \ne 0$).
Furthermore, since $\widetilde u_\infty \ge 0$, it follows from
Lemma~\ref{princ_eval.lem} that
$\lambda_\infty = \mu_0(\gamma_\infty) = \lambda_{\rm max}$.

Now suppose that
$\gamma_\infty \ne 0$.
Choose a closed subinterval $I \subset (0,1)$ such that the
restriction
$\gamma_{\infty,I} := \gamma_\infty |_I \ne 0$ in $L^\infty(I)$.
Applying Lemma~\ref{princ_eval.lem} on the interval
$I$, there exists a principal eigenvalue $\mu_{\infty,I} > 0$ and a
positive eigenfunction
$w_{\infty,I} \in C^1(I)$ satisfying the  eigenvalue problem
\begin{equation}  \label{inf_eval_prb.eq}
\begin{gathered}
 - \phi_p(w_{\infty,I}')'
 = \mu_{\infty,I} \gamma_{\infty,I} \phi_p(w_{\infty,I}) ,
 \quad \text{on }I,
\\
w_{\infty,I} = 0 ,
 \quad \text{on }\partial I.
\end{gathered}
\end{equation}
Now, for any $n \ge 1$, $g(u_n) \ge \gamma_{\infty,I}$ on $I$,
so comparing  \eqref{w_n.eq} and \eqref{inf_eval_prb.eq}
shows that
if $\lambda_n > \mu_{\infty,I}$ then $w_n$ must have a zero in $I$
(by the $p$-Laplacian form of the Sturm comparison theorem,
see  \cite[Theorem~6]{HM}).
However, this contradicts the positivity of $u_n$,
so in fact $\lambda_n \le \mu_{\infty,I}$ for all $n \ge 1$.
Thus, we have shown that  $\gamma_\infty \ne 0 \implies \lambda_\infty <
\infty$,
which completes the proof of Proposition~\ref{lan_lamax.prop}.
\end{proof}

\begin{proposition} \label{g_dec_globcurve.prop}
There exists a $C^1$ function $u : (0,\lambda_{\rm max}) \to C^0[0,1]$
having the
properties $(a)$-$(b)$ described in
Theorem~\ref{un_sol_ex.thm}.
\end{proposition}

\begin{proof}
We first search for solutions of  \eqref{Flau.eq} near to
$(\lambda,u) = (0,0)$.
Writing $\lambda^{p^*} = \widetilde \lambda$, equation \eqref{Flau.eq}
becomes
$$
\widetilde F(\widetilde \lambda,u) := u + \widetilde \lambda
\Delta_p^{-1}(f(u)) = 0 , \quad (\widetilde \lambda,u) \in \mathbb{R}
\times C^0[0,1].
$$
By Theorem~\ref{diff_K.thm}, $\widetilde F$ is $C^1$ near to $(0,0)$
(since $f(0)>0$),
and  clearly
$\widetilde F(0,0) = 0$, $D_u \widetilde F(0,0) = I$.
Hence, by the implicit function theorem, there exists $\epsilon>0$ and a
$C^1$ function $\widetilde u : (-\epsilon,\epsilon) \to C^0[0,1]$, such
that
$\widetilde F(\widetilde \lambda,\widetilde u(\widetilde \lambda)) = 0$,
$\widetilde \lambda \in (-\epsilon,\epsilon)$.
Also, it is clear that $\widetilde u'(0) = - \Delta_p^{-1}(f(0))$.
Defining $u : (0,\epsilon)  \to C^0[0,1]$ by
$u(\lambda) = \widetilde u(\lambda^{p^*})$, this function has the
required
properties
on the interval $(0,\epsilon)$.
We now need to extend this function to $(0,\lambda_{\rm max})$.

Suppose that $(\lambda_0,u_0) \in \mathcal{S}$.
By Proposition~\ref{surjective.prop} the derivative $D_u
F(\lambda_0,u_0)$
is
non-singular so, by the implicit function theorem,
there exists a $C^1$ function $u$, defined on a
neighbourhood of $\lambda_0$, such that
$u(\lambda_0) = u_0$ and $F(\lambda,u(\lambda)) \equiv 0$.
Hence, the function $u$ constructed above near to $\lambda = 0$ extends
in a
$C^1$ manner to a maximal interval
$(0,\bar\lambda)$, for some $\bar\lambda \in (0,\infty]$.

\begin{lemma}  \label{ulan_nnfty.lem}
If $(\lambda_n)$ is a sequence in $(0,\bar\lambda)$ with
$\lambda_n \to \bar\lambda$,
then $|u(\lambda_n)|_0 \to \infty$.
\end{lemma}

\begin{proof}
Suppose that the sequence $(|u(\lambda_n)|_0)$ is bounded
(the same argument deals with the case where a subsequence is bounded).
First suppose that $\bar\lambda < \infty$.
Then, by \eqref{Flau.eq} and the compactness of the operator
$\Delta_p^{-1} : L^1(0,1) \to C^0[0,1]$
(Theorem~\ref{De_inv.thm}),
we may suppose that $u(\lambda_n) \to u_\infty$ in $C^0[0,1]$,
and hence, by the continuity of
$F$, $(\bar\lambda,u_\infty) \in \mathcal{S}_0$.
But now, by the implicit function theorem, the
parametrisation can be extended to the right of $\bar\lambda$, which
contradicts
the maximality of the interval $(0,\bar\lambda)$ and shows that
the case $\bar\lambda < \infty$ cannot occur.

Now suppose that $\bar\lambda = \infty$.
Then by \eqref{posf.eq}  and the boundedness of the
sequence $(|u(\lambda_n)|_0)$, there exists
$\delta > 0 $ such that $f(u_n) \ge 4\delta$ on $[0,1]$.
For each $n$, the function $u_n$ attains its maximum at a point
$x_n$.
Suppose that $x_n \ge 1/2$
(the case $x_n < 1/2$ is similar).
Then, since $u_n'(x_n) = 0$ and $(\lambda_n,u_n)$ satisfies
\eqref{main.eq},
we have
$$
\phi_p(u_n') \ge \lambda_n \delta, \quad \text{on $[0,\tfrac14]$,}
$$
and hence, by \eqref{bc.eq},
$
u_n(\tfrac14) \ge \tfrac14 \phi_{p^*}(\lambda_n \delta) .
$
Since $\lambda_n \to \infty$, this contradicts the boundedness of
$(|u(\lambda_n)|_0)$,
and so proves the lemma.
\end{proof}


Combining Proposition~\ref{lan_lamax.prop} with
Lemma~\ref{ulan_nnfty.lem}
shows that $\bar\lambda = \lambda_{\rm max}$, which completes the proof
of
Proposition~\ref{g_dec_globcurve.prop}.
\end{proof}


Let  $\mathcal{S}_0$ denote the curve of solution
$\{(\lambda,u(\lambda)) : \lambda \in (0,\lambda_{\rm max})\}$.
We now show that this is the only curve of solutions.

\begin{proposition} \label{S_eq_Sz.prop}
$\mathcal{S} = \mathcal{S}_0$.
\end{proposition}

\begin{proof}
Suppose, on the contrary, that there exists a solution
$(\lambda_1,u_1) \in \mathcal{S} \setminus \mathcal{S}_0$,
and let  $\mathcal{S}_1$ denote the connected component of $\mathcal{S}$ containing
$(\lambda_1,u_1)$.
Then the proof of
Proposition~\ref{g_dec_globcurve.prop}
shows that $\mathcal{S}_1$ is a $C^1$ curve in $(0,\infty) \times
C^0[0,1]$,
having a parametrisation of the form
$\lambda \to (\lambda,\widetilde u(\lambda))$ defined on a
maximal open interval
$I_1 \subset (0,\infty)$,
and the only possible non-zero end point of $I_1$ is $\lambda_{\rm
max}$.
Hence, the left end point of $I$ must be at $\lambda = 0$.
However, Proposition~\ref{lan_lamax.prop} shows
that the function $\widetilde u$ is bounded near to $\lambda = 0$, and
so, by
\eqref{Flau.eq}, $\lim_{\lambda \to 0} \widetilde u(\lambda) = 0$.
However, by the implicit function theorem, this would imply that
$\mathcal{S}_1 = \mathcal{S}_0$
(recall the first part of the proof of
Proposition~\ref{g_dec_globcurve.prop}),
contradicting the assumption that
$(\lambda_1,u_1) \in \mathcal{S} \setminus \mathcal{S}_0$.
\end{proof}

Combining Propositions~\ref{g_dec_globcurve.prop}
and~\ref{S_eq_Sz.prop}  completes the proof of
Theorem~\ref{un_sol_ex.thm}.

\section{Stability}  \label{stability.sec}

We now consider the following time-dependent,
initial value problem
\begin{gather}
v_t(t) = \Delta_p(v(t)) + \lambda f(v(t))  ,   \quad t \ge 0,
\label{time_dept.eq}
\\
v(0) = v_0 \ge 0 ,           \label{time_dept_nval.eq}
\end{gather}
for $\lambda \ge 0$.
Clearly, any positive solution of \eqref{reform.eq} is an equilibrium
solution of
\eqref{time_dept.eq}-\eqref{time_dept_nval.eq},
so by Theorem~\ref{un_sol_ex.thm}, for each $\lambda \in [0,\lambda_{\rm
max})$,
\eqref{time_dept.eq}-\eqref{time_dept_nval.eq}
has a unique equilibrium solution $u(\lambda)$
(here, we write $u(0)=0$).
In this section we will briefly consider the stability of this
solution.


To state precisely what we mean by a solution of
\eqref{time_dept.eq}-\eqref{time_dept_nval.eq}, we define the space
$$
S_p := C([0,\infty),W_{0,\rm w}^{1,p}(0,1))
\cap
W^{1,2}((0,\infty),L^2(0,1)) ,
$$
where $C([0,\infty),W_{0,\rm w}^{1,p}(0,1))$ denotes the space of
$W_0^{1,p}(0,1)$-valued, weakly
continuous functions on $[0,\infty)$.
The compactness of the embedding
$W_0^{1,p}(0,1) \to C^0[0,1]$ implies that
$C([0,\infty),W_{0,\rm w}^{1,p}(0,1)) \subset C([0,\infty),C^0[0,1])$.

\begin{definition} \rm
For any $\lambda \ge 0$ and
$v_0 \in \mathcal{D}_p$,
a {\em solution}  of
\eqref{time_dept.eq}-\eqref{time_dept_nval.eq}
is a function
$v(\lambda,v_0) \in S_p$
satisfying \eqref{time_dept_nval.eq},
such that,
for almost all $t \ge 0$,  $v(\lambda,v_0)(t) \in \mathcal{D}_p$ and
\eqref{time_dept.eq} holds
(in the $L^2(0,1)$ sense).
\end{definition}

We  now prove that $u(\lambda)$ is globally asymptotically stable,
with respect to positive solutions,
in the sense of the following theorem.

\begin{theorem}  \label{stab.thm}
For each $\lambda \in (0,\lambda_{\rm max})$ and $0 \le v_0 \in \mathcal{D}_p$,
the problem
\eqref{time_dept.eq}-\eqref{time_dept_nval.eq}
has a unique solution
$v(\lambda,v_0) \in S_p$.
In addition, $v(\lambda,v_0)$ satisfies:\\
$(a)\ 0 \le v(\lambda,v_0)(t) \in \mathcal{D}_p$, for each $t \ge 0;$\\
$(b)\ \lim_{t \to \infty} v(\lambda,v_0)(t) = u(\lambda)$, in
$C^0[0,1]$.
\end{theorem}


\begin{proof}
We consider fixed $\lambda \in (0,\lambda_{\rm max})$ and $0 \le v_0 \in
\mathcal{D}_p$,
and suppose, for now, that a solution $v(\lambda,v_0) \in S_p$ exists.
Then the non-negativity property in (a) is well-known, but for
convenience we give a short proof.
For any function $w$, let $w^-(x) := \min\{w(x),0\}$.
By the definition of $S_p$, the function
$\|v(\lambda,v_0)^-(\cdot)\|_2$
is absolutely continuous on $[0,\infty)$, with
$\|v(\lambda,v_0)^-(0)\|_2 = 0$
and,
writing $v$ for $v(\lambda,v_0)(t)$,
\begin{align*}
\frac12 \frac{d}{dt} \|v^- \|_2^2
&= \int_0^1 v_t  v^-
= \int_0^1 \Delta_p(v)  v^-  +  \lambda \int_0^1 f(v)  v^-
\le 0 ,
\quad \text{a.e. $t \ge 0$}
\end{align*}
(the differentiability properties follow from \cite[Prob 30.9]{ZEI}
and a slight modification of \cite[Lemma 7.6]{GT};
the final inequality follows from
an integration by parts and \eqref{posf.eq}).
Hence, $\|v(\lambda,v_0)^-(\cdot)\|_2 \equiv 0$, which proves the
non-negativity.

By the preceding argument, we need only search for non-negative
solutions.
We now define the set $W_{0,+}^{1,p}(0,1) := \{ w \in W_0^{1,p}(0,1) : w
\ge 0\}$, and
the functions
\begin{gather*}
F(x,\xi) := \int_0^\xi f(x,s)\,ds,
\quad  (x,\xi) \in [0,1] \times [0,\infty),
\\
E(w) := \frac{1}{p} \int_0^1 |w'|^p - \lambda \int_0^1 F(w),
\quad w \in W_{0,+}^{1,p}(0,1) .
\end{gather*}
The function $E : W_{0,+}^{1,p}(0,1) \to \mathbb{R}$ is continuous.

\begin{lemma}  \label{Ebdd.lem}
There exists an increasing function $M : \mathbb{R} \to (0,\infty)$ such
that,
$$
|w|_0  + \|w\|_{1,p} < M(E(w)) ,
\quad w \in W_{0,+}^{1,p}(0,1) .
$$
\end{lemma}

\begin{proof}
Suppose that there exists $R \in \mathbb{R}$ and
$0 \ne w_n \in W_{0,+}^{1,p}(0,1)$, $n = 1,2,\dots,$
such that $E(w_n) \le R$ and $\|w_n\|_{1,p} \to \infty$;
let
$\widetilde w_n := w_n / \|w_n\|_{1,p}$.
By the compactness of the embedding $W_0^{1,p}(0,1) \to C^0[0,1]$,
we may assume that $\widetilde w_n \to \widetilde w_\infty$ in
$C^0[0,1]$, for some
$0 \le \widetilde w_\infty \in C^0[0,1]$,
and it suffices to show that this leads to a contradiction.

A similar argument to the proof of \eqref{fdivmodcvges.eq}
shows that
\begin{equation}  \label{F_lim.eq}
\int_0^1
\frac{F(w_n)}{\|w_n\|_{1,p}^p}
\to
\frac{1}{p}  \int_0^1 \gamma_\infty \widetilde w_\infty^p .
\end{equation}
Now suppose that $\int_0^1 \gamma_\infty \widetilde w_\infty^p > 0$.
Then, by Lemma~\ref{princ_eval.lem}, for each $n \ge 1,$
\begin{equation}  \label{est_pos.eq}
\int_0^1 |\widetilde w_n'|^p
\ge \mu_0(\gamma_\infty) \int_0^1 \gamma_\infty \widetilde w_n^p
\to  \mu_0(\gamma_\infty) \int_0^1 \gamma_\infty \widetilde w_\infty^p
> 0,
\end{equation}
and combining \eqref{F_lim.eq} with \eqref{est_pos.eq} shows that
$E(w_n) \to \infty$
(since $\lambda < \lambda_{\rm max} = \mu_0(\gamma_\infty)$).
However, this  contradicts the initial assumption that
$E(w_n) \le R$, $n \ge 1.$
Next, suppose that $\int_0^1 \gamma_\infty \widetilde w_\infty^p = 0$,
but
$\widetilde w_\infty^p \ne 0$.
Then,  by Lemma~\ref{princ_eval.lem},
\begin{equation}  \label{est_zero.eq}
\int_0^1 |\widetilde w_n'|^p
\ge \mu_0({\bf 1}) \|\widetilde w_n\|_p^p
\to  \mu_0(1) \|\widetilde w_\infty\|_p^p
> 0
\end{equation}
(where $\bf 1$ denotes the weight function that is identically 1 on
$[0,1]$),
and combining \eqref{F_lim.eq} with \eqref{est_zero.eq} again yields
the contradiction
$E(w_n) \to \infty$.
Finally, suppose that $\widetilde w_\infty^p = 0$.
Since $\|\widetilde w_n\|_{1,p} = 1$, $n \ge 1$, this implies that
$
\int_0^1 |\widetilde w_n'|^p \to 1,
$
which  again yields a contradiction, and so  completes the proof of the
lemma.
\end{proof}

By the definition of $S_p$, \cite[Lemma~3.3]{BRE}, and the
argument in the proof of \cite[Lemma~3.1]{CF}, the function
$E(v(\lambda,v_0)(\cdot))$
is absolutely continuous and decreasing,
so by Lemma~\ref{Ebdd.lem},
any solution $v(\lambda,v_0)$ of
\eqref{time_dept.eq}-\eqref{time_dept_nval.eq} satisfies
\begin{equation}  \label{v_le_MEv.eq}
|v(\lambda,v_0)(t)|_0 < M(E(v_0)),  \quad  t \ge 0.
\end{equation}

Now,  choosing a decreasing function $\theta_0 \in C^\infty(\mathbb{R},\mathbb{R})$,
with
$$
\theta_0(s) =
\begin{cases}
1, & s \le M(E(v_0)) ,
\\
0, & s \ge 2 M(E(v_0)) ,
\end{cases}
$$
and defining
$\hat f_0 : [0,1] \times \mathbb{R} \to (0,\infty)$ by
$$
\hat f_0(x,\xi) = f(x,\theta_0(\xi)\xi), \quad (x,\xi) \in [0,1] \times
\mathbb{R},
$$
we consider the initial value problem
\begin{gather}
\hat v_t(t) = \Delta_p(\hat v(t)) + \lambda \hat f_0(\hat v(t))  ,
\quad t \ge 0,
\label{h_time_dept.eq}
\\
\hat v(0) = v_0 .           \label{h_time_dept_nval.eq}
\end{gather}
It follows from the form of $\theta_0$ and $\hat f_0$,
and  \eqref{v_le_MEv.eq},
that any solution of
\eqref{time_dept.eq}-\eqref{time_dept_nval.eq}
is a solution of
\eqref{h_time_dept.eq}-\eqref{h_time_dept_nval.eq},
and vice versa.
Now, since $\hat f_0$ is bounded and Lipschitz,
\eqref{h_time_dept.eq}-\eqref{h_time_dept_nval.eq}
has a unique solution $\hat v(\lambda,v_0) \in S_p$,
see \cite[Theorem 3.4]{AO}
(the existence and properties of solutions under the hypothesis that
$\hat f_0$ is bounded and Lipschitz  is also discussed in
\cite[Remark~2.2]{CF}).
Hence, $\hat v(\lambda,v_0)$
is in fact a solution of
\eqref{time_dept.eq}-\eqref{time_dept_nval.eq},
which we will denote by $v(\lambda,v_0)$.
In addition, \cite[Theorem 3.4]{AO} also shows that
$v(\lambda,v_0)(t) \in \mathcal{D}_p$, for all $t \ge 0$.

Finally, suppose that there exists $\epsilon > 0$ and  $t_n > 0$,
$n = 1,2,\dots,$
such that
$t_n \to \infty$ and $|v(\lambda,v_0)(t_n) - u(\lambda)|_0 \ge
\epsilon$.
By Lemma~\ref{Ebdd.lem} and compactness we may suppose that
$v(\lambda,v_0)(t_n) \to v_\infty$ in $C^0[0,1]$,
for some $v_\infty \ge 0$.
By the argument in the proof of \cite[Lemma~3.1]{CF}, this now implies
that $v_\infty$ must be an equilibrium solution of \eqref{time_dept.eq},
so that $v_\infty = u(\lambda)$.
However, this contradicts the choice of the sequence $(t_n)$ and so
completes the proof of Theorem~\ref{stab.thm}.
\end{proof}


\section{The case when $f(\cdot,0)=0$}  \label{bif_curve.sec}

In this final section we briefly consider the situation when we
replace the condition  \eqref{posf.eq} with the condition
\begin{equation}  \label{bif_posf.eq}
f(x,0) = 0, \quad f(x,\xi) > 0, \quad
(x,\xi) \in [0,1] \times (0,\infty) .
\end{equation}
We continue to assume that \eqref{fcond.eq} holds so that,
for each $x \in [0,1]$,
the function $g(x,\cdot)$ is still decreasing on $(0,\infty)$
and the following limit exists
\begin{equation}  \label{gaz_def.eq}
\gamma_0(x) := \lim_{\xi\to 0^+}  g(x,\xi)  \ge \gamma_\infty(x) ,
 \quad x \in [0,1]
\end{equation}
(recall that $g$ was defined in \eqref{g_defn.eq}).
We also impose the additional assumption:
\begin{equation}  \label{gaz_Lone.eq}
\text{$\gamma_0 \in L^1(0,1)$ \quad and \quad
$\lim_{\xi\to 0^+}  g_\xi(x,\xi)  < 0$ for $x$ in a set of positive
measure.}
\end{equation}
This implies that $\gamma_0 > \gamma_\infty$
in a set of positive measure, and so by Lemma~\ref{princ_eval.lem},
$$
\lambda_{\rm min} := \mu_0(\gamma_0) < \lambda_{\rm max}
= \mu_0(\gamma_\infty) .
$$

Under the assumption \eqref{bif_posf.eq},
$(\lambda,0)$ is a solution of \eqref{reform.eq} for all
$\lambda \in \mathbb{R}$,
and solutions need not be positive.
Hence, we define the set of positive solutions
\[
\mathcal{S}^+ := \{ \text{$(\lambda,u) \in \mathcal{S}$: $u \ne 0$ and $u \ge 0$}\}.
\]
We will show that with these assumptions a $C^1$ curve of
positive solutions of \eqref{reform.eq} exists,
similar to that constructed in Theorem~\ref{un_sol_ex.thm},
but now defined over the interval
$(\lambda_{\rm min},\lambda_{\rm max})$.

It will be useful to observe that, by the above hypotheses on $f$
(and hence $g$) any solution
$(\lambda,u) \in \mathcal{S}^+$ of \eqref{reform.eq} satisfies the equation
\begin{equation}  \label{g_in_g_form.eq}
 - \Delta_p(u) = \lambda g(u)  \phi_p (u) ,
\end{equation}
with $g(u) \in L^1(0,1)$ and
\begin{equation}  \label{bif_bnds_g.eq}
\gamma_0 \ge g(u) \ge \gamma_\infty  ,
\end{equation}
and these inequalities are strict on (possibly different) sets of
positive measure.

\begin{lemma} \label{bif_pos.lem}
The results of Lemma~\ref{pos.lem} hold for every
$(\lambda,u) \in \mathcal{S}^+$.
\end{lemma}

\begin{proof}
This follows from the arguments in the proof of
Lemma~\ref{pos.lem}  and the uniqueness of solutions of the initial
value problem for the differential equation \eqref{g_in_g_form.eq},
see for example \cite[Lemma 3.1]{BD}.
\end{proof}

We now have the following analogue of
Theorem~\ref{un_sol_ex.thm}


\begin{theorem} \label{bif_un_sol_ex.thm}
Suppose that $f$ satisfies
\eqref{fcond.eq}, \eqref{bif_posf.eq} and \eqref{gaz_Lone.eq}.
\\
There exists a $C^1$ function
$u : (\lambda_{\rm min},\lambda_{\rm max}) \to C^0[0,1]$
such that:
\begin{itemize}
\item[(a)]
$\lim_{\lambda\to \lambda_{\rm min}} u(\lambda) = 0$
and
$\lim_{\lambda \to \lambda_{\rm max}} |u(\lambda)|_0 = \infty$;
\item[(b)]
if  $\lambda \in (\lambda_{\rm min},\lambda_{\rm max})$ then
$u(\lambda)$
is the unique positive solution of \eqref{reform.eq};
\item[(c)]
if
$\lambda \in (0,\infty) \setminus (\lambda_{\rm min},\lambda_{\rm max})$
then \eqref{reform.eq} has no positive solution.
\end{itemize}
Hence,  $\mathcal{S}^+$ consists precisely of the curve of solutions
$\{(\lambda,u(\lambda)) : \lambda \in
(\lambda_{\rm min},\lambda_{\rm max})\}$.
\end{theorem}

\begin{proof}
The proof is similar to the proof of
Theorem~\ref{un_sol_ex.thm},
so we will merely sketch the necessary modifications to the preceding
arguments.
By \eqref{bif_posf.eq} and Lemma~\ref{bif_pos.lem},
if $(\lambda,u) \in \mathcal{S}^+$ then $f(u) > 0$ in $(0,1)$ and
$f(u)(0) = f(u)(1) = 0$, so it follows from
Theorem~\ref{diff_K.thm} and Lemma~\ref{pos.lem} that
Lemma~\ref{F_diff.lem}
holds here.
Also,
Propositions~\ref{surjective.prop}, \ref{lan_lamax.prop}
and their proofs hold unchanged,
while a similar argument to the proof of
Proposition~\ref{lan_lamax.prop}
also proves the following result.


\begin{proposition} \label{bif_lan_lammin.prop}
Suppose that $(\lambda_n,u_n) \in \mathcal{S}^+$, $n = 1,2,\dots,$ satisfies
$|u_n|_0 \to 0$.
Then $\lambda_n \to \lambda_{\rm min}$.
\end{proposition}


\begin{proposition} \label{bif_g_dec_globcurve.prop}
There exists a $C^1$ function
$u : (\lambda_{\rm min},\lambda_{\rm max}) \to C^0[0,1]$
having the
properties $(a)$-$(c)$ described in
Theorem~\ref{bif_un_sol_ex.thm}.
\end{proposition}

\begin{proof}
For any $\lambda > 0$, it follows from \cite[Lemma 3.3]{KLS}
that a positive solution $u$ of \eqref{g_in_g_form.eq} (if it exists)
must be unique
(using the fact that for each $x \in [0,1]$,
the function $g(x,\cdot)$ is decreasing).
Now, if $\lambda \not\in (\lambda_{\rm min},\lambda_{\rm max})$
then by  \eqref{bif_bnds_g.eq} and Lemma~\ref{princ_eval.lem},
\eqref{g_in_g_form.eq} has no positive solution.
On the other hand, if
$\lambda \in (\lambda_{\rm min},\lambda_{\rm max})$
then \cite[Theorem 6.4]{RYN} shows that \eqref{reform.eq} has a positive
solution.
The preceding argument then constructs a suitable $C^1$ function $u$
with the required properties.
\end{proof}

This completes the proof of Theorem~\ref{bif_un_sol_ex.thm}.
\end{proof}

\begin{thebibliography}{99}

\bibitem{AW2}
{\sc I. Addou, S.-H. Wang},
Exact multiplicity results for a $p$-Laplacian problem with
concave-convex-concave nonlinearities,
{\em Nonlinear Analysis}
{\bf 53} (2003), 111--137.

\bibitem{AW1}
{\sc I. Addou, S.-H. Wang},
Exact multiplicity results for a $p$-Laplacian positone problem with
concave-convex-concave nonlinearities,
{\em Electron. J. Diff. Equns.}
{\bf 2004}, No. 72.

\bibitem{AO}
{\sc G. Akagai, M. Otani},
Evolution inclusions governed by subdifferentials in reflexive Banach
spaces,
{\em J. Evolution Equations}
{\bf 4} (2004), 519--541.

\bibitem{BA}
{\sc S. Bouguima, A. Lakmeche},
Multiple solutions of a nonlinear problem involving the $p$-Laplacian,
{\em Comm. Appl. Nonlinear Anal.}
{\bf 7} (2000), 83--96.

\bibitem{BD}
{\sc P. Binding, P. Dr\'abek},
Sturm-Liouville theory for the $p$-Laplacian,
{\em  Studia Sci. Math. Hungar.}
{\bf 40} (2003), 375--396.

\bibitem{BR}
{\sc P. A. Binding, B. P. Rynne},
The spectrum of the periodic $p$-Laplacian,
{\em J. Differential Equations}
{\bf 235} (2007), 199--218.

\bibitem{BK}
{\sc H. Brezis, S. Kamin},
Sublinear elliptic equations in $\mathbb{R}^n$,
{\em  Manuscripta Math.}
{\bf 74} (1992), 87--106.

\bibitem{BRE}
{\sc H. Brezis,}
{\em Op\'erateurs Maximaux Monotones et Semi-groupes de Contraction
dans les Espace de Hilbert,}
North-Holland, Vol. 5, North-Holland, Amsterdam (1973).

\bibitem{BIS}
{\sc K. J. Brown, M. M. A. Ibrahim, R. Shivaji},
$S$-shaped bifurcation curves,
{\em Nonlinear Anal.}
{\bf 5} (1981), 475--486.

\bibitem{CS}
{\sc J. Cheng, Y. Shao},
The positive solutions of boundary value problems for a class of
one-dimensional $p$-Laplacians,
{\em Nonlinear Analysis}
{\bf 68} (2008),  883--891.

\bibitem{CF}
{\sc R. Chill, A. Fiorenza},
Convergence and decay rate to equilibrium of bounded solutions of
quasilinear parabolic equations,
{\em J. Differential Equations}
{\bf 228} (2006),  611--632.

\bibitem{CR}
{\sc M. G. Crandall, P. H. Rabinowitz},
Bifurcation, perturbation of simple eigenvalues, and linearized
stability,
{\em Arch. Rat. Mech. Anal.}
{\bf 52} (1973), 161--180.

\bibitem{CUE}
{\sc M. Cuesta,}
Eigenvalue problems for the $p$-Laplacian with indefinite weights,
{\em Electron. J. Differential Equations}
{\bf 2001} No. 33.

\bibitem{DANSTR}
{\sc E. N. Dancer},
On the structure of solutions of an equation in catalysis theory when
a parameter is large,
{\em J. Differential Equations}
{\bf 37} (1980), 404--437.

\bibitem{GS}
{\sc J. Garcia-Melian, J. Sabina de Lis,}
A local bifurcation theorem for degenerate elliptic equations with
radial symmetry,
{\em J. Differential Equns.}
{\bf 179} (2002), 27--43.

\bibitem{GT}
{\sc D. Gilbarg, N. S. Trudinger},
{\em Elliptic Partial Differential Equations of Second Order},
Springer, 1983.

\bibitem{HER}
{\sc J. Hernandez},
Qualitative methods for nonlinear diffusion equations,
{\em  Nonlinear diffusion problems (Montecatini Terme, 1985),}
47--118, Lecture Notes in Math., 1224,
Springer, Berlin, 1986.

\bibitem{HM}
{\sc Y.-X. Huang, G. Metzen,}
The existence of solutions to a class of semilinear differential
equations,
{\em Differential Integral Equations} {\bf 8} (1995), 429--452.

\bibitem{KLS}
{\sc R. Kajikiya, Y-H. Lee, I. Sim},
One-dimensional p-Laplacian with a strong singular indeﬁnite weight, I.
Eigenvalue,
{\em J. Differential Equations }
{\bf 244} (2008),  1985--2019.

\bibitem{KS}
{\sc J. Karatson, P. L. Simon},
Exact multiplicity for degenerate two-point boundary value problems
with $p$-convex nonlinearity,
{\em Nonlinear Analysis}
{\bf 52} (2003), 1569--1590.


\bibitem{K1}
{\sc P. Korman},
A global solution curve for a class of semilinear equations,
{\em  Electron. J. Differ. Equ. Conf.}
{\bf 1} (1997), 119--127.

\bibitem{K2}
{\sc P. Korman, Y. Li},
On the exactness of an S-shaped bifurcation curve,
{\em  Proc. Amer. Math. Soc.}
{\bf 127} (1999), 1011--1020.

\bibitem{K3}
{\sc P. Korman, J. Shi},
Instability and exact multiplicity of solutions of semilinear equations,
{\em  Electron. J. Differ. Equ. Conf.}
{\bf 5} (2000), 311--322.

\bibitem{LAE}
{\sc T. Laetsch},
The number of solutions of a nonlinear two point boundary value problem,
{\em Indiana Univ. Mat. J.}
{\bf 20} (1970), 1--13.

\bibitem{RYN}
{\sc B. P. Rynne},
$p$-Laplacian problems with jumping nonlinearities,
{\em J. Differential Equations}
{\bf 226} (2006), 501--524.

\bibitem{SHI}
{\sc J. Shi},
Blow up points of solution curves for a semilinear problem,
{\em Topol. Methods Nonlinear Anal.}
{\bf 15} (2000), 251--266.

\bibitem{SS}
{\sc J. Shi,  R. Shivaji},
Global bifurcations of concave semipositone problems,
{\em Evolution equations},  385--398, Lecture Notes in Pure and Appl.
Math., 234, Dekker, New York, 2003.

\bibitem{WY}
{\sc S.-H. Wang, T.-Z. Yeh},
Exact multiplicity of solutions and S-shaped bifurcation curves for the
$p$-Laplacian perturbed Gelfand problem in one space variable,
{\em J. Math. Anal. Appl.}
{\bf 342}  (2008), 1175--1191.

\bibitem{ZEI}
{\sc E. Zeidler},
{\em Nonlinear Functional Analysis and its Applications II/B},
Springer, New York (1990).


\end{thebibliography}


\end{document}

