\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2016 (2016), No. 56, 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}
\thanks{\copyright 2016 Texas State University.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2016/56\hfil Existence of solutions]
{Existence of solutions for Dirichlet quasilinear systems including a nonlinear 
function of the derivative}

\author[S. Heidarkhani, M. Ferrara, G. A. Afrouzi, G. Caristi, S. Moradi 
\hfil EJDE-2016/56\hfilneg]
{Shapour Heidarkhani, Massimiliano Ferrara,
Ghasem A. Afrouzi, Giuseppe Caristi, Shahin Moradi}

\address{Shapour Heidarkhani \newline
Department of Mathematics, Faculty of Sciences,
Razi University,  67149 Kermanshah, Iran}
\email{s.heidarkhani@razi.ac.ir}

\address{Massimiliano Ferrara \newline
 Department of Law and Economics,
University Mediterranea of Reggio Calabria,
Via dei Bianchi, 2 - 89127 Reggio Calabria, Italy}
\email{massimiliano.ferrara@unirc.it}

\address{Ghasem A. Afrouzi \newline
Department of Mathematics,
Faculty of Mathematical Sciences,
University of Mazandaran, Babolsar, Iran}
\email{afrouzi@umz.ac.ir}

\address{Giuseppe Caristi \newline
Department of Economics, University of Messina,
via dei Verdi, 75, Messina, Italy}
\email{gcaristi@unime.it}

\address{Shahin Moradi \newline
Department of Mathematics,
Faculty of Mathematical Sciences,
University of Mazandaran, Babolsar, Iran}
\email{shahin.moradi86@yahoo.com}

\thanks{Submitted  November 17, 2015. Published February 25, 2016.}
\subjclass[2010]{34B15, 58E05}
\keywords{Classical solution; Dirichlet quasilinear system; critical point theory;
\hfill\break\indent variational methods}

\begin{abstract}
 In this article we establish the existence of at least one
 non-trivial classical solution for Dirichlet quasilinear systems
 with a nonlinear dependence on the derivative.
 We use variational methods for smooth functionals defined 
 on reflexive Banach spaces, and  assumptions on the asymptotic behaviour 
 of the nonlinear data.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{remark}[theorem]{Remark}
\newtheorem{example}[theorem]{Example}
\allowdisplaybreaks


\section{Introduction}

In this article we consider the  quasilinear system
\begin{equation}\label{1}
\begin{gathered}
-(p_i-1)|u_i'(x)|^{p_i-2}u_i''(x)
=\lambda F_{u_i}(x,u_1,\dots,u_n)h_i(x,u_i'),\quad x\in(a,b),\\
u_i(a)=u_i(b)=0
\end{gathered}
\end{equation}
for $1\leq i\leq n$ where $p_i>1$ for $1\leq i\leq n$,
$\lambda>0$,  $a,b\in\mathbb{R}$ with $a<b$,
$F:[a,b]\times \mathbb{R}^{n}\to \mathbb{R}$ is measurable with respect to
$x$, for every $(t_1,\dots,t_n)\in\mathbb{R}^n$,
continuously differentiable in $(t_1,\dots,t_n)$, for almost every
$x\in [a,b]$, and $F_{t_i}(x,0,\dots,0)=0$ for all $x\in [a,b]$ and for
$1\leq i\leq n$, $h_i:[a,b]\times \mathbb{R}\to [0,+\infty[$ is a bounded and
continuous function with $m_i:=\inf_{(x,t)\in [a,b]\times
\mathbb{R}} h_i(x,t)>0$ for $1\leq i\leq n$. Here, $F_{t_i}$
denotes the partial derivative of $F$ with respect to $t_i$.

 Owing to the importance of second-order differential equations
with nonlinear derivative
dependence in physics, in the last decade or so, many authors
applied the variational method to study the existence of solutions
of the problems of the form of \eqref{1} or its variations, see,
for example, \cite{AH,AB,DHS,GHK,HH,HM,L}. We note that the main
tools in these cited papers are critical points theorems due to
Ricceri and their variants.

Our goal is to establish some new criteria for system \eqref{1}
to have at least one non-trivial classical solution by
applying the following critical points theorem due to Ricceri
\cite[Theorem 2.1]{R1}.

\begin{theorem}\label{th21}
Let $X$ be a reflexive real Banach
space, let $\Phi, \Psi:X \to \mathbb{R}$ be two
G\^ateaux differentiable functionals such that $\Phi$ is
sequentially weakly lower semicontinuous, strongly continuous and
coercive in $X$ and $\Psi$ is sequentially weakly upper
semicontinuous in $X$. Let $I_\lambda$ be the functional defined
as $I_\lambda:=\Phi-\lambda\Psi$, $\lambda\in\mathbb{R}$, and for
every $r>\inf_X \Phi$, let $\varphi$ be the function defined as
$$
\varphi(r):=\inf_{u\in\Phi^{-1}(-\infty,r)}
\frac{\sup_{v\in\Phi^{-1}(-\infty,r)}\Psi(v)
-\Psi(u)}{r-\Phi(u)}.
$$
Then, for every
$r > \inf_X \Phi$ and every $\lambda\in(0,\frac{1}{\varphi(r)})$, the
restriction of the functional $I_\lambda$ to
 $\Phi^{-1}(-\infty,r)$ admits a global minimum, which is a
critical point (precisely a local minimum) of $I_\lambda$ in $X$.
\end{theorem}

 We refer the interested reader to the papers
\cite{AHB,FeMol,GM,MolSer1,MolSer2,MolRad} in
 which Theorem \ref{th21} has been used for the existence of at least
one nontrivial solution for boundary value problem.

\section{Preliminaries}

Let $X$ be the Cartesian product of $n$ Sobolev
spaces $W^{1,p_1}_0([a,b]),\dots$,
and $W^{1,p_n}_0([a,b])$, i.e.,
$X=W^{1,p_1}_0([a,b])\times\cdots\times W^{1,p_n}_0([a,b])$,
equipped with the norm
$$
\|(u_1,\dots,u_n)\| =\sum_{i=1}^{n}\|u'_i\|_{p_i}
$$
where
$$
\|u'_i\|_{p_i} = \Big(\int_{a}^{b}|u'_i(x)|^{p_i}dx\Big)^ {1/p_i}, \quad i= 1,
\dots, n.
$$
Since $p_i > 1$ for $i = 1,\dots, n$, $X$ is
compactly embedded in $(C([a,b]) )^n$.

 By a classical solution of \eqref{1}, we mean a function
$u=(u_1,\dots,u_n)$ such
that, for $i =1,\dots,n$, $u_i\in C^1[a,b]$, $u'_i\in AC[a,b]$,
and $u_i(x)$ satisfies \eqref{1} a.e. on $[a,b]$. We say that a
function $u=(u_1,\dots,u_n)\in X$ is a weak solution of the
system \eqref{1} if
\begin{align*}
&\sum_{i=1}^{n}\int_{a}^{b}\Big(\int_0^{u_i'(x)}
\frac{(p_i-1)|\tau|^{p_i-2}}{h_i(x,\tau)}d\tau\Big)v_i'(x)dx\\
&-\lambda\int_{a}^{b}\sum_{i=1}^{n}F_{u_i}(x,u_1(x),\dots,u_n(x))v_i(x)dx=0
\end{align*}
for every $v=(v_1,\dots,v_n)\in X$. Using standard methods, we
see that a weak solution of system \eqref{1} is indeed a
classical solution (see \cite[Lemma 2.2]{GHK}).

Let
\begin{gather*}
\underline{p}:=\min\{p_i: 1\leq i\leq n\},\quad
\overline{p}:=\max\{p_i: 1\leq i\leq n\},\\
m_i:=\inf_{(x,t)\in [a,b]\times \mathbb{R}} h_i(x,t)>0, \quad
 \text{for } 1\leq i\leq n,\\
M_i:=\sup_{(x,t)\in [a,b]\times \mathbb{R}} h_i(x,t), \quad \
 text{for } 1\leq i\leq n,\\
\overline{M}:=\max\{M_i: 1\leq i\leq n\},\quad
\underline{M} :=\min\{m_i: 1\leq i\leq n\}.
\end{gather*}
Then, $\overline{M}\geq M_i\geq m_i \geq \underline{M}>0$ for each
$i = 1,\dots, n$.
Put
$$
H_i(x,t)=\int_0^{t}\Big(\int_0^{\tau}\frac{(p_i-1)|
\delta|^{p_i-2}}{h_i(x,\delta)}d\delta\Big)d\tau
$$
 for all $(x,t)\in [a,b]\times \mathbb{R}$, $1\leq i \leq n$.

\section{Main results}

We formulate our min result as follows.
Put
$$
p^{\ast}=\begin{cases}
\overline{p} & \text{if } b-a\geq 1,\\
\underline{p}& \text{if } 0<b-a< 1.
\end{cases}
$$

\begin{theorem}\label{t3}
Assume that
\begin{align}\label{new eq7}
\sup_{r>0}\frac{r}{\int_{a}^{b}\max_{(t_1,\dots,t_n)\in
Q(\frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}r}{2^{\underline{p}}})}
F(x,t_1,\dots,t_n)dx}>1
\end{align}
where
\[
Q(\frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}r}{2^{\underline{p}}})
=\Big\{(t_1,\dots,t_n)\in
\mathbb{R}^n:\ \sum_{i=1}^{n}\frac{| t_i|^{p_i}}{p_i}\leq
\frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}r}{2^{\underline{p}}}\Big\}
\]
and there are a non-empty open set $D\subseteq(a,b)$
and $B\subset D$ of positive Lebesgue measure such that
\begin{gather*}
\limsup_{(\xi_1,\dots,\xi_n)\to(0^+,\dots,0^+)}\frac{\operatorname{ess\,inf}
_{x\in B}F(x,\xi_1,\dots,\xi_n)}{\sum_{i=1}^{n}
 |\xi_i|^{\underline{p}}}=+\infty,\\
\liminf_{(\xi_1,\dots,\xi_n)\to(0^+,\dots,0^+)}
\frac{\operatorname{ess\,inf}_{x\in D}F(x,\xi_1,\dots,\xi_n)}{\sum_{i=1}^{n}
 |\xi_i|^{\underline{p}}}>-\infty.
\end{gather*}
Then, for each
$$
\lambda\in \Lambda
=\Big(0,\sup_{r>0}\frac{r}{\int_{a}^{b}\max_{(t_1,\dots,t_n)\in
Q(\frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}r}{2^{\underline{p}}})
}F(x,t_1,\dots,t_n)dx}\Big),
$$
system \eqref{1} admits at
least one non-trivial classical solution
 $u_\lambda=(u_{1\lambda},\dots,u_{n\lambda})$ in $X$.
Moreover, we have
$$
\lim_{\lambda\to0^+}\|u_\lambda\|=0
$$
and the real function
\[
\lambda\to \sum_{i=1}^{n}\int_{a}^{b}H_i(x,u_{i\lambda}'(x))dx
 -\lambda \int_{a}^{b}F(x,u_{1\lambda}(x),\dots,u_{n\lambda}(x))dx
\]
is negative and strictly decreasing in the open interval $\Lambda$.
\end{theorem}

\begin{proof}
Our aim is to apply Theorem \ref{th21} to system \eqref{1}.
Let the functionals
$\Phi, \Psi$ be defined by
\begin{gather}\label{e8}
\Phi(u)=\sum_{i=1}^{n}\int_{a}^{b}H_i(x,u_i'(x))dx ,\\
\label{e9}
\Psi(u)=\int_{a}^{b}F(x,u_1(x),\dots,u_n(x))dx
\end{gather}
for $u=(u_1,\dots,u_n)\in X$, and put
$$
I_\lambda(u)= \Phi(u)-\lambda\Psi(u)
$$
for $u\in X$. Let us prove that the functionals $\Phi$ and $\Psi$ satisfy the
required conditions in Theorem \ref{th21}. It is well known
that $\Psi$ is a differentiable functional whose differential at
the point $u\in X$ is
$$
\Psi'(u)(v)=\int_{a}^{b}\sum_{i=1}^{n}F_{u_i}(x,u_1(x),\dots,u_n(x))v_i(x)dx
$$
for every $v=(v_1,\dots,v_n)\in X$,
and $\Psi$ is sequentially weakly upper semicontinuous. Moreover,
according to the definition $\Phi$, we see that $\Phi$ is
continuous. Since
 $0<\underline{M}\leq h_i(x,t)\leq  \overline{M}$ for each $(x,t)\in
[a,b]\times \mathbb{R}$ and $i=1,\dots,n$, from \eqref{e8} we see
that
\begin{equation}\label{6}
\frac{1}{\overline{M}}
\sum_{i=1}^{n}\frac{\|u'_i\|^{p_i}_{p_i}}{p_i}
\leq  \Phi(u)\leq\frac{1}{\underline{M}}
 \sum_{i=1}^{n}\frac{\|u'_i\|^{p_i}_{p_i}}{p_i}\
\quad \text{for all}\ u\in X.
\end{equation}
From \eqref{6}, it follows $\lim_{\|u\|\to+\infty}\Phi(u)=+\infty$,
namely $\Phi$ is coercive. Moreover, $\Phi$ is continuously differentiable
 whose differential at the point $u\in X$ is
\[
\Phi'(u)(v)= \sum_{i=1}^{n}\int_{a}^{b}\Big(\int_0^{u_i'(x)}
\frac{(p_i-1)|\tau|^{p_i-2}}{h_i(x,\tau)}d\tau\Big)v_i'(x)dx
\]
for every $v\in X$. Furthermore, $\Phi$ is
sequentially weakly lower semicontinuous (see
\cite[Proposition25.20]{Zid}).
 Since for $1 \leq i\leq n$,
$$
\max_{x\in [a,b]}|u_i(x)|\leq
\frac{(b-a)^\frac{p_i-1}{p_i}}{2}\|u'_i\|_{p_i}
$$
for each $u_i\in W_0^{1,p_i}([a,b])$ (see \cite{T}), we have
\begin{equation*}
\max_{x\in [a,b]}\sum_{i=1}^{n}\frac{| u_i(x)|^{p_i}}{p_i}
\leq \frac{(b-a)^{p^{\ast}-1}}{2^{\underline{p}}}
 \sum_{i=1}^{n}\frac{\|u'_i\|^{p_i}_{p_i}}{p_i}
\end{equation*}
for each $u\in X$. This, for each $r>0$, along with \eqref{6}, ensures
that
\begin{equation} \label{7}
\begin{aligned}
 &\Phi^{-1}(-\infty,r)\\
&=\{u\in X;\Phi(u)< r\} \\
&\subseteq \Big\{ u\in X: \max\sum_{i=1}^{n}\frac{|u_i(x)|^{p_i}}{p_i}
\leq  \frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}r}{2^{\underline{p}}}
 \text{ for each }  x\in [a,b]\Big\}.
\end{aligned}
\end{equation}
 So,
\begin{equation}\label{8}
\max_{u\in\Phi^{-1}(-\infty,r)}\Psi(u)
\leq \int_{a}^{b}\max_{(t_1,\dots,t_n)\in
Q(\frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}r}{2^{\underline{p}}})}
F(x,t_1,\dots,t_n)dx.
\end{equation}
From the definition of $\varphi(r)$, since
$0\in \Phi^{-1}(-\infty, r)$ and $\Phi(0)=\Psi(0)=0$, one has
\begin{equation}\label{new eq11}
\begin{aligned}
\varphi(r)
&= \inf_{u\in\Phi^{-1}(-\infty,r)}
\frac{(\sup_{v\in\Phi^{-1}(-\infty,r)}\Psi(v))-\Psi(u)}{r-\Phi(u)}\\
&\leq \frac{\sup_{v\in\Phi^{-1}(-\infty,r)}\Psi(v)}{r}\\
& \leq\frac{\int_{a}^{b}\max_{(t_1,\dots,t_n)\in
Q(\frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}r}{2^{\underline{p}}})
} F(x,t_1,\dots,t_n)dx}{r}.
\end{aligned}
\end{equation}
Hence, putting
$$
\lambda^*=
\sup_{r>0}\frac{r}{\int_{a}^{b}\max_{(t_1,\dots,t_n)\in
Q(\frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}r}{2^{\underline{p}}})
}F(x,t_1,\dots,t_n)dx},
$$
Theorem \ref{th21} ensures that for every
$\lambda \in (0, \lambda^*)\subseteq(0,\frac{1}{\varphi(r)})$,
the functional $I_\lambda$ admits at least
one critical point (local minima)
$u_\lambda\in \Phi^{-1}(-\infty,r)$.
Now for every fixed $\lambda \in (0,\lambda^*)$ we show that
$u_\lambda=(u_{1\lambda},\dots,u_{n\lambda})\neq0$ and the map
$$
(0, \lambda^*)\ni \lambda \mapsto I_\lambda(u_\lambda),
$$
is negative. To this end, let us verify that
\begin{equation}\label{e10}
\limsup_{\|u\|\to 0^+}\frac{\Psi(u)}{\Phi(u)}=+\infty.
\end{equation}
Owing to our assumptions at zero, we can fix a sequence
$\{\xi^k=(\zeta^k,\dots,\zeta^k) \}\subset(\mathbb{R}^{+})^{n}$
converging to zero and two constants $\sigma,\ \kappa$ (with
$\sigma>0$) such that
\begin{gather*}
\lim_{k\to +\infty}\frac{\operatorname{ess\,inf}_{x\in B}F(x,\xi^k)}
{\sum_{i=1}^{n}|\zeta^k|^{\underline{p}}}=+\infty,\\
\operatorname{ess\,inf}_{x\in D}F(x, \xi)\geq\kappa
\sum_{i=1}^{n}|\xi_i|^{\underline{p}}
\end{gather*}
where
$\xi=(\xi_1,\dots,\xi_n)$, $\xi_i\in [0, \sigma]$,
 $1\leq i \leq n$.
Now, fix a set $C\subset B$ of positive measure and a function
$v=(v_1,\dots,v_n)\in X$ such that:
\begin{itemize}
\item[(i)] $v_i(x)\in [0, 1]$, for $1\leq i \leq n$ and for every $x\in[a,b]$,
\item[(ii)] $v(x)=(1,\dots,1)$ for every $x\in C$,
\item[(iii)] $v(x)=(0,\dots,0)$ for every $x\in(a,b)\setminus D$.
\end{itemize}
Hence, fix $M>0$ and consider a real positive number $\eta$ with
$$
M<\frac{n\eta \operatorname{meas}(C) +\kappa\int_{D \setminus C
}\sum_{i=1}^{n}|v_i(x)|^{\underline{p}} dx}
{\frac{1}{\underline{p}\underline{M}}\sum_{i=1}^{n}
\|v'_i\|^{\underline{p}}_{\underline{p}}}.
$$
Then, there is $k_0\in \mathbb{N}$ such that $\zeta^{k}<\sigma$ and
$$
\operatorname{ess\,inf}_{x\in B}F(x, \xi^{k})
\geq\eta\sum_{i=1}^{n}|\zeta^{k}|^{\underline{p}}
$$
for every $k> k_0$. Now, for every $k>k_0$, recalling the
properties of the function $v$ (that is
$0\leq\zeta^k v_i(x)<\sigma$, $1\leq i \leq n$ for $k$ sufficiently large),
one has
\begin{align*}
\frac{\Psi(\xi^{k} v)}{\Phi(\xi^{k}v)}
&=\frac{\int_CF(x,\zeta^k,\dots,\zeta^k)dx+\int_{D \setminus C
}F(x,\zeta^k v_1(x),\dots,\zeta^kv_n(x))dx}{\Phi(\xi^{k} v)}
\\
&>\frac{n\eta \operatorname{meas}(C)
+\kappa\int_{D \setminus C }\sum_{i=1}^{n}|v_i(x)|^{\underline{p}} dx}
{\frac{1}{\underline{p}\underline{M}}\sum_{i=1}^{n}
\|v'_i\|^{\underline{p}}_{\underline{p}}}>M
\end{align*}
where $\xi^{k} v=(\zeta^k v_1,\dots,\zeta^k v_n)$. Since $M$
could be consider arbitrarily large, it follows that
$$
\lim_{k\to\infty}\frac{\Psi(\xi^{k} v)}{\Phi(\xi^{k} v)}=+\infty,
$$
from which \eqref{e10} follows.
 Hence, there exists a
sequence $\{w_{k}=(w_{1k},\dots,w_{nk})\}\subset X$ strongly converging
 to zero, $w_k\in\Phi^{-1}(-\infty, r)$ and
$$
I_\lambda(w_{k})=\Phi(w_{k})-\lambda\Psi(w_{k})<0.
$$
Since $u_\lambda=(u_{1\lambda},\dots,u_{n\lambda})$ is a global
minimum of the restriction of $I_\lambda$ to the set $\Phi^{-1}(-\infty,r)$,
we deduce that
\begin{equation}\label{e11}
I_\lambda(u_\lambda)<0,
\end{equation}
so that $u_\lambda$ is not trivial. From \eqref{e11} we easily
see that the map
\begin{equation}\label{e12}
(0, \lambda^*)\ni \lambda \mapsto I_\lambda(u_\lambda)
\end{equation}
is negative.

Now we prove that $\lim_{\lambda\to 0^+}\|u_\lambda\|=0$.
 Since $\Phi$ is coercive and for
$\lambda\in(0,\lambda^*)$ the solution $u_\lambda\in
\Phi^{-1}(-\infty, r)$, one has that there exists a positive
constant $L$ such that $\|u_\lambda\|\leq L$ for every
$\lambda\in(0,\lambda^*)$. Therefore, there exists a positive
constant $\mathcal{M}$ such that
\begin{equation}\label{e13}
\Big|\int_a^b\sum_{i=1}^{n}F_{u_{\lambda_i}}
(x,u_{1\lambda}(x),\dots,u_{n\lambda}(x))u_\lambda(x)dx\Big|
\leq\mathcal{M}\|u_\lambda\|\leq\mathcal{M}L
\end{equation}
for every $\lambda\in(0,\lambda^*)$. Since $u_\lambda$ is a
critical point of $I_\lambda$, we have
$I'_\lambda(u_\lambda)(v)=0$ for any $v\in X$ and every
$\lambda\in(0,\lambda^*)$. In particular
$I'_\lambda(u_\lambda)(u_\lambda)=0$;  that is,
\begin{equation}\label{e14}
\Phi'(u_\lambda)(u_\lambda)
=\lambda\int_a^b\sum_{i=1}^{n}F_{u_{\lambda_i}}(x,u_{1\lambda}(x),
\dots,u_{n\lambda}(x))u_\lambda(x)dx
\end{equation}
for every $\lambda\in(0,\lambda^*)$. Then, since
$$ 
0\leq \frac{1}{\overline{M}}\sum_{i=1}^{n}\|u'_i\|^{p_i}_{p_i} 
\leq \Phi'(u_\lambda)(u_\lambda)
$$
by \eqref{e14} it follows that
$$
0\leq \frac{1}{\overline{M}}\sum_{i=1}^{n}\|u'_i\|^{p_i}_{p_i} 
\leq \lambda\int_a^b\sum_{i=1}^{n}F_{u_{\lambda_i}}
(x,u_{1\lambda}(x),\dots,u_{n\lambda}(x))u_\lambda(x)dx
$$ 
for any $\lambda\in(0,\lambda^*)$. Letting $\lambda\to 0^+$, by
\eqref{e13} we have $\lim_{\lambda\to0^+}\|u_\lambda\|=0$.
Further, taking $X\hookrightarrow (C([a,b]))^n$ into account, one
has
\begin{equation}\label{e15}
\lim_{\lambda\to0^+}\|u_\lambda\|_\infty=0.
\end{equation}
Finally, we show that the map $\lambda\mapsto
I_\lambda(u_\lambda)$ is strictly decreasing in $(0, \lambda^*)$.
For our aim we see that for any $u\in X$,
\begin{equation}\label{e16}
I_\lambda(u)=\lambda\Big(\frac{\Phi(u)}{\lambda}-\Psi(u)\Big).
\end{equation}
Now, let us fix $0<\lambda_1<\lambda_2<\lambda^*$ and let
$\bar{u}_{\lambda_ 1}=(u_{1\lambda_1},\dots,u_{n\lambda_1})$
and $\bar{u}_{\lambda_2}=(u_{1\lambda_{2}},\dots,u_{n\lambda_{2}})$ 
be the global minimum of the functional $I_{\lambda_ i}$ restricted to 
$\Phi(-\infty, r)$ for $i=1,2$. Also, let
$$
m_{\lambda_ i}=\Big(\frac{\Phi(\bar{u}_{\lambda_ i})}
{\lambda_ i}-\Psi(\bar{u}_{\lambda_ i})\Big)
=\inf_{v\in \Phi^{-1}(-\infty, r)}\Big(\frac{\Phi(v)}{\lambda_i}-\Psi(v)\Big)
$$
for $i=1,2$.
Clearly, \eqref{e12} and \eqref{e16}, since $\lambda>0$,
imply that
\begin{equation}\label{e17}
m_{\lambda i}<0, \quad \text{for } i=1 \ 2.
\end{equation}
Moreover, since  $0<\lambda_1<\lambda_2$, we have
\begin{equation}\label{e18}
m_{\lambda_ 2}\leq m_{\lambda_ 1}.
\end{equation}
Hence, from \eqref{e16}-\eqref{e18} and again since
$0<\lambda_1<\lambda_2$, we obtain
$$
I_{\lambda _2} (\bar{u}_{\lambda_ 2})=\lambda_2 m_{\lambda_ 2}
\leq \lambda_2 m_{\lambda_ 1}
<\lambda_1 m_{\lambda_ 1}=I_{\lambda _1 }(\bar{u}_{\lambda_ 1}),
$$
which means the map $ \lambda \mapsto I_\lambda(u_\lambda)$ is
strictly decreasing in $\lambda \in (0, \lambda^*)$. Since
$\lambda<\lambda^*$ is arbitrary, we observe $ \lambda \mapsto
I_\lambda(u_\lambda)$ is strictly decreasing in $(0, \lambda^*)$.
The proof is complete.
\end{proof}

 \begin{remark} \label{rmk3.1} \rm
 Here employing Ricceri's variational principle we are looking for
the existence of critical points of the functional $I_\lambda$
naturally associated to  system \eqref{1}.
 We emphasize that
by direct minimization, we can not argue, in general for finding
the critical points of $I_\lambda$. Because, in general,
$I_\lambda$ can be unbounded from the following in $X$.
 Indeed, for example, in the case when 
\[
F(x,t_1,\dots,t_n)=(\sum_{i=1}^{n}(|t_i|+|t_i|^{q_i}))
\]
 for all $(x,t_1,\dots,t_n)\in [a,b]\times \mathbb{R}^{n}$ with 
$q_i\in(\overline{p},+\infty)$ for 
$1\leq i\leq n$, for any fixed $u\in X\backslash\{0\} $ and 
$\iota\in \mathbb{R}$, we obtain
\begin{align*}
I_\lambda(\iota u)
&=\Phi(\iota u)-\lambda\int_a^bF(x,\iota u_1(x),\dots,\iota u_n(x))dx\\
&\leq\frac{\iota^{\overline{p}}}
{\underline{p}\underline{M}}\sum_{i=1}^{n}\|u'_i\|^{p_i}_{p_i}-\lambda\iota
\sum_{i=1}^{n}\|u_i\|_{L^1}-\lambda
\sum_{i=1}^{n}\iota^{q_i}\|u_i\|_{L^{q_i}}^{q_i}
\to-\infty
\end{align*}
as $\iota\to+\infty$.
\end{remark}

\begin{remark} \label{rmk3.2} \rm
 We want to point out that the energy functional $I_\lambda$ associated with 
 system \eqref{1} is not coercive. Indeed, fix $u\in X\backslash\{0\} $
 and $\iota\in \mathbb{R}$, then, 
for $F(x,t_1,\dots,t_n)=\sum_{i=1}^{n}|t_i|^{s_i}$ 
for all $(x,t_1,\dots,t_n)\in [a,b]\times
\mathbb{R}^{n}$ with $s_i\in(\overline{p},+\infty)$. We
have
\begin{align*}
I_\lambda(\iota u)&=\Phi(\iota u)-\lambda\int_a^bF(x,\iota
u_1(x),\dots,\iota u_n(x))dx\\
&\leq\frac{\iota^{\overline{p}}}{\underline{p}
\underline{M}}\sum_{i=1}^{n}\|u_i\|^{p_i}_{p_i}-\lambda\sum_{i=1}^{n}\iota^{s_i}
\|u_i\|_{L^{s_i}}^{s_i}\to-\infty
\end{align*}
as $\iota\to+\infty$.
 \end{remark}

 \begin{remark} \label{rmk3.3} \rm
 If in Theorem \ref{t3}, $F_{t_i}(x,t_1,\dots,t_n)\geq0$ 
 for a.e. $x\in [a, b]$, $1\leq i\leq n$, the condition \eqref{new eq7} 
becomes to the more simple and significative form
 \begin{equation}\label{new eq8}
\sup_{r>0}\frac{r}{\int_{a}^{b}
F(x,\sqrt[\underline{p}]{\frac{(b-a)^{p^{\ast}-1}
\overline{p}^2\overline{M}r}{2^{\underline{p}}}},
 \dots,\sqrt[\underline{p}]{\frac{(b-a)^{p^{\ast}-1}\overline{p}^2
\overline{M}r}{2^{\underline{p}}}})dx}>1.
 \end{equation}
 Moreover, if 
 \[
 \limsup_{r\to+\infty}\frac{r}{\int_{a}^{b} 
F(x,\sqrt[\underline{p}]{\frac{(b-a)^{p^{\ast}-1}\overline{p}^2
\overline{M}r}{2^{\underline{p}}}},
 \dots,\sqrt[\underline{p}]{\frac{(b-a)^{p^{\ast}-1}\overline{p}^2
\overline{M}r}{2^{\underline{p}}}})dx}>1,
\]
 then,  \eqref{new eq8} automatically holds.
 \end{remark}

\begin{remark} \label{rmk3.4} \rm
For fixed $\bar{r}>0$ let
$$
\frac{\bar{r}}{\int_{a}^{b}\max_{(t_1,\dots,t_n)\in
Q(\frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}\bar{r}}{2^{\underline{p}}})}
F(x,t_1,\dots,t_n)dx}>1.
$$ 
Then the result of Theorem \ref{t3} holds. 
In fact, it ensures the existence of least one non-trivial
classical solution
 $u_\lambda=(u_{1\lambda},\dots,u_{n\lambda})$ in $X$
such that
$$
\max\sum_{i=1}^{n}\frac{|u_{i\lambda}(x)|^{p_i}}{p_i}\leq
 \frac{(b-a)^{p^{\ast}-1}\overline{p}\overline{M}\bar{r}}
{2^{\underline{p}}}\quad \text{for each }  x\in  [a,b].
$$
 \end{remark}

 \begin{remark} \label{rmk3.5} \rm
 We observe that Theorem \ref{t3} is a bifurcation result in the sense that the pair
 $(0,0)$ belongs to the closure of the set
 $\{ (u_\lambda, \lambda)\in X\times (0, +\infty): u_\lambda=(u_{1\lambda},
\dots,u_{n\lambda})$ is a non-trivial classical solution of  \eqref{1} $\}$
 in $X\times \mathbb{R}$.
Indeed, by Theorem \ref{t3} we have that
 $$
\|u_\lambda\|\to 0 \quad \text{as } \lambda \to 0.
$$
 Hence, there exist two sequences $\{u_j=(u_{1j}\dots,u_{nj})\}$ in $X$ 
and $\{(\lambda_1,\dots,\lambda_n)\}$ in $(\mathbb{R}^+)^n$
 (here $u_{ij}=u_{\lambda_i}$) such that for $1\leq i\leq n$
 $$ 
\lambda_i \to 0^+\quad \text{and}\quad\|u_{ij}\|\to 0,
$$
 as $j\to+\infty$. Moreover, we want to emphasis that because the map
$(0, \lambda^*)\ni \lambda \mapsto I_\lambda(u_\lambda)$
is strictly decreasing, for every
 $\lambda_1,\lambda_2\in(0,\lambda^\star)$, with
 $\lambda_1\neq\lambda_2$, the solutions $\bar{u}_{\lambda_1}$ and
 $\bar{u}_{\lambda_2}$ ensured by Theorem \ref{t3} are different.
 \end{remark}

 We present the following example in which the hypotheses of Theorem
\ref{t3} are satisfied.

\begin{example} \label{examp3.1} \rm 
Consider the system
\begin{equation}\label{e29}
\begin{gathered}
-u''_1(x)=\lambda F_{u_1}(u_1,u_2)h_1(u_1'), \quad x\in(0,1),\\
-u''_2(x)=\lambda F_{u_2}(u_1,u_2)h_2(u_2'),\quad x\in(0,1),\\
u_1(0)=u_1(1)=0,\quad u_2(0)=u_2(1)=0
\end{gathered}
\end{equation}
where $F(t_1,t_2)=(t_1^2+t_2^2)\sin^2(t_1)+\cos(t_2)$ for $t_1,t_2\in\mathbb{R}$,
$$
h_1(t_1)=  \begin{cases}
2, &\text{if } t_1\in(-\infty,0),\\
\frac {2}{1+\sqrt{t_1}}, &\text{if }  t_1\in[0,1],\\
1,  &\text{if }  t_1\in(1,+\infty)
\end{cases} 
$$
 and $h_2(t_{2})=\frac{4}{3+sint_{2}}$ for $t_{2}\in \mathbb{R}$. 
We observe that $F$, $h_1$ and $h_2$ are continuous, $h_1$ and $h_2$ 
are bounded with  $m_1=\inf h_1=1$, $m_2=\inf h_2=1$, $M_1=\sup h_1=2$ 
and $M_2=\sup h_2=2$.
 By the definitions of $h_1$ and $h_2$ we have
 $$ 
H_1(t_1)= \begin{cases}
 \frac{1}{4}t_1^2, &\text{if } t_1\in(-\infty,0),\\
 \frac{1}{4}t_1^2+\frac{2}{15}t_1^{\frac{5}{2}}, &\text{if }  t_1\in[0,1],\\
 \frac{1}{2}t_1^2-\frac{1}{4}t_1+\frac{2}{15},  &\text{if }  t_1\in(1,+\infty)
 \end{cases} 
$$
and $H_2(t_{2})=\frac{3}{8}t_{2}^2-\frac{\sin t_{2}}{4}$ for
$t_{2}\in \mathbb{R}$. It is easy to see that all assumptions of
Theorem \ref{t3} are satisfied. By Theorem \ref{t3} for each
$\lambda\in (0,\frac{1}{4 }$, system \eqref{e29} admits at
least one nontrivial classical solution
$u_\lambda=(u_{1\lambda},u_{2\lambda})$ in $X$. Moreover, we have
$\lim_{\lambda\to 0^+}\|u_\lambda\|=0$
and the real function
\begin{align*}
\lambda\mapsto &\int_0^{1}\Big(H_1(u_{1\lambda}'(x))+H_{2}
 (u_{2\lambda}'(x))\Big)dx\\
&-\lambda \int_0^{1}\big((u^2_{1\lambda}(x)+u^2_{2\lambda}(x))
 \sin^2(u_{1\lambda}(x))+\cos(u_{2\lambda}(x))\big)dx
\end{align*}
is negative and strictly decreasing in $(0,\frac{1}{4 })$.
\end{example}

 As an application of Theorem \ref{t3}, we consider the  problem
\begin{equation}\label{new e22}
\begin{gathered}
-(p-1)|u'(x)|^{p-2}u''(x)=\lambda f(x,u)h(u'),\quad x\in(a,b),\\
u(a)=u(b)=0
\end{gathered}
\end{equation}
where $p>1$, $\lambda>0$, $f:[a,b]\times\mathbb{R}\to\mathbb{R}$ is an
$L^1$-Carath\'{e}odory function such that $f(x,0)=0$ for a.e. $x\in [a, b]$, 
and $h:\mathbb{R}\to [0,+\infty)$ is a
bounded and continuous function with
$m=\inf_{t\in\mathbb{R}}h(t)>0$.

Let $M=\sup_{t\in\mathbb{R}}h(t)$,
\begin{gather*}
F(x,t)=\int_0^{t}f(x,\xi)d\xi \quad \text{for every } (x,t)\in[a,b]\times\mathbb{R},\\
H(t)=\int_0^{t}\Big(\int_0^{\tau}\frac{(p-1)|\delta|^{p-2}}{h(\delta)}d\delta\Big)d\tau
\end{gather*}
 for all $t\in \mathbb{R}$.
For $\gamma>0$, define
$$
W(\gamma)=\{t\in\mathbb{R}:|t|^p\leq p\gamma\}.
$$
 We conclude this section by giving the following consequence of Theorem \ref{t3}.

 \begin{theorem}\label{t34}
Assume that
 \[
 \sup_{r>0}\frac{r}{\int_{a}^{b}\max_{t\in W(\frac{(b-a)^{p-1}pMr}{2^p})}
 F(x,t)dx}>1
 \]
 and there are a non-empty open sets $D\subseteq(0, T)$ and $B\subset D$ 
of positive Lebesgue measure such that 
\begin{gather*}
\limsup_{\xi\to0^+}\frac{\operatorname{ess\,inf}_{x\in B}F(x,\xi)}{\xi^{p}}=+\infty,\\
\liminf_{\xi\to0^+}\frac{\operatorname{ess\,inf}_{x\in D}F(x,\xi)}{\xi^{p}}>-\infty.
\end{gather*}
Then, for each 
\[
\lambda\in\Lambda=
\Big(0,\sup_{r>0}\frac{r}{\int_{a}^{b}\max_{t\in
W(\frac{(b-a)^{p-1}pMr}{2^p})} F(x,t)dx}\Big),
\]
  problem
\eqref{new e22} admits at least one nontrivial classical solution
$u_\lambda\in W^{1,p}_0([a,b])$ such that
$\lim_{\lambda\to0^+}\|u_\lambda\|=0$ and the real function
\[
\lambda\mapsto \int_{a}^{b}H(u'_{\lambda}(x))dx-\lambda \int_{a}^b
F(x,u_{\lambda}(x))dx
\]
is negative and strictly decreasing in the open interval
$\Lambda$.
 \end{theorem}

  Now, we point out the following result, as a consequence of Theorem 
\ref{t34} in which the function $f$ has
separated variables.

\begin{theorem}\label{th33}
Let $\alpha\in L^\infty([a,b])$ with
$\operatorname{ess\,inf}_{x\in[a,b]}\alpha(x)>0$.
Further, let $f:\mathbb{R}\to\mathbb{R}$ be a continuous function
such that $f(0)=0$. Put $F(\xi)=\int_0^{\xi}f(t)dt$ for all
$\xi\in\mathbb{R}$, and assume that   
\begin{gather*}
\frac{1}{\int_{a}^{b}\alpha(x)dx}\sup_{r>0}\frac{r}{\max_{\xi\in
W(\frac{(b-a)^{2}Mr}{2})} F(\xi)}>1, \\
\lim_{\xi\to0^+}\frac{F(\xi)}{\xi^2}=+\infty.
\end{gather*}
Then, for each
\[
\lambda\in\Lambda=\Big(0,\frac{1}{\int_{a}^{b}\alpha(x)dx}
\sup_{r>0}\frac{r}{\max_{\xi\in
W(\frac{(b-a)^{2}Mr}{2})} F(\xi)}\Big),
\]
 the problem
\begin{equation}
\begin{gathered}
-u''=\lambda \alpha(x)f(u)h(u'), \quad x\in(a,b),\\
u(a)=u(b)=0
\end{gathered}
\end{equation}
admits at least one nontrivial classical solution 
$u_\lambda\in W^{1,2}_0([a,b])$ such that
$$
\lim_{\lambda\to 0^+}\|u_\lambda\|=0
$$ 
and the real function
\begin{align*}
\lambda\to \int_{a}^{b}H(u'_{\lambda}(x))dx
-\lambda \int_{a}^b\alpha(x)F(u_\lambda(x))dx
\end{align*}
is negative and strictly decreasing in $\Lambda$.
\end{theorem}

 Next we given an example to illustrate Theorem \ref{th33}.

\begin{example} \label{examp3.2} \rm
Consider the  problem
\begin{equation}\label{new e20}
\begin{gathered}
-u''=\lambda(1+e^x)f(u)h(u'),\quad x\in(0,1),\\
 u(0)=u(1)=0
 \end{gathered}
 \end{equation}
 where $f(t)=t^2\sin^2(t)-\sin(t)$ for $t\in\mathbb{R}$ and
 $$
h(t)= \begin{cases}
 1, &\text{if } t\in(-\infty,0),\\
 \frac {1}{1+t^2}, &\text{if } t\in[0,1],\\
 \frac{1}{2},  &\text{if }  t\in(1,+\infty).
 \end{cases}
$$
 We observe that $f$ and $h$ are continuous and $h$ is bounded with
 $m=\inf_{t\in \mathbb{R}} h(t)=\frac{1}{2}$ and 
$M=\sup_{t\in \mathbb{R}} h(t)=1$.
 By the expression of $f$ and $h$ we have
 $$
F(t)=\frac{1}{6}t^3-\frac{1}{4}\Big(t^2\sin(2t)+t\cos(2t)
-\frac{1}{2}\sin(2t)\Big)+\cos(t)-1
$$
 for every $t\in\mathbb{R}$ and
 $$
 H(t)=  \begin{cases}
 \frac{1}{2}t^2, &\text{if } t\in(-\infty,0),\\
 \frac{1}{2}t^2+\frac{1}{12}t^{4}, &\text{if }  t\in[0,1],\\
 t^2-\frac{1}{2}t+\frac{1}{12},  &\text{if }  t\in(1,+\infty).
  \end{cases}
$$
 It is easy to see that all assumptions of Theorem \ref{th33} are
 satisfied. By Theorem \ref{th33} for each $\lambda\in(0,\frac{1}{4e})$, 
 problem \eqref{new e20} admits at
 least one nontrivial classical solution $u_\lambda$ in $W^{1,2}_0([0,1])$.
  Moreover, $\lim_{\lambda\to 0^+}\|u_\lambda\|=0$
  and the real function
  \[
  \lambda\mapsto \int_0^{1}H(u_\lambda'(x))dx-\lambda
  \int_0^{1}(1+e^x)F(u_\lambda(x))dx
  \]
  is negative and strictly decreasing in $(0,1/(4e))$.
  \end{example}

\begin{thebibliography}{99}

\bibitem{AHB} G. A. Afrouzi, A. Hadjian, G. Molica Bisci;
  \emph{Some remarks for one-dimensional mean curvature
problems through a local minimization principle,} Adv. Nonlinear
Anal. 2 (2013) 427-441.

\bibitem{AH} G. A. Afrouzi, S. Heidarkhani;
 \emph{Three solutions for a quasilinear boundary value problem,} Nonlinear
Anal. TMA 69 (2008) 3330-3336.

\bibitem{AB} D. Averna, G. Bonanno;
 \emph{Three solutions for quasilinear two-point boundary value problem involving 
the one-dimensional $p$-Laplacian,} Proc. Edinb. Math. Soc. 47 (2004) 257-270.

\bibitem{Br}H. Br\'{e}zis;
 \emph{Analyse Functionelle- Theorie et Applications,} Masson, Paris (1983).

\bibitem{DHS} G. D'Agu\`{i}, S. Heidarkhani, A. Sciammetta;
\emph{Infinitely many solutions for a class of quasilinear two-point boundary
 value systems,} Electron. J. Qual. Theory Differ. Equ. 2015 (2015), No. 8, 1-15.

\bibitem{FeMol} M. Ferrara, G. Molica Bisci;
\emph{Existence results for elliptic problems with Hardy potential,} 
Bull. Sci. Math. 138 (2014) 846-859.

\bibitem{GM} M. Galewski, G. Molica Bisci;
 \emph{Existence results for one-dimensional fractional equations},
Math. Meth. Appl. Sci., to appear.

\bibitem{GHK} J. R. Graef, S. Heidarkhani, L. Kong;
 \emph{A critical points approach for the existence of multiple solutions 
of a Dirichlet quasilinear system,} J. Math. Anal. Appl. 388 (2012) 1268-1278.

\bibitem{HH} S. Heidarkhani, J. Henderson;
 \emph{Multiple solutions for a Dirichlet quasilinear system
containing a parameter,} Georgian Math. J. 21 (2014) 187-197.

\bibitem{HM} S. Heidarkhani, D. Motreanu;
 \emph{Multiplicity results for a two-point boundary value problem,} 
Panamer. Math. J. 19 (2009) 69-78.

\bibitem{L} R. Livrea;
 \emph{Existence of three solutions for a quasilinear two point boundary 
value problem,} Arch. Math. (Basel) 79 (2002) 288-298.

\bibitem{MolSer1} G. Molica Bisci, R. Servadei;
 \emph{A bifurcation result for non-local
fractional equations,} Anal. Appl. 13, No. 4 (2015) 371-394.

\bibitem{MolSer2} G. Molica Bisci, R. Servadei;
 \emph{Lower semicontinuity of functionals
of fractional type and applications to nonlocal equations with
critical Sobolev exponent,} Adv. Differential Equations 20 (2015)
635-660.

\bibitem{MolRad} G. Molica Bisci, V. R\u{a}dulescu;
 \emph{Bifurcation analysis of a singular elliptic problem modelling the 
equilibrium of anisotropic continuous media,} Topol. Methods Nonlinear Anal. 
45 (2015) 493-508.

\bibitem{R1} B. Ricceri;
\emph{A general variational principle and some of its
applications,} J. Comput. Appl. Math. 113 (2000) 401-410.

\bibitem{T} G. Talenti;
 \emph{Some inequalities of Sobolev type on
two-dimensional spheres, in: W. Walter (Ed.),} General
inequalities, 5 (Oberwolfach, 1986), 401-408, Internat.
Schriftenreihe Numer. Math., 80, Birkhauser, Basel (1987).

\bibitem{Zid} E. Zeidler;
 \emph{Nonlinear functional analysis and its applications,} 
Vol. II, Springer, Berlin-Heidelberg-New York (1985).

\end{thebibliography}

\end{document}
