\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2018 (2018), No. 189, pp. 1--12.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2018 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2018/189\hfil Solutions to a singular quasilinear problem]
{Exactness of the number of positive solutions to a singular
quasilinear problem}

\author[G. Anello, L. Vilasi \hfil EJDE-2018/189\hfilneg]
{Giovanni Anello, Luca Vilasi}

\address{Giovanni Anello \newline
Department of Mathematics and Computer Sciences,
Physical Sciences and Earth Sciences,
University of Messina,
Viale F. Stagno d'Alcontres 31,
98166 Messina, Italy}
\email{ganello@unime.it}

\address{Luca Vilasi \newline
Department of Mathematics and Computer Sciences,
Physical Sciences and Earth Sciences,
University of Messina,
Viale F. Stagno d'Alcontres 31,
98166 Messina, Italy}
\email{lvilasi@unime.it}

\thanks{Submitted April 19, 2018. Published November 20, 2018.}
\subjclass[2010]{34B08, 34B16, 34B18}
\keywords{Exactness; singular problem; positive solution; quadrature method}

\begin{abstract}
 We study the exact multiplicity of positive solutions to the one-dimensional
 Dirichlet problem
 \begin{gather*}
 -(|u'|^{p-2}u')' = \lambda u^{s-1} - \mu u^{r-1} \quad \text{in }
 \mathopen{]}0,1\mathclose{[}\\
 u(0) = u(1) = 0,
 \end{gather*}
 where $r\in\mathopen{]}0,1\mathclose{[}$,
 $p\in\mathopen{]}1,+\infty\mathclose{[}$, $r<s<p$ and
 $\lambda,\mu\in\mathopen{]}0,+\infty\mathclose{[}$. We shed light,
 in particular, on the case $r\in\mathopen{]}0,\min\{s, p/(p+1)\}\mathclose{[}$,
 completely determining the bifurcation diagram and solving some related open
 problems. Our approach relies upon quadrature methods.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{remark}[theorem]{Remark}
\allowdisplaybreaks

\section{Introduction}

In this article we focus on the singular one-dimensional Dirichlet problem
\begin{equation} \label{problema} %\tag{$P_{\lambda,\mu}$}
\begin{gathered}
-(|u'|^{p-2}u')' = \lambda u^{s-1} - \mu u^{r-1} \quad \text{in }
\mathopen{]}0,1\mathclose{[}\\
u(0) = u(1) = 0,
\end{gathered}
\end{equation}
where $r\in\mathopen{]}0,1\mathclose{[}$,
$p\in\mathopen{]}1,+\infty\mathclose{[}$, $r<s<p$ and
$\lambda,\mu$ are positive real parameters. This model has captured a lot
of attention over the latest years, both in the semilinear and in the
quasilinear case (see for instance
\cite{cholazmck1998some,diaher2014positive,diaherman2009branches,
horsim2009on,sim2009exact} and their references).
In the framework of chemical reactions the nonlinearity in
\eqref{problema} assumes the meaning of
a strong absorption (or endothermic) process with respect to the diffusion (see
\cite{diaherman2009branches} and the references therein).

Our aim is to explore the exact number of positive classical solutions
to \eqref{problema} by means of quadrature methods. Such a number,
established as one of the parameters of the problem varies,
depends solely on the profile of a suitable integral function,
as clarified later (see, besides the above references, also
\cite{anevil2016uniqueness,anevil2015positive,chesha2008the,fenzhage2008exact}
for applications of this procedure to several classes of one-dimensional problems).

Along this line of research, in \cite{horsim2009on} the authors studied
the exact multiplicity of positive solutions to
\begin{equation}\label{problHorSim} % \tag{$P_{r,L}$}
\begin{gathered}
-u'' = \beta - u^{r-1} \quad \text{in } \mathopen{]}-L,L\mathclose{[}\\
u(-L) = u(L) = 0,
\end{gathered}
\end{equation}
with $\beta,L>0$ and $0<r<1$. This can be viewed as a particular case
of \eqref{problema} in which $p=2$ and $s=1$ (observe that, after a suitable
change of variables, problem \eqref{problema} can be always rescaled so that
the parameters appear only in the endpoints of the interval, see
Theorem \ref{teoresattezza}). In this paper the authors complemented the
analysis carried out in \cite{cholazmck1998some}, answering some open issues
posed there and proving that, for an explicitly determined $L_0>0$,
if $0<r\leq 1/2$ there is no solution to \eqref{problHorSim}
when $L<L_0$ and there is a unique solution for $L\geq L_0$.
For $1/2<r<1$ there exists
$L_{\rm min}<L_0$ such that, for $L<L_{\rm min}$ there is no solution to
\eqref{problHorSim}, for
$L=L_{\rm min}$ there is a unique solution, for $L_{\rm min}<L\leq L_0$ there
are exactly two solutions and for $L>L_0$ there exists exactly one
solution (cf. \cite[Theorem 1]{horsim2009on}). The
picture depicted above basically reveals the existence of a critical
exponent $r^*:=1/2$ below which one has uniqueness and above which one
has multiplicity, and that two is the maximum
number of possible positive solutions to \eqref{problHorSim}.

The semilinear case with $s\neq 1$ was investigated in \cite{sim2009exact},
where the author determined all possible bifurcation diagrams that can
occur for different values of $s$ and $r$,
giving a short survey of the known results and pointing out the situations
where the exact multiplicity of the solutions is not known. The main result
is quite similar to the one obtained
for problem \eqref{problHorSim}; in this context the ``limit''
situation occurs when $s=2-2r$: for
$s\in\mathopen{]}1,2-2r]$ the bifurcation diagram is a monotone curve, while for
$s\in\mathopen{]}2-2r,2\mathclose{[}$ there exist an upper branch of
positive solutions for $L>L_{\rm min}$ and a lower one stopping at
$L=L_0>L_{\rm min}$ (see \cite[Theorem 3.2]{sim2009exact}).
The case $0<s<1$ is left open; the monotonicity of the time-map associated
with the problem is in fact proved only for $s\in\mathopen{]}1,2\mathclose{[}$
(see \cite[Lemma 3.1]{sim2009exact}).

The quasilinear counterpart \eqref{problema} was the object of an in-depth study in
\cite{diaherman2009branches}, which generalizes the results of
 \cite{diaher1999global} to the
singular framework. As clearly shown in this paper, when approaching
\eqref{problema} via quadrature methods, two different situations arise,
depending on whether
\begin{equation}\label{integraleFinfinito}
\int_0^\delta \frac{dt}{(-F(t))^\frac{p+1}{p}}= +\infty
\end{equation}
or
\begin{equation}\label{integraleFfinito}
\int_0^\delta \frac{dt}{(-F(t))^\frac{p+1}{p}}< +\infty,
\end{equation}
for some $\delta>0$ small enough, where
\begin{equation}\label{defF}
F(t):=\frac{1}{s}t^s-\frac{1}{r}t^r \quad \text{ for all } t>0.
\end{equation}
The validity of either \eqref {integraleFinfinito} or \eqref{integraleFfinito}
depends on the relationship between $r$ and $\frac{p}{p+1}$;
indeed, \eqref{integraleFinfinito} (respectively,
\eqref{integraleFfinito}) holds if and only if $r\geq\frac{p}{p+1}$
(respectively, $0<r<\frac{p}{p+1}$)
(see \cite[Lemma 2 and Proposition 4]{diaherman2009branches}).

The analysis in the first case was performed and completed in
\cite[Theorem 2(a)--(e)]{diaherman2009branches}.
 It can be summarized as follows: there exist $\lambda_1>\lambda_0>0$
such that problem \eqref{problema} admits:
\begin{itemize}
\item[(i)] no positive solutions if $\lambda<\lambda_0$;
\item[(ii)] a unique positive solution if $\lambda=\lambda_0$;
\item[(iii)] two positive solutions if
$\lambda\in\mathopen{]}\lambda_0,\lambda_1]$ and one of those corresponding
to $\lambda_1$ is a compacton, i.e. its first derivative vanishes at $t=0$
and $t=1$;
\item[(iv)] a positive solution and a continuum of non-negative solutions
compactly supported in $\mathopen{]}0,1\mathopen{[}$ if $\lambda>\lambda_1$.
\end{itemize}
The item (iv) is particularly interesting as it states that the lower branch
of the positive solutions can be continued to obtain non-negative solutions
with compact support in $\mathopen{]}0,1\mathopen{[}$.
Such a family of solutions, which vanishes on a positive-measured
subset of $\mathopen{]}0,1\mathopen{[}$ called the ``dead core'',
 turns out to be generated by the
compacton arising when $\lambda=\lambda_1$ (see
\cite[Theorem 2(e)]{diaherman2009branches} for
an explicit expression of the members of this family).

We also point out that the positive solutions $u$ of (ii), (iii) and (iv),
excluding of course the compacton, satisfy Hopf's condition on the boundary,
that is $u'(0)<0$ and $u'(1)>0$.

Condition \eqref{integraleFfinito} is more delicate to handle and the authors
obtained only the following partial result: the fact $(iv)$ continuing to hold, for
$\lambda\in\mathopen{]}\lambda_0,\lambda_1\mathclose{[}$
there exists $\delta>0$ such that there
is a unique positive solution if $r\in\mathopen{]}0,\delta\mathclose{[}$
and there are exactly two positive solutions if
$r\in\mathopen{\big]}\frac{p}{p+1}-\delta,\frac{p}{p+1}\mathclose{\big[}$.
They were able to establish the exact number of positive solutions only
when $s=2r$ and this amounts to one if $p\leq\frac{2r}{1-r}$, two if
$p>\frac{2r}{1-r}$ (see \cite[Theorem 2(f)]{diaherman2009branches}).
The reason of the success in the latter case relies on the fact that
the derivative of the auxiliary function $\psi$ defined by \eqref{defpsi} at
the unique zero $t_0$ of $F$ can be expressed in closed form in terms of the
Euler Beta function (see \cite[Proposition 4]{diaherman2009branches}).
And the knowledge of the sign of $\psi'(t_0)$, as explained later,
is crucial to draw the bifurcation diagram of \eqref{problema} when
$0<r<\frac{p}{p+1}$.

Driven by the analogy with the semilinear case, it is quite reasonable to
conjecture the existence
of a unique exponent $r^*\in\mathopen{\big]}0,\frac{p}{p+1}\mathclose{\big[}$
such that one has
exactly one solution for $0<r\leq r^*$ and exactly two for $r^*<r<\frac{p}{p+1}$
(see, in particular, the survey paper \cite{diaher2014positive}).
In this article we give a positive answer to the above conjecture,
proving the uniqueness of the zero of the equation
$\psi'(t_0(r))=0$ when $r\in\mathopen{]}0,\min\lbrace s,\frac{p}{p+1}\rbrace
\mathclose{[}$ (Theorem \ref{derivataint0}) and determining therefore the exact
multiplicity of solutions to \eqref{problema} when \eqref{integraleFfinito}
holds (Theorem \ref{teoresattezza}).
Our results are illustrated in the next section, preceded by some preliminary facts.

\section{Results}

Let $F$ be the function defined by \eqref{defF}.
Denoting by $F^{(i)}$, $i\in\mathbb{N}\cup\{0\}$, the
$i$-th derivative of $F$ (as usual, $F^{(0)}:=F$), it is easy to prove the
following property (cf. \cite[Lemma 2.1]{anevil2015positive}).

\begin{lemma}\label{segnoFederivate}
There exist $t_0,t_1\in \mathopen{]}0,+\infty\mathclose{[}$, with $t_0 > t_1$,
such that, for each $i=0,1$ one has
\begin{equation}\label{segnoFi}
 \begin{gathered}
 F^{(i)}(t)<0 \quad \text{if } t\in \mathopen{]}0,t_i\mathclose{[},\\
 F^{(i)}(t)>0 \quad \text{if } t\in \mathopen{]}t_i,+\infty\mathclose{[}.
 \end{gathered}
 \end{equation}
\end{lemma}

For the sequel, it will be useful to know the behavior in
$\mathopen{]}0,1\mathclose{[}$ of the
real functions of the next three lemmas.

\begin{lemma}\label{funzionelogtx}
Let $x>0$. The function $j(t):=\frac{\log t}{1-t^x}$ is increasing in
$\mathopen{]}0,1\mathclose{[}$, with
$$
\lim_{t\to 1^-}j(t)=-\frac{1}{x}.
$$
\end{lemma}

\begin{proof}
 One has
$$
j'(t)=\frac{t^{-1}(1-t^x)+xt^{x-1}\log t}{(1-t^x)^2}=t^x\frac{(t^{-x}-1)
+x\log t}{t(1-t^x)^2},
\quad \text{for all } t\in \mathopen{]}0,1\mathclose{[}.
$$
Putting
$$
\tilde j(t):=(t^{-x}-1)+x\log t, \quad \text{for all } t\in \mathopen{]}0,1],
$$
one has $\tilde j(1)=0$ and
$$
\tilde j'(t)=-xt^{-x-1}+xt^{-1}=xt^{-1}(1-t^{-x})<0, \quad \text{for all } t\in
\mathopen{]}0,1\mathclose{[}.
$$
Hence, $\tilde j(t)>0$ for all $t\in \mathopen{]}0,1\mathclose{[}$,
and so is $j'(t)$. The limit
follows by an immediate application of L'H\^{o}pital's rule.
\end{proof}

\begin{lemma}\label{funzionetxty}
 Let $0<x<y$. The function $k(t){:}=\frac{1-t^y}{t^x-t^y}$ is decreasing
in $\mathopen{]}0,1\mathclose{[}$, with
 $$
 \lim_{t\to 1^-}k(t)=\frac{y}{y-x}.
 $$
\end{lemma}

\begin{proof}
One has
\begin{align*}
k'(t) &=\frac{-yt^{y-1}(t^x-t^y)-(1-t^y)(xt^{x-1}-yt^{y-1})}{(t^x-t^y)^2}\\
&= -\frac{t^{x-1}}{(t^x-t^y)^2}((y-x)t^y+x-yt^{y-x}),
\end{align*}
for any $t\in \mathopen{]}0,1\mathclose{[}$. If we define
$$
\tilde k(t):=(y-x)t^y+x-yt^{y-x}, \quad \text{for any } t\in \mathopen{]}0,1],
$$
we get $\tilde k(1)=0$ and
$$
\tilde k'(t)=y(y-x)t^{y-1}-y(y-x)^{y-x-1}=y(y-x)t^{y-x-1}(t^x-1)<0,
$$
for all $t\in \mathopen{]}0,1\mathclose{[}$. As a result, $\tilde k(t)>0$
for all $t\in \mathopen{]}0,1\mathclose{[}$ and so $k'(t)<0$ for all
$t\in \mathopen{]}0,1\mathclose{[}$.
\end{proof}

\begin{lemma}\label{funzionetxty1}
Let $y\in \mathopen{]}0,+\infty[$ and $t\in\mathopen{]}0,1[$. The function
$l(x):=\frac{y-x}{t^x-t^y}$ is increasing in $[0,y[$ and
$$
\lim_{x\to y}l(x)=-\frac{1}{t^y\log t}.
$$
\end{lemma}

\begin{proof}
One has
\begin{align*}
l'(x) & =(t^x-t^y)^{-2}(-(t^x-t^y)-(y-x)t^x\log t)\\
 & = (t^x-t^y)^{-2}t^x(t^{y-x}-1-(y-x)\log t),
\end{align*}
for all $x\in \mathopen{]}0,y\mathclose{[}$. Fix $x\in[0,y[$ and put
$\tilde l(\tau):= \tau^{y-x}-1-(y-x)\log \tau$, for all $\tau\in \mathopen{]}0,1[$.
One has $\tilde l(1)=0$ and
$\tilde l'(\tau)=\tau^{-1}(y-x)(\tau^{y-x}-1)\leq 0$ for all
$\tau\in \mathopen{]}0,1[$. Therefore
$\tilde l(\tau)\geq 0$ for all $\tau\in \mathopen{]}0,1[$ and so
$l'(x)=(t^x-t^y)^{-2}t^x\tilde l(t)\geq 0$ for all $x\in \mathopen{]}0,y[$.
The limit is a consequence of L'H\^{o}pital's rule.
\end{proof}

The essential part of our approach is the analysis of the following function
 associated with \eqref{problema}:
\begin{equation}\label{defpsi}
\psi(\sigma):=\int_0^\sigma \frac{dt}{(F(\sigma)-F(t))^{1/p}}, \quad
\text{for each }
\sigma\geq t_0.
\end{equation}

In \cite{diaherman2009branches} the authors proved the following basic facts
concerning the profile of $\psi$
(cf. \cite[Propositions 1 and 4]{diaherman2009branches}).

\begin{lemma}\label{proprietapsi}
Let $r\in\mathopen{\big]}0,\frac{p}{p+1}\mathclose{\big[}$,
$p\in\mathopen{]}1,+\infty\mathclose{[}$, $r<s<p$ and let
$\psi{:}[t_0,+\infty\mathclose{[}\to\mathbb{R}$ be
the function defined by \eqref{defpsi}. Then:
\begin{itemize}
\item[(i)] $\psi\in C^0([t_0,+\infty))\cap C^1((t_0,+\infty))$;
\item[(ii)] $\lim_{\sigma\to +\infty}\psi(\sigma)=+\infty$;
\item[(iii)] if $s\geq \frac{p}{p+1}$, there exists $\delta>0$ small
 enough such that
\begin{gather*}
0 < \psi'(t_0) < +\infty \quad \text{if } r\in \mathopen{]}0, \delta\mathclose{[};\\
-\infty < \psi'(t_0) < 0 \quad \text{if } r\in \mathopen{\Big]} \frac{p}{p+1}-\delta,
\frac{p}{p+1}\mathclose{\Big[};
\end{gather*}
\item[(iv)] if $s=2r$ then
\begin{gather*}
0 < \psi'(t_0) < +\infty \quad\text{if } p<\frac{2r}{1-r};\\
\psi'(t_0)=0 \quad\text{if } p=\frac{2r}{1-r};\\
-\infty < \psi'(t_0) < 0 \quad\text{if } p>\frac{2r}{1-r}.
\end{gather*}
\end{itemize}
\end{lemma}

For $s<p/(p+1)$, Lemma \ref{proprietapsi} gives no information about the sign of
$\psi'(t_0)$ when $r$ is near $s$. The next result provides this missing detail.

\begin{lemma}\label{proprietapsi1}
Assume $s<\frac{p}{p+1}$. Then, there exists $\delta>0$ such that
$\psi'(t_0)>0$, for each $r\in \mathopen{]}s-\delta,s[$.
\end{lemma}

\begin{proof}
For any $r\in \mathopen{]}0,s\mathclose{[}$, we can write
\[
\psi'(t_0)=\frac{1}{p}\big(\frac{s}{r}\big)^{-\frac{r}{p(s-r)}}
r^{1/p}(s-r)^{-\frac{1}{p}}\hat{\phi}(r),
\]
where
\[
\hat{\phi}(r):=\int_0^1 \Big(\frac{s-r}{t^r-t^s}\Big)^{1/p}\frac{(p-s)(1-t^s)-
(p-r)(1-t^r)}{t^r-t^s}dt.
\]
For each $t\in \mathopen{]}0,1[$, using L'H\^{o}pital's rule, one has
\begin{align*}
&\lim_{r\to s^-}\Big(\frac{s-r}{t^r-t^s}\Big)^{1/p}\frac{(p-s)(1-t^s)
 - (p-r)(1-t^r)}{t^r-t^s}\\
&=\Big(-\frac{1}{t^s\log t}\Big)^{1/p}\frac{1-t^s+(p-s)t^s\log t}{t^s\log t}\\
&=(-\log t)^{-\frac{p+1}{p}}t^{-s/p}[1-t^{-s}-(p-s)\log t].
\end{align*}
Moreover, thanks to Lemma \ref{funzionetxty1}, we have the estimate
\begin{equation}\label{ineq1}
\begin{split}
&\Big|\Big(\frac{s-r}{t^r-t^s}\Big)^{1/p}
 \frac{(p-s)(1-t^s)- (p-r)(1-t^r)}{t^r-t^s}\Big|\\
&=\Big(\frac{s-r}{t^r-t^s}\Big)^{1/p}\big|p-r-\frac{s-r}{t^r-t^s}(1-t^s)\big| \\
&\leq t^{-s/p}(-\log t)^{-\frac{1}{p}}\Big(p-t^{-s}\frac{1-t^s}{\log t}\Big),
\end{split}
\end{equation}
and the function
\[
\mathopen{]}0,1\mathclose{[}\ni t\mapsto t^{-s/p}(-\log
t)^{-\frac{1}{p}}\Big(p-t^{-s}\frac{1-t^s}{\log t}\Big)
\]
is summable in $[0,1]$.

Then, by Lebesgue's dominated convergence theorem, the function $\hat{\phi}$
admits limit as $r\to s^-$ and
\[
\lim_{r\to s^-}\hat{\phi}(r)=\int_0^1(-\log
t)^{-\frac{p+1}{p}}t^{-s/p}[1-t^{-s}-1-(p-s)\log t]dt:=\bar{\phi}(s).
\]
Let us to evaluate $\bar{\phi}(s)$. An integration by parts yields
\begin{align*}
\bar{\phi}(s)
& =\int_0^1(-\log t)^{-\frac{p+1}{p}}t^{-s/p}(1-t^{-s}-(p-s)\log
t)dt\\
&=\int_0^1(-\log t)^{-\frac{p+1}{p}}t^{-1}
 \big(t^{1-\frac{s}{p}}-t^{1-s-\frac{s}{p}}\big)dt
 +(p-s)\int_0^1(-\log t)^{-\frac{1}{p}}t^{-s/p}dt\\
&=\big[p(-\log t)^{-\frac{1}{p}}\big(t^{1-\frac{s}{p}}-t^{1-s-\frac{s}{p}}\big)
 \big]_0^1\\
&\quad -p\int_0^1(-\log
 t)^{-\frac{1}{p}}\big[\big(1-\frac{s}{p}\big)t^{-s/p}
 -\big(1-s-\frac{s}{p}\big)t^{-s-\frac{s}{p}}\big]dt\\
&\quad +(p-s)\int_0^1(-\log t)^{-\frac{1}{p}}t^{-s/p}dt\\
&=\int_0^1(-\log t)^{-\frac{1}{p}}(p-s(p+1))t^{-s\frac{p+1}{p}}dt>0.
\end{align*}
This completes the proof.
\end{proof}

Now, we add a new piece of information on the shape of $\psi$ showing that
it possesses at most one critical point in $\mathopen{]}t_0, +\infty\mathclose{[}$.

\begin{lemma}\label{uniczeroderivata}
The equation $\psi'(\sigma)=0$ admits at most one solution in
 $\mathopen{]}t_0,+\infty[$.
\end{lemma}

\begin{proof}
After a change of variables, we can write $\psi$ as
\begin{equation}\label{cambiovarpsi}
\psi(\sigma) = \int_0^1 \frac{\sigma}{(F(\sigma)-F(\sigma t))^{1/p}}dt
= \int_0^1 \frac{dt}{(\eta_t(\sigma))^{1/p}},
\end{equation}
for any $\sigma\in[t_0,+\infty\mathclose{[}$, where
$$
\eta_t(\sigma):=\frac{1}{s}(1-t^s)\sigma^{s-p}-\frac{1}{r}(1-t^r)\sigma^{r-p}
$$
for all $t\in[0,1]$ and for all $\sigma\in [t_0,+\infty[$.
Then, the first two derivatives of $\psi$ take the form
\begin{align*}
\psi'(\sigma)&=-\frac{1}{p}\int_0^1(\eta_t(\sigma))^{-\frac{1}{p}-1}
\Big(\frac{s-p}{s}(1-t^s)\sigma^{s-p-1}
-\frac{r-p}{r}(1-t^r)\sigma^{r-p-1}\Big)dt,\\
\psi''(\sigma)
&= \frac{1}{p}\big(1+\frac{1}{p}\big)
 \int_0^1(\eta_t(\sigma))^{-\frac{1}{p}-2}\Big(\frac{s-p}{s}(1-t^s)\sigma^{s-p-1}\\
&\quad -\frac{r-p}{r}(1-t^r)\sigma^{r-p-1}\Big)^2 dt\\
&\quad -\frac{1}{p}\int_0^1(\eta_t(\sigma))^{-\frac{1}{p}-1}
 \Big(\frac{s-p}{s}(s-p-1)(1-t^s)\sigma^{s-p-2} \\
&\quad - \frac{r-p}{r}(r-p-1)(1-t^r)\sigma^{r-p-2}\Big)dt.
\end{align*}
Now, fix $K\in\mathopen{]}p+1-s, p+1-r\mathclose{[}$ and put
\begin{equation}\label{phi}
\chi(\sigma):=\sigma \psi''(\sigma) + K \psi'(\sigma),
\end{equation}
for all $\sigma\in[t_0,+\infty[$. For any $\sigma\in [t_0,+\infty[$ we
can deduce the  lower estimate
\begin{align*}
\chi(\sigma)
&>-\frac{1}{p}\int_0^1(\eta_t(\sigma))^{-\frac{1}{p}-1}
\Big(\frac{s-p}{s}(s-p-1)(1-t^s)\sigma^{s-p-1} \\
&\quad -\frac{r-p}{r}(r-p-1)(1-t^r)\sigma^{r-p-1}\Big) dt\\
&\quad-\frac{K}{p}\int_0^1(\eta_t(\sigma))^{-\frac{1}{p}-1}
\Big( \frac{s-p}{s}(1-t^s)\sigma^{s-p-1}
 - \frac{r-p}{r}(1-t^r)\sigma^{r-p-1}\Big) dt\\
&=\frac{1}{p}\int_0^1(\eta_t(\sigma))^{-\frac{1}{p}-1}
 \Big( \frac{p-s}{s}(K-p-1+s)(1-t^s)\sigma^{s-p-1} \\
&\quad - \frac{p-r}{r}(K-p-1+r)(1-t^r)\sigma^{r-p-1}\Big) dt\;
\end{align*}
and the choice of $K$ forces $\chi(\sigma)>0$. At this point,
solving \eqref{phi} for $\psi'$, one has
$$
\psi'(\sigma)=\frac{1}{\sigma^K}\Big(\int_{t_0}^\sigma \tau^{K-1}\chi(\tau)d\tau
+ c\Big)
$$
for all $\sigma\in [t_0,+\infty[$ and for a suitable constant $c\in\mathbb{R}$.
In view of the positivity of $\chi$, $\psi'$ admits at most one zero in
$]t_0,+\infty[$ and the conclusion is achieved.
\end{proof}

Lemma \ref{proprietapsi} leaves the question whether
$r\mapsto \psi'(t_0(r))$ changes sign at most
once in $]0,\min\{s,\frac{p}{p+1}\}[$ open. In the next theorem, we answer
this question affirmatively. This allows us to completely compute the number
of solutions to $(P_{\lambda, \mu})$ as the exponent $r$ varies and,
in particular, to extend the result in
\cite{sim2009exact} to the case $s>0$.

\begin{theorem}\label{derivataint0}
Let $p\in\mathopen{]}1,+\infty\mathclose{[}$ and
$s\in\mathopen{]}0,p\mathclose{[}$. Denote by
$\phi{:}]0,\min\{s,\frac{p}{p+1}\}[\to\mathbb{R}$ the function
\begin{equation}\label{defphi}
\phi(r):=\psi'(t_0(r)), \quad \text{for each }
r\in ]0,\min\{s,\frac{p}{p+1}\}[.
\end{equation}
Then
\begin{itemize}
\item[(i)] if $s\geq \frac{p}{p+1}$, the equation $\phi(r)=0$ admits
exactly one zero, say $r^*$, in $]0,\frac{p}{p+1}[$,
with $\phi$ positive in $]0,r^*[$ and negative in $]r^*,\frac{p}{p+1}[$;
\item[(ii)] if $s<\frac{p}{p+1}$, the function $\phi$ is positive in $]0,s[$.
\end{itemize}
\end{theorem}

\begin{proof}
To obtain conclusions (i) and (ii) it is sufficient to show that $\phi$ changes sign
at most once in $]0,\min\{s,\frac{p}{p+1}\}[$. Indeed, if this is true,
since $\phi\in C^0(]0,\min\{s,\frac{p}{p+1}\}[)$, conclusion
(i) follows from (iii) of Lemma \ref{proprietapsi}, and conclusion (ii)
 is a consequence of Lemma \ref{proprietapsi1}.

Taking \eqref{cambiovarpsi} into account we obtain
\begin{align*}
\psi'(\sigma)
& = \frac{\psi(\sigma)}{\sigma}
 +\int_0^1 \frac{\sigma(F'(\sigma)-tF'(\sigma t))}
 {p(F(\sigma)-F(\sigma t))^\frac{p+1}{p}}dt\\
& = \frac{1}{p}\int_0^1 \frac{p(F(\sigma)-F(\sigma t))
 - \sigma( F'(\sigma)-tF'(\sigma t))}{(F(\sigma)-F(\sigma t))^\frac{p+1}{p}}dt\\
& = \frac{1}{p}\int_0^1
\frac{\frac{p-s}{s}\sigma^s(1-t^s)-\frac{p-r}{r}\sigma^r(1-t^r)}{(F(\sigma)
 -F(\sigma t))^\frac{p+1}{p}}dt
\end{align*}
for any $\sigma\in\mathopen{]}t_0,+\infty\mathclose{[}$. Taking the limit as
$\sigma\to t_0(r)$ in the above expression and bearing in mind that
$t_0(r)=( \frac{s}{r})^{1/(s-r)}$, we therefore obtain
\begin{align*}
\phi(r)
& = \frac{1}{p}\int_0^1 \frac{\frac{p-s}{s} t_0(r)^s (1-t^s)
 -\frac{p-r}{r}t_0(r)^r(1-t^r)}{(-F(t_0(r) t))^\frac{p+1}{p}}dt\\
& = \frac{1}{p} t_0(r)^{-\frac{r}{p}}
 \int_0^1 \frac{\frac{p-s}{s} t_0(r)^{s-r} (1-t^s)
 -\frac{p-r}{r} (1-t^r)}{\left( \frac{t^r}{r}-t_0(r)^{s-r}\frac{t^s}{s}\right)
 ^\frac{p+1}{p}}dt\\
& = \frac{1}{p}\left(\frac{s}{r}\right)^{-\frac{r}{p(s-r)}}r^{1/p}
 \int_0^1 \frac{(p-s) (1-t^s)- (p-r)(1-t^r)}{\left(t^r-t^s\right)
 ^\frac{p+1}{p}}dt\\
& = \frac{1}{p}\left(\frac{s}{r}\right)^{-\frac{r}{p(s-r)}}r^{1/p}
 (s-r) \int_0^1 \frac{\frac{p-s}{s-r} (1-t^s)
 - \frac{p-r}{s-r}(1-t^r)}{\left(t^r-t^s\right) ^\frac{p+1}{p}}dt.
\end{align*}

In light of the above manipulations, to conclude the proof it is equivalent
to show that $\tilde{\phi}(r)=0$ has at most one solution in
$]0,\min\{s,\frac{p}{p+1}\}[$, where
\begin{equation}\label{defphitilde}
\tilde{\phi}(r):= \int_0^1 \frac{\frac{p-s}{s-r} (1-t^s)-
\frac{p-r}{s-r}(1-t^r)}{\left(t^r-t^s\right) ^\frac{p+1}{p}}dt
\end{equation}
for any $r\in ]0,\min\{s,\frac{p}{p+1}\}[$.
 Now, the integrand function
in \eqref{defphitilde} can be expressed as
$$
\frac{\frac{p-s}{s-r}(1-t^s)-\frac{p-r}{s-r}(1-t^r)}
{(t^r-t^s)^{\frac{p+1}{p}}}=(1-t^s)^{-\frac{1}{p}}g(t,r)
$$
for all $t\in\mathopen{]}0,1\mathclose{[}$, where
$$
g(t,r):=h(t,r)\big(\frac{p-r}{s-r}-h(t,r)^p\big)
$$
and
\begin{equation}\label{defh}
h(t,r):=\Big(\frac{1-t^s}{t^r-t^s}\Big)^{1/p}.
\end{equation}
By an easy calculation we get, for all $t\in\mathopen{]}0,1\mathclose{[}$,
$$
\frac{\partial h}{\partial r}(t,r)
=-\frac{h(t,r)}{p}\cdot\frac{\log t}{1-t^{s-r}}
$$
and therefore
\begin{align*}
 \frac{\partial g}{\partial r}(t,r)
& =\frac{\partial h}{\partial r}(t,r) \Big(\frac{p-r}{s-r}-h(t,r)^p\Big)\\
&\quad + h(t,r)\Big(\frac{p-s}{(s-r)^2} - p
 h(t,r)^{p-1}\frac{\partial h}{\partial r}(t,r)\Big)\\
&= -\frac{h(t,r)}{p}\cdot\frac{\log t}{1-t^{s-r}}
\Big(\frac{p-r}{s-r}-h(t,r)^p\Big)\\
 &\quad + h(t,r)\Big(\frac{p-s}{(s-r)^2}+ h(t,r)^p\frac{\log t}{1-t^{s-r}}\Big)\\
 & =\Big(\frac{p+1}{p}\Big)\frac{\log t}{1-t^{s-r}} h(t,r)^{p+1}\\
 &\quad - \Big(\frac{p-r}{p(s-r)}\cdot\frac{\log t}{1-t^{s-r}}
 -\frac{p-s}{(s-r)^2}\Big)h(t,r).
\end{align*}

Now, consider the function
$$
\alpha(t):=\big(1-\frac{r}{p}\big)\frac{1}{s-r}
 -\big(1+\frac{1}{p}\big)\frac{1-t^s}{t^r-t^s},
\quad t\in \mathopen{]}0,1\mathclose{[}.
$$
In view of Lemma \ref{funzionetxty}, $\alpha$ is strictly increasing in
$\mathopen{]}0,1\mathclose{[}$ and one has
\begin{gather*}
\lim_{t\to 0^+}\alpha(t) =-\infty,\\
\lim_{t\to 1^-}\alpha(t)
=\big(1-\frac{r}{p}\big)\frac{1}{s-r}-\big(1+\frac{1}{p}\big)\frac{s}{s-r}.
\end{gather*}
Two possible situations, related to the range of $s$, occur.
\begin{itemize}
\item[(a)] If $s\geq \frac{p}{p+1}$, then $\alpha(t)<0$ for all
$t\in \mathopen{]}0,1\mathclose{[}$;

\item[(b)] if $0<s<\frac{p}{p+1}$, then $\alpha(t)<0$ for all
$t\in \mathopen{]}0,1\mathclose{[}$ if $p-(p+1)s<s$ and $p-(p+1)s<r<s$,
 while the equation $\alpha(t)=0$ admits a unique root in
$\mathopen{]}0,1\mathclose{[}$, say $t(r)$, provided that
$0<r<\min\{s,p-(p+1)s\}$.
Moreover, $\alpha$ is negative in $\mathopen{]}0,t(r)\mathclose{[}$ and positive in
$\mathopen{]}t(r),1\mathclose{[}$.
\end{itemize}
Notice that in case (b), by the implicit function theorem the function
$\mathopen{]}0,s\mathclose{[}\ni r\mapsto t(r)$ is continuous in
$\mathopen{]}0,\min\{s,p-(p+1)s\}\mathclose{[}$. Moreover, when $p-(p+1)s\in
\mathopen{]}0,s\mathclose{[}$, it is easy to see that
\begin{equation} \label{5}
\lim_{r\to p-(p+1)s^-}t(r)=1.
\end{equation}

Next, define $M(r):=1$, for all $r\in ]0,s[$,
if $s\geq \frac{p}{p+1}$; for $s<\frac{p}{p+1}$ set instead
$$
M(r):=\begin{cases}
 -\frac{(s-r)\log t(r)}{1-t(r)^{s-r}} & \text{if } 0<r< \min\{s,p-(p+1)s\}\\
 1 & \text{if } p-(p+1)s<s \text{ and } p-(p+1)s\leq r<s.
 \end{cases}
$$
As a function of the variable $r$, $M$ turns out to be defined in the whole interval
$\mathopen{]}0,s\mathclose{[}$ and it is continuous there as, by \eqref{5} and Lemma
\ref{funzionelogtx}, one has
$$
\lim_{r\to p-(p+1)s^-}M(r)=1, \quad \text{whenever } p-(p+1)s\in
\mathopen{]}0,s\mathclose{[}.
$$
Also, observe that again by Lemma \ref{funzionelogtx} we can infer that
\begin{equation}\label{M}
M(r)\geq 1, \quad \text{for all } r\in \mathopen{]}0,s\mathclose{[},
\end{equation}
and that the function
$$
\beta(t):=\frac{M(r)}{s-r}+\frac{\log t}{1-t^{s-r}}, \quad
t\in \mathopen{]}0,1\mathclose{[},
$$
satisfies the same properties (a) and (b) as the function $\alpha$. Hence,
\begin{equation}\label{phipsi}
\alpha(t)\beta(t)> 0, \quad\text{for all but at most one }
t\in \mathopen{]}0,1\mathclose{[}.
\end{equation}
Therefore, taking \eqref{defh} into account and denoting for brevity
$h(t,r)$ by $h$, we obtain
\begin{align*}
& M(r)\frac{p+1}{p(s-r)}g(t,r)-\frac{\partial g}{\partial r}(t,r)\\
&=-\frac{p+1}{p}\Big(\frac{M(r)}{s-r}+\frac{\log t}{1-t^{s-r}}\Big)h^{p+1} \\
&\quad +\Big(\frac{p-r}{p(s-r)}\cdot\frac{\log t}{1-t^{s-r}}
  - \frac{p-s}{(s-r)^2} + M(r)\frac{(p+1)(p-r)}{p(s-r)^2}\Big)h\\
&=\Big[-\big(1+\frac{1}{p}\big)\Big(\frac{M(r)}{s-r}+\frac{\log t}{1-t^{s-r}}\Big)
 \frac{1-t^s}{t^r-t^s} + \big(1-\frac{r}{p}\big)
 \frac{1}{s-r}\cdot\frac{\log t}{1-t^{s-r}} \\
&\quad +\frac{1}{(s-r)^2}\big( M(r)(p-r)-(p-s)\big) 
 + \big(1-\frac{r}{p}\big)\frac{M(r)}{(s-r)^2}\Big] h\\
&=\Big[-\big(1+\frac{1}{p}\big)\Big(\frac{M(r)}{s-r}+\frac{\log t}{1-t^{s-r}}\Big)
 \frac{1-t^s}{t^r-t^s} +\big(1-\frac{r}{p}\big)
 \frac{1}{s-r}\Big(\frac{\log t}{1-t^{s-r}} \\
&\quad +\frac{M(r)}{s-r}\Big) +\frac{M(r)(p-r)-(p-s)}{(s-r)^2}\Big] h\\
&=\Big\{\Big(\frac{M(r)}{s-r}+\frac{\log t}{1-t^{s-r}}\Big)
 \Big[-\big(1+\frac{1}{p}\big)\frac{1-t^s}{t^r-t^s}
 +\big(1-\frac{r}{p}\big)\frac{1}{s-r}\Big] \\
&\quad +\frac{M(r)(p-r)-(p-s)}{(s-r)^2}\Big\}h.
\end{align*}
Noticing that by \eqref{M} one has
$$
\frac{M(r)(p-r)-(p-s)}{(s-r)^2}\geq \frac{1}{s-r}>0,
$$
and recalling \eqref{phipsi}, we can continue the previous chain of 
computations to finally arrive at
\begin{equation}\label{catenadisfinale}
\begin{split}
& M(r)\frac{p+1}{p(s-r)}g(t,r)-\frac{\partial g}{\partial r}(t,r)\\
&\geq\Big(\frac{M(r)}{s-r}+\frac{\log
 t}{1-t^{s-r}}\Big)\Big[-\big(1+\frac{1}{p}\big)
 \frac{1-t^s}{t^r-t^s}+\big(1-\frac{r}{p}\big)\frac{1}{s-r}\Big] h\\
&=\beta(t)\alpha(t)h
> 0
\end{split}
\end{equation}
for all but at most one $t\in \mathopen{]}0,1\mathclose{[}$.
 Multiplying \eqref{catenadisfinale}
by $(1-t^s)^{-\frac{1}{p}}$ and integrating between 0 and 1, we obtain
$$
 M(r)\frac{p+1}{p(s-r)}\tilde{\phi}(r) - \tilde{\phi}'(r)> 0
$$
for any $r\in ]0,\min\{s,\frac{p}{p+1}\}[$.
 If we set
\begin{equation}\label{defomega}
\omega(r):= M(r)\frac{p+1}{p(s-r)}\tilde{\phi}(r) - \tilde{\phi}'(r)
\end{equation}
and solve \eqref{defomega} for $\tilde{\phi}$ we get, for a fixed 
$r_0\in ]0,\min\{s,\frac{p}{p+1}\}[$,
$$
\tilde{\phi}(r)= e^{\frac{p+1}{p}\int_{r_0}^r \frac{M(\sigma)}{s-\sigma}d\sigma}
\Big(-\int_{r_0}^r \omega(\sigma) e^{
-\frac{p+1}{p}\int_{r_0}^\sigma\frac{M(\tau)}{s-\tau}d\tau} d\sigma + c\Big)
$$
for some $c\in\mathbb{R}$. Since $\omega> 0$, $\tilde{\phi}$ 
may vanish at most once in $]0,\min\{s,\frac{p}{p+1}\}[$
 and the proof is complete. 
\end{proof}

Using all the information obtained up to now, it is not difficult to 
derive the following result concerning the exactness of the number of positive 
solutions to \eqref{problema}.

\begin{theorem}\label{teoresattezza}
Let $r^*\in\mathopen{\big]}0,\frac{p}{p+1}\mathclose{\big[}$ be as in 
Theorem \ref{derivataint0} and let
\begin{equation}\label{deflambda1}
\lambda_1:=\mu^\frac{p-s}{p-r}\Big(2\Big(\frac{p-1}{p}\Big)^{1/p}
\int_0^{t_0}\frac{dt}{(-F(t))^{1/p}}\Big)^\frac{p(s-r)}{p-r}.
\end{equation}
Then:
\begin{itemize}
\item[(i)] if $s\geq \frac{p}{p+1}$ and $r\in \mathopen{]}0,r^*]$, or
if $s<\frac{p}{p+1}$, problem \eqref{problema} admits no positive solution 
for $\lambda<\lambda_1$ and a unique positive solution for $\lambda\geq\lambda_1$;
\item[(ii)] if $s\geq \frac{p}{p+1}$ and 
$r\in \mathopen{\big]}r^*,\frac{p}{p+1}\mathclose{\big[}$, 
there exists $\lambda_0\in \mathopen{]}0,\lambda_1\mathclose{[}$ 
such that problem \eqref{problema} admits no positive solution for 
$\lambda<\lambda_0$; a unique solution for $\lambda=\lambda_0$; 
two solutions for $\lambda_0<\lambda\leq\lambda_1$ and one solution for 
$\lambda>\lambda_1$.
\end{itemize}
\end{theorem}

\begin{remark} \rm 
As pointed out in the introduction, the solution corresponding to $\lambda_1$ 
in case (i), and one of the solutions corresponding to $\lambda_1$ in case (ii), 
are of compacton-type. In both cases the compacton generates, for 
$\lambda>\lambda_1$, a continuum of non-negative solutions
to \eqref{problema} compactly supported in $\mathopen{]}0,1\mathclose{[}$. 
In addition, the positive solutions, excluding the compactons, satisfy Hopf's 
condition at $t=0$ and $t=1$.
\end{remark}

\begin{proof}[Proof of Theorem \ref{teoresattezza}]
 Let $L>0$ and consider the problem
\begin{equation} \label{problemaL} % \tag{$P_{L}$}
\begin{gathered}
-(|u'|^{p-2}u')' = u^{s-1} - u^{r-1} \quad \text{in } 
 \mathopen{]}0,L\mathclose{[}\\
u(0) = u(L) = 0. 
\end{gathered}
\end{equation}
It is not difficult to verify that the rescaling
\begin{equation}\label{rescaling}
u(x)=\left( \frac{\mu}{\lambda}\right)^\frac{1}{s-r} v\Big(\mu^\frac{s-p}{p(s-r)}
\lambda^\frac{p-r}{p(s-r)}x\Big)
\end{equation}
turns any solution u to \eqref{problema} into a solution $v$ to \eqref{problemaL}, 
with
\begin{equation}\label{definzL}
L:=\mu^\frac{s-p}{p(s-r)} \lambda^\frac{p-r}{p(s-r)},
\end{equation}
and vice versa. Therefore let us focus on problem \eqref{problemaL}. 
It is well-known that any positive solution $u$ to \eqref{problemaL} 
is symmetric in $[0,L]$ and therefore
$$
\sigma:=\|u\|_\infty=u\big(\frac{L}{2}\big).
$$
Multiplying the first equation in \eqref{problemaL} by $u'$ and integrating from
$x\in\mathopen{]}0,L/2[$ to $L/2$ we arrive at
\begin{equation}\label{relazcompacton}
\frac{p-1}{p}|u'(x)|^p = F(\sigma)-F(u(x)).
\end{equation}
If we plug the value $x=0$ in the above relation, we deduce that $F(\sigma)\geq 0$, 
and, in  light of Lemma \ref{segnoFederivate}, that $\sigma \geq t_0$.
Being $u$ increasing in $]0,L/2[$, from \eqref{relazcompacton} one has
\begin{equation}\label{relazcompacton1}
u'(x)=\Big(\frac{p}{p-1}\Big)^{1/p}\left(F(\sigma)-F(u(x))\right)^{1/p}
\end{equation}
and therefore
$$
\int_0^x \frac{u'(y)}{(F(\sigma)-F(u(y)))^{1/p}}dy 
= \Big(\frac{p}{p-1}\Big)^{1/p} x.
$$
Putting $u(y)=t$ we get
\begin{equation}\label{udefimplicita}
\int_0^{u(x)} \frac{dt}{(F(\sigma)-F(t))^{1/p}} 
= \Big(\frac{p}{p-1}\Big)^{1/p} x,
\end{equation}
which evaluated at $x=L/2$ yields
\begin{equation}\label{legamesigmaLp}
\psi(\sigma) = \frac{L}{2}\Big(\frac{p}{p-1}\Big)^{1/p}.
\end{equation}
Hence the number of positive solutions to \eqref{problemaL} amounts to 
the number of solutions to \eqref{legamesigmaLp}, as $L$ varies in 
$\mathopen{]}0,+\infty\mathclose{[}$.

Starting from these premises, let us prove (i). 
If $s\geq \frac{p}{p+1}$ and $r\in \mathopen{]}0,r^*]$,
or if $s<\frac{p}{p+1}$, thanks to Theorem \ref{derivataint0}, $\psi$ is
increasing in a right neighbourhood of $t_0$ and, by Lemmas \ref{proprietapsi} and
\ref{uniczeroderivata}, it cannot admit any maximum. 
So the number of solutions to \eqref{legamesigmaLp} is zero if 
$L<L_1:=2(\frac{p-1}{p})^{1/p}\psi(t_0)$, one if
$L\geq L_1$. Now, let us pass to (ii). 
If $s\geq \frac{p}{p+1}$ and 
$r\in\mathopen{]}r^*,\frac{p}{p+1}\mathclose{[}$,
again due to Theorem \ref{derivataint0}
and Lemmas \ref{proprietapsi} and \ref{uniczeroderivata}, $\psi$ has a
 global minimum in $\mathopen{]}t_0, +\infty\mathclose{[}$. 
So \eqref{legamesigmaLp} will admit no solution for
$L<L_0:=2\left(\frac{p-1}{p}\right)^{1/p}\min_{\mathopen{]}t_0,
+\infty\mathclose{[}}\psi$,
one solution for $L=L_0$ and $L>L_1$ and two solutions for $0<L_0\leq L_1$. 
The conclusion then
follows using the rescaling \eqref{rescaling}-\eqref{definzL}. 
\end{proof}

\subsection*{Acknowledgements} 
The authors are members of the Gruppo Nazionale per l'Analisi
Matematica, la Probabilit\`a e le loro Applicazioni (GNAMPA) of
the Istituto Nazionale di Alta Matematica (INdAM).

\begin{thebibliography}{00}

\bibitem{anevil2016uniqueness}
G.~Anello, L.~Vilasi;
 \emph{Uniqueness of positive and compacton-type solutions for a
resonant quasilinear problem}. 
Topol. Methods Nonlinear Anal., \textbf{49(2)} (2017), 565--575.

\bibitem{anevil2015positive}
G.~Anello, L.~Vilasi;
 \emph{Positive and compacton-type solutions for a quasilinear two-point
boundary value problem}. J. Math. Anal. Appl., \textbf{431} (2015), 429--439.

\bibitem{chesha2008the}
J.~Cheng, Y.~Shao;
\emph{The positive solutions of boundary value problems for a class of
one-dimensional $p$-Laplacians}. Nonlinear Anal., \textbf{68}(\textbf{4}) 
(2008), 883--891.

\bibitem{cholazmck1998some}
Y. S. Choi, A. C. Lazer,  P. J. Mckenna;
\emph{Some remarks on a singular elliptic boundary value problem}. 
Nonlinear Anal., \textbf{32} (1998), 305--314.

\bibitem{diaher2014positive}
J. I. D\'iaz, J.~Hern\'andez;
\emph{Positive and free boundary solutions to singular
nonlinear elliptic problems with absorption: an overview and open problems}. 
Electron. J. Differ. Equ., Conf. 21 (2014), pp. 32--44.

\bibitem{diaher1999global}
J. I. Diaz, J.~Hern\'andez;
\emph{Global bifurcation and continua of nonnegative
solutions for a quasilinear elliptic problem}. C. R. Acad. Sci. Paris, S\'{e}rie I,
\textbf{329} (1999), 587--592.

\bibitem{diaherman2009branches}
J. I. D\'iaz,  J.~Hern\'andez,  F.J. Mancebo;
\emph{Branches of positive and free boundary
solutions for some singular quasilinear elliptic problems}. 
J. Math. Anal. Appl., \textbf{352} (2009), 449--474.

\bibitem{fenzhage2008exact}
M.~Feng, X.~Zhang,  W.~Ge;
\emph{Exact number of solutions for a class of two-point boundary
value problems with one-dimensional $p$-Laplacian}. J. Math. Anal. Appl.,
\textbf{338} (\textbf{2}) (2008), 784--792.

\bibitem{horsim2009on}
T. L. Horv\'at, P. L. Simon;
\emph{On the exact number of solutions of a singular boundary-value problem}. 
Diff. Int. Eq., \textbf{22} (\textbf{7-8}) (2009), 787--796.

\bibitem{sim2009exact}
P. L. Simon;
\emph{Exact multiplicity of positive solutions for a class of singular semilinear
equations}. Differ. Equ. and Dyn. Syst., \textbf{17} (\textbf{1-2}) (2009), 147--161.

\end{thebibliography}

\end{document}
