\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2016 (2016), No. 272, pp. 1--8.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2016 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2016/272\hfil Hessian systems with convection terms]
{Existence of positive entire radial solutions to a $(k_1,k_2)$-Hessian
systems with convection terms}

\author[D.-P. Covei \hfil EJDE-2016/272\hfilneg]
{Dragos-Patru Covei}

\address{Dragos-Patru Covei \newline
Department of Applied Mathematics,
The Bucharest University of Economic Studies,
Piata Romana, 1st district, postal code 010374,
postal office 22, Romania}
\email{coveidragos@yahoo.com}

\thanks{Submitted  August 23, 2016. Published October 10, 2016.}
\subjclass[2010]{35J25, 35J47, 35J96}
\keywords{Entire solution; large solution; elliptic system}

\begin{abstract}
 In this article, we prove two new results on the existence
 of positive entire large and bounded radial solutions for nonlinear
 system with gradient terms
 \begin{gather*}
 S_{k_1}(\lambda (D^{2}u_1) )+b_1(| x| ) | \nabla u_1|^{k_1}
 =p_1(| x| ) f_1(u_1,u_2) \quad\text{for }x\in \mathbb{R}^{N},  \\
 S_{k_2}(\lambda (D^{2}u_2) ) +b_2(| x| ) | \nabla u_2|^{k_2}
 =p_2(| x| ) f_2(u_1,u_2) \quad\text{for }x\in \mathbb{R}^{N},
 \end{gather*}
 where $S_{k_i}(\lambda (D^{2}u_i) ) $ is the
 $k_i$-Hessian operator, $b_1$, $p_1$, $f_1$, $b_2$, $p_2$ and $f_2$
 are continuous functions satisfying certain properties. Our results
 expand those by Zhang and Zhou \cite{ZZ}. The main difficulty in dealing
 with our system is the presence of the convection term.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{remark}[theorem]{Remark}
\allowdisplaybreaks

\section{Introduction}

The purpose of this article is to present new results concerning the nonlinear
Hessian system with convection terms
\begin{equation}
\begin{gathered}
S_{k_1}(\lambda (D^{2}u_1) ) +b_1(| x| ) | \nabla u_1|^{k_1}
=p_1(| x| ) f_1(u_1,u_2) , x\in \mathbb{R}^{N},  \\
S_{k_2}(\lambda (D^{2}u_2) ) +b_2(| x| ) | \nabla u_2|^{k_2}
=p_2(| x| ) f_2(u_1,u_2) , x\in \mathbb{R}^{N},
\end{gathered}  \label{11}
\end{equation}
where $N\geq 3$, $b_1$, $p_1$, $f_1$, $b_2$, $p_2$, $f_2$ are
continuous functions satisfying certain properties,
$k_1,k_2\in \{1,2,\dots ,N\} $ and
\begin{equation*}
S_{k_i}(\lambda ) =\sum_{1\leq i_1<\dots <i_{k_i}\leq
N}\lambda _{i_1}\dots \lambda _{i_{k_i}},\quad
\lambda =(\lambda _1,\dots ,\lambda _{N}) \in \mathbb{R}^{N},\; i=1,2,
\end{equation*}
denotes the $k_i$-th elementary symmetric function.
In the literature $S_{k_i}(\lambda (D^{2}u_i) ) $ it is called the
$k_i$-Hessian operator. For instance, the following well known operators
are included in this class:
$S_1(\lambda (D^{2}u_i)) =\sum_{i=1}^{N}\lambda _i=\Delta u_i$
the Laplacian operator and
$S_{N}(\lambda (D^{2}u_i) ) =\prod _{i=1}^N \lambda _i=\det (D^{2}u_i) $
the Monge-Amp\`{e}re.

In recent years equations of the type \eqref{11} have been the subject of
rather deep investigations since appears from many branches of mathematics
and applied mathematics. For more surveys on these questions we refer the
paper to: Alves and Holanda \cite{CA}, Bao-Ji and Li \cite{BAOII},
 Bandle and Giarrusso \cite{BA}, Cl\'{e}ment-Man\'{a}sevich and Mitidieri \cite{CL},
De Figueiredo and Jianfu \cite{FI}, Galaktionov and V\'{a}zquez \cite{GA},
Jiang and Lv \cite{GR}, Salani \cite{SA}, Ji and Bao \cite{BAO}, Jian
\cite{J}, Peterson and Wood \cite{PW}, Pripoae \cite{PRI}, Quittner \cite{Q},
Li and Yang \cite{LY}, Li-Zhang and Zhang \cite{Z}, Viaclovsky \cite{V,VI},
Zhang and Zhou \cite{ZZ} and not the last Zhang \cite{ZEJ}.

The motivation for studying \eqref{11} comes from the work of Jiang and Lv
\cite{GR} where they study the system
\begin{gather*}
\Delta u_1+| \nabla u_1|
 =p_1(|x| ) f_1(u_1,u_2) \quad\text{for }x\in \mathbb{R}^{N}\; (N\geq 3), \\
\Delta u_2+| \nabla u_2|
 =p_2(|x| ) f_2(u_1,u_2) \quad\text{for }x\in \mathbb{R}^{N}\; (N\geq 3),
\end{gather*}
and from the recently work of Zhang and Zhou \cite{ZZ} where the authors
considered the system
\begin{gather*}
S_k(\lambda (D^{2}u_1) )
=p_1(|x| ) f_1(u_2) \quad\text{for }x\in \mathbb{R}^{N}\; (N\geq 3), \\
S_k(\lambda (D^{2}u_2) )
=p_2(|x| ) f_2(u_1) \quad\text{for }x\in \mathbb{R}^{N}\; (N\geq 3).
\end{gather*}
Our purpose is to expand and improve  the results in \cite{ZZ} for the more general
system \eqref{11}. By analogy with the work of Zhang and Zhou \cite{ZZ} we
introduce the following notations
\begin{gather*}
C_{0} =(N-1)!/[ k_1!(N-k_1)!] ,C_{00}=(N-1)!/[
k_2!(N-k_2)!] , \\
B_1^{-}(\xi ) =\frac{\xi ^{k_1-N}}{C_{0}}
\exp\Big(-\int_{0}^{\xi }\frac{1}{C_{0}}t^{k_1-1}b_1(t) dt\Big),\\
B_1^{+}(\xi ) =\xi ^{N-1}
\exp\big(\int_{0}^{\xi }\frac{1}{C_{0}}t^{k_1-1}b_1
(t) dt\Big)p_1(\xi ) , \\
B_2^{-}(\xi ) =\frac{\xi ^{k_2-N}}{C_{00}}
\exp\Big(-\int_{0}^{\xi }\frac{1}{C_{00}}t^{k_2-1}b_2(t) dt\Big),\\
B_2^{+}(\xi ) =\xi ^{N-1}\exp\Big(\int_{0}^{\xi }\frac{1}{C_{00}}
t^{k_2-1}b_2(t) dt\Big)p_2(\xi ) , \\
P_1(r) =\int_{0}^{r}\Big(B_1^{-}(r)
\int_{0}^{r}B_1^{+}(t) dt\Big) ^{1/k_1}dr, \\
P_2(r) =\int_{0}^{r}\Big(B_2^{-}(r)
\int_{0}^{r}B_2^{+}(t) dt\Big) ^{1/k_2}dr, \\
F_{1,2}(r) =\int_{a_1+a_2}^{r}\frac{1}{f_1^{1/k_1}
(t,t) +f_2^{1/k_2}(t,t) }dt\\
 \text{for }r\geq a_1+a_2>0,\;  a_1\geq 0,\; a_1\geq 0, \\
P_1(\infty ) =\lim_{r\to \infty }P_1(r) ,\quad
P_2(\infty ) =\lim_{r\to \infty }P_2(r) , \quad
F_{1,2}(\infty ) =\lim_{s\to \infty }F_{1,2}(s) .
\end{gather*}
We will always assume that the variable weights functions
$b_1,b_2,p_1,p_2$ and the nonlinearities $f_1$, $f_2$ satisfy:
\begin{itemize}
\item[(A1)] $b_1,b_2:[ 0,\infty ) \to [ 0,\infty) $ and 
$p_1,p_2:[ 0,\infty ) \to [0,\infty ) $ are spherically symmetric continuous 
functions (i.e., $p_i(x) =p_i(| x| ) $ and $b_i(x) =b_i(| x|) $ 
for $i=1,2$);

\item[(A2)]  $f_1$, $f_2:[ 0,\infty ) \times [ 0,\infty
) \to [ 0,\infty ) $ are continuous and increasing.
\end{itemize}
Here is a first result.

\begin{theorem} \label{thm1}
We assume that $F_{1,2}(\infty ) =\infty $
and \textrm{(A1)}, hold. Furthermore, if $f_1$ and $f_2$
satisfy \textrm{(A2)}, then  system \eqref{11} has at least one positive
radial solution $(u_1,u_2) \in C^{2}([ 0,\infty) ) \times C^{2}([ 0,\infty ) ) $
 with central value in $(a_1,a_2) $. Moreover, the following hold:
\begin{itemize}
\item[(1)]  If $P_1(\infty ) +P_2(\infty ) <\infty $,
then $\lim_{r\to \infty}u_1(r) <\infty $
and $\lim_{r\to \infty}u_2(r) <\infty $.

\item[(2)] If $P_1(\infty ) =\infty $ and $P_2(\infty) =\infty $, then 
\[
\lim_{r\to \infty}u_1(r) =\infty \quad\text{and}\quad
\lim_{r\to \infty}u_2(r) =\infty .
\]
\end{itemize}
\end{theorem}

In the same spirit we also have, our next result.

\begin{theorem}\label{thm2}
Assume that the hypotheses \textrm{(A1)} and \textrm{(A2)} are
satisfied.  If $P_1(\infty ) +P_2(\infty ) <F_{1,2}(\infty ) <\infty $, then
system \eqref{11} has one positive bounded radial solution 
$(u_1,u_2) \in C^{2}([ 0,\infty ) ) \times C^{2}([ 0,\infty ))$,
 with central value in $(a_1,a_2)$, such that
\begin{gather*}
a_1+f_1^{1/k_1}(a_1,a_2) P_1(r) 
\leq u_1(r) \leq F_{1,2}^{-1}(P_1(r) +P_2(r) ) , \\
a_2+f_2^{1/k_2}(a_1,a_2) P_2(r)
\leq u_2(r) \leq F_{1,2}^{-1}(P_1(r) +P_2(r) ) .
\end{gather*}
\end{theorem}

\section{Proofs of main results}

In this section we give the proof of Theorems \ref{thm1} and \ref{thm2}.
For the readers' convenience, we recall the radial form of the $k$-Hessian
operator.

\begin{remark}[see  \cite{BAO,SA}] \rm
Assume $\varphi \in C^{2}[ 0,R) $ is radially symmetric with 
$\varphi '(0) =0$. Then, for $k\in \{1,2,\dots ,N\} $ and $u(x) =\varphi (r) $ 
where $r=| x| <R$, we have that $u\in C^{2}(B_{R}) $, and
\begin{gather*}
\lambda (D^{2}u(r) ) 
=\begin{cases}
(\varphi ''(r) ,\frac{\varphi '(r) }{r},\dots ,\frac{\varphi '(r) }{r}
) & \text{for }r\in (0,R) , \\
(\varphi ''(0) ,\varphi ''(0) ,\dots ,\varphi ''(0) ) & \text{for }r=0\,;
\end{cases}
\\
S_k(\lambda (D^{2}u(r) ) ) 
=\begin{cases}
C_{N-1}^{k-1}\varphi ''(r)\big(\frac{\varphi '(r)}{
r}\big) ^{k-1}+C_{N-1}^{k-1}\frac{N-k}{k}\big(\frac{\varphi '(r)
}{r}\big) ^{k}, & r\in (0,R) , \\
C_{N}^{k}(\varphi ''(0) ) ^{k} &\text{for }r=0,
\end{cases}
\end{gather*}
where the prime denotes differentiation with respect to $r=|x| $ and 
$C_{N-1}^{k-1}=(N-1)!/[ (k-1)!(N-k)!] $.
\end{remark}


\subsection*{Proof of the Theorems \ref{thm1} and \ref{thm2}}

We start by showing that  system \eqref{11} has positive radial
solutions. For this purpose, we show that the system of ordinary differential
equations 
\begin{equation}
\begin{gathered}
\begin{aligned}
&\frac{C_{N-1}^{k_1-1}}{r^{N-1}}[\frac{r^{N-k_1}}{k_1}e^{\int_{0}^{r}
\frac{1}{C_{0}}t^{k_1-1}b_1(t) dt}(u_1')^{k_1}]'\\
&=e^{\int_{0}^{r}\frac{1}{C_{0}}t^{k_1-1}b_1(
t) dt}p_1(r) f_1(u_1,u_2) ,\quad r>0,
\\
&\frac{C_{N-1}^{k_2-1}}{r^{N-1}}[\frac{r^{N-k_2}}{k_2}e^{\int_{0}^{r}
\frac{1}{C_{00}}t^{k_2-1}b_2(t) dt}(u_2')^{k_2}]'\\
&=e^{\int_{0}^{r}\frac{1}{C_{00}}t^{k_2-1}b_2(t) dt}p_2(r) f_2(u_1,u_2) , r>0,
\end{aligned}\\
u_1'(r) \geq 0\quad \text{and}\quad
u_2'(r) \geq 0\quad\text{for }r\in [ 0,\infty ) , \\
u_1(0) =a_1\quad \text{and}\quad u_2(0) =a_2,
\end{gathered}  \label{ss}
\end{equation}
has a solution. Therefore, at least one solution of \eqref{ss} can be obtained
using successive approximation by defining the sequences 
$\{u_1^{m}\} _{m\geq 0}$ and $\{ u_2^{m}\} _{m\geq 0}$ on 
$[ 0,\infty ) $ in the following way
\begin{equation}
\begin{gathered}
u_1^{0}=a_1, \quad  u_2^{0}=a_2\quad\text{for }r\geq 0 \\
u_1^{m}(s) =a_1+\int_{0}^{r}\Big[ B_1^{-}(t)
\int_{0}^{t}B_1^{+}(s) f_1(u_1^{m-1}(s)
,u_2^{m-1}(s) ) ds\Big] ^{1/k_1}dt, \\
u_2^{m}(s) =a_2+\int_{0}^{r}\Big[ B_2^{-}(t)
\int_{0}^{t}B_2^{+}(s) f_2(u_1^{m-1}(s)
,u_2^{m-1}(s) ) ds\Big] ^{1/k_2}dt.
\end{gathered}  \label{recs}
\end{equation}
It is easy to see that
$\{ u_1^{m}\} _{m\geq 0}$  and $\{ u_2^{m}\}_{m\geq 0}$
are non-decreasing on $[ 0,\infty ) $. Indeed, we consider
\begin{align*}
u_1^{1}(r) 
&= a_1+\int_{0}^{r}\Big[B_1^{-}(t)
\int_{0}^{t}B_1^{+}(s) f_1(u_1^{0}(s),u_2^{0}(s) ) ds\Big]^{1/k_1}dt \\
&= a_1+\int_{0}^{r}\Big[B_1^{-}(t) \int_{0}^{t}B_1^{+}(s) f_1(a_1,a_2) ds
\Big]^{1/k_1}dt \\
&\leq a_1+\int_{0}^{r}\Big[B_1^{-}(t)
\int_{0}^{t}B_1^{+}(s) f_1(u_1^{1}(s),u_2^{1}(s) ) ds\Big]^{1/k_1}dt=u_1^{2}(r) .
\end{align*}
This implies that $u_1^{1}(r) \leq u_1^{2}(r) $
which further produces $u_1^{2}(r) \leq u_1^{3}(r) $. 
Continuing, an induction argument applied to \eqref{recs} show
that for any $r\geq 0$ we have
\begin{equation*}
u_1^{m}(r) \leq u_1^{m+1}(r) \quad\text{and}\quad
u_2^{m}(r) \leq u_2^{m+1}(r) \quad \text{for any } m\in \mathbb{N}
\end{equation*}
i.e., $\{ u_1^{m}\} _{m\geq 0}$ and $\{ u_2^{m}\}_{m\geq 0}$ are non-decreasing 
on $[ 0,\infty ) $. By their monotonicity, 
we have the inequalities
\begin{gather}
C_{N-1}^{k_1-1}\{\frac{r^{N-k_1}}{k_1}e^{\int_{0}^{r}\frac{1}{C_{0}}
t^{k_1-1}b_1(t) dt}[ (u_1^{m}) '] ^{k_1}\}' 
\leq B_1^{+}(r) f_1(u_1^{m},u_2^{m}) ,  \label{gen1} 
\\
C_{N-1}^{k_2-1}\{\frac{r^{N-k_2}}{k_2}e^{\int_{0}^{r}\frac{1}{C_{00}}
t^{k_2-1}b_2(t) dt}[(u_2^{^{m}}) ']^{k_2}\}' 
\leq B_2^{+}(r) f_2(u_1^{m},u_2^{m}) .  \label{gen2}
\end{gather}
After integration from $0$ to $r$, an easy calculation yields
\begin{equation}
\begin{aligned}
&(u_1^{m}(r) ) '\\
&\leq \Big(B_1^{-}( r) \int_{0}^{r}B_1^{+}(t) f_1(u_1^{m}(
t) ,u_2^{m}(t) ) dt\Big)^{1/k_1}   \\
&\leq \Big(B_1^{-}(r) \int_{0}^{r}B_1^{+}(t) f_1(u_1^{m}(t) +u_2^{m}(t)
,u_1^{m}(t) +u_2^{m}(t) ) dt\Big)^{1/k_1}   \\
&\leq (f_1^{1/k_1}+f_2^{1/k_2})(u_1^{m}(r)
+u_2^{m}(r) ,u_1^{m}(r) +u_2^{m}(r) ) (B_1^{-}(r)
\int_{0}^{r}B_1^{+}(t) dt)^{1/k_1}.
\end{aligned} \label{exin}
\end{equation}
As before, exactly the same type of conclusion holds for
$(u_2^{m}(r) ) '$:
\begin{equation}
\begin{aligned}
(u_2^{m}(r) ) '
&\leq \Big(B_2^{-}(r) \int_{0}^{r}B_2^{+}(z) f_2(u_1^{m}(z) ,u_2^{m}(z) ) dz\Big)^{1/k_2}   \\
&\leq (f_1^{1/k_1}+f_2^{1/k_2})(
u_1^{m}+u_2^{m},u_1^{m}+u_2^{m}) \Big(B_2^{-}(r)
\int_{0}^{r}B_2^{+}(t) dt\Big)^{1/k_2}.
\end{aligned} \label{exin2}
\end{equation}
Summing the  inequalities \eqref{exin} and \eqref{exin2}, we
obtain
\begin{equation}
\frac{(u_1^{m}(r) +u_2^{m}(r) )'}{(f_1^{1/k_1}+f_2^{1/k_2})(u_1^{m}(r)
+u_2^{m}(r) ,u_1^{m}(r) +u_2^{m}(r) ) }
\leq P_1'(r) +P_2'(r) .  \label{mat}
\end{equation}
Integrating from $0$ to $r$, we obtain
\begin{equation*}
\int_{a_1+a_2}^{u_1^{m}(r) +u_2^{m}(r) }
\frac{1}{f_1^{1/k_1}(t,t) +f_2^{1/k_2}(t,t) }
dt\leq P_1(r) +P_2(r) .
\end{equation*}
We now have
\begin{equation}
F_{1,2}(u_1^{m}(r) +u_2^{m}(r) )
\leq P_1(r) +P_2(r) ,  \label{ints}
\end{equation}
which will play a basic role in the proof of our main results. The
inequalities \eqref{ints} can be reformulated as
\begin{equation}
u_1^{m}(r) +u_2^{m}(r) \leq F_{1,2}^{-1}(
P_1(r) +P_2(r) ) .  \label{int}
\end{equation}
This can be easily seen from the fact that $F_{1,2}$ is a bijection with the
inverse function $F_{1,2}^{-1}$ strictly increasing on $[ 0,\infty) $.
So, we have found upper bounds for
$\{ u_1^{m}\} _{m\geq 0}$  and $\{ u_2^{m}\}_{m\geq 0}$
which are dependent of $r$. We are now ready to give a complete proof of the
Theorems \ref{thm1} and \ref{thm2}.

\begin{proof}[Proof of Theorem \ref{thm1} completed]
 When $F_{1,2}(\infty) =\infty $ it follows that the sequences 
$\{ u_1^{m}\}_{m\geq 0}$ and $\{ u_2^{m}\} _{m\geq 0}$ are bounded and
equicontinuous on $[ 0,c_{0}] $ for arbitrary $c_{0}>0$. By the
Arzela-Ascoli theorem, $\{ (u_1^{m},u_2^{m}) \} _{m\geq 0}$ has
a subsequence converging uniformly to $(u_1,u_2) $ on 
$[0,c_{0}] \times [ 0,c_{0}] $. Since $\{u_1^{m}\} _{m\geq 0}$ and 
$\{ u_2^{m}\} _{m\geq 0}$ are
non-decreasing on $[ 0,\infty ) $ we see that 
$\{ (u_1^{m},u_2^{m}) \} _{m\geq 0}$ itself converges uniformly
to $(u_1,u_2) $ on $[ 0,c_{0}] \times [0,c_{0}] $. At the end of this process, 
we conclude by the arbitrariness of $c_{0}>0$, that $(u_1,u_2) $ is positive
entire solution of system \eqref{ss}. The solution constructed in this way
will be radially symmetric. Since the radial solutions of the ordinary
differential equations system \eqref{ss} are solutions \eqref{11} it follows
that the radial solutions of \eqref{11} with $u_1(0) =a_1$,
 $u_2(0) =a_2$ satisfy
\begin{gather}
u_1(r) = a_1+\int_{0}^{r}\Big(B_1^{-}(y)
\int_{0}^{y}B_1^{+}(t) f_1(u_1(t)
,u_2(t) ) dt\Big)^{1/k_1}dy,   \label{eq1} \\
u_2(r) = a_2+\int_{0}^{r}\Big(B_2^{-}(y)
\int_{0}^{y}B_2^{+}(t) f_2(u_1(t),u_2(t) ) dt\Big)^{1/k_2}dy,  \label{eq2}
\end{gather}
for all $r\geq 0$. Next, it is easy to verify that the Cases 1. and 2. occur.
\smallskip

\noindent\textbf{Case 1.} When
$P_1(\infty ) +P_2(\infty ) <\infty$,
it is not difficult to deduce from \eqref{eq1} and \eqref{eq2} that
\begin{equation*}
u_1(r) +u_2(r) \leq F_{1,2}^{-1}(P_1(\infty ) +P_2(\infty ) ) <\infty \quad
\text{for all }r\geq 0,
\end{equation*}
and so $(u_1,u_2) $ is bounded. We next consider:
\smallskip

\noindent\textbf{Case 2.}
When $P_1(\infty ) =P_2(\infty ) =\infty$,
we observe that
\begin{equation}
\begin{aligned}
u_1(r) &= a_1+\int_{0}^{r}\Big(B_1^{-}(t)
\int_{0}^{t}B_1^{+}(s) f_1(u_1(s),u_2(s) ) ds\Big)^{1/k_1}dt   \\
&\geq a_1+f_1^{1/k_1}(a_1,a_2)
\int_{0}^{r}\Big(B_1^{-}(t) \int_{0}^{t}B_1^{+}(s)
ds\Big)^{1/k_1}dt  \label{i1} \\
&= a_1+f_1^{1/k_1}(a_1,a_2) P_1(r) .
\end{aligned}
\end{equation}
The same computations as in \eqref{i1} yields
\begin{equation*}
u_2(r) \geq a_2+f_2^{1/k_2}(a_1,a_2) P_2(r) .
\end{equation*}
and passing to the limit as $r\to \infty $ in \eqref{i1} and in the
above inequality we conclude that
\begin{equation*}
\lim_{r\to \infty }u_1(r) =\lim_{r\to \infty
}u_2(r) =\infty ,
\end{equation*}
which yields the result.
\end{proof}


\begin{proof}[Proof of Theorem \ref{thm2} completed]
In view of the above analysis, the proof can be easily deduced from
\begin{equation*}
F_{1,2}(u_1^{m}(r) +u_2^{m}(r) )
\leq P_1(\infty ) +P_2(\infty ) <F_{1,2}(
\infty ) <\infty ,
\end{equation*}
Indeed, since $F_{1,2}^{-1}$ is strictly increasing on 
$[ 0,\infty) $, we find that
\begin{equation*}
u_1^{m}(r) +u_2^{m}(r) 
\leq F_{1,2}^{-1}(P_1(\infty ) +P_2(\infty ) ) <\infty ,
\end{equation*}
and then the non-decreasing sequences
$\{ u_1^{m}\} _{m\geq 0}$  and $\{ u_2^{m}\}_{m\geq 0}$
are bounded above for all $r\geq 0$ and all $m$. The final step, is to
conclude that
\begin{equation*}
(u_1^{m}(r) ,u_2^{m}(r) )
\to (u_1(r) ,u_2(r) ) \quad \text{as }m\to \infty
\end{equation*}
and the limit functions $u_1$ and $u_2$ are positive entire bounded
radial solutions of system \eqref{11}.
\end{proof}

\begin{remark} \rm
Make the same assumptions as in Theorem \ref{thm1} or Theorem \ref{thm2} on 
$b_1$, $p_1$, $f_1$, $b_2$, $p_2$, $f_2$. If, in addition,
\begin{gather}
p_1(| x| )  \geq \Big(C_{N-1}^{k_1-1}\frac{
N-k_1}{k_1| x| ^{N}}-\frac{b_1(|x| ) }{| x| ^{N-k_1}}
\Big)\int_{0}^{| x| }\frac{s^{N-1}}{C_{0}}p_1(s) ds,\quad x\in \mathbb{R}^{N}, 
 \label{ing1} \\
p_2(| x| )  \geq \Big(C_{N-1}^{k_2-1}\frac{
N-k_2}{k_2| x| ^{N}}-\frac{b_2(|x| ) }{| x| ^{N-k_2}}
\Big)\int_{0}^{| x| }\frac{s^{N-1}}{C_{00}}p_2(s) ds,\quad x\in \mathbb{R}^{N}, 
 \label{ing2}
\end{gather}
then the solution $(u_1,u_2) $ is convex.
\end{remark}

\begin{proof} 
It is clear that
\begin{equation}
\frac{C_{N-1}^{k_1-1}}{r^{N-1}}\Big[\frac{r^{N-k_1}}{k_1}e^{\int_{0}^{r}
\frac{1}{C_{0}}t^{k_1-1}b_1(t) dt}(u_1')^{k_1}\Big]'
=e^{\int_{0}^{r}\frac{1}{C_{0}}t^{k_1-1}b_1(t) dt}p_1(r) f_1(u_1,u_2) ,  \label{jo}
\end{equation}
and integrating from $0$ to $r$ yields
\begin{align*}
&r^{N-k_1}e^{\int_{0}^{r}\frac{1}{C_{0}}t^{k_1-1}b_1(t)
dt}(u_1'(r) )^{k_1} \\
&= \int_{0}^{r}\frac{s^{N-1}}{C_{0}}e^{\int_{0}^{s}\frac{1}{C_{0}}
t^{k_1-1}b_1(t) dt}p_1(s) f_1(
u_1(s) ,u_2(s) ) ds \\
&\leq  f_1(u_1(r) ,u_2(r) )
\int_{0}^{r}\frac{s^{N-1}}{C_{0}}e^{\int_{0}^{s}\frac{1}{C_{0}}
t^{k_1-1}b_1(t) dt}p_1(s) ds,
\end{align*}
which yields
\begin{equation}
\begin{aligned}
(\frac{u_1'(r) }{r})^{k_1}
&\leq \frac{f_1(u_1(r) ,u_2(r) ) }{
r^{N}e^{\int_{0}^{r}\frac{1}{C_{0}}t^{k_1-1}b_1(t) dt}}
\int_{0}^{r}\frac{s^{N-1}}{C_{0}}e^{\int_{0}^{s}\frac{1}{C_{0}}
t^{k_1-1}b_1(t) dt}p_1(s) ds   \\
&\leq \frac{f_1(u_1(r) ,u_2(r) )}{r^{N}}\int_{0}^{r}\frac{s^{N-1}}{C_{0}}p_1(s) ds.
\end{aligned}\label{conv}
\end{equation}
On the other hand  inequality \eqref{jo} can be written in the form
\begin{equation}
\begin{aligned}
&C_{N-1}^{k_1-1}u_1''(r) (\frac{
u_1'}{r})^{k_1-1}+C_{N-1}^{k_1-1}\frac{N-k_1}{k_1}(
\frac{u_1'}{r})^{k_1}+b_1(r) (u_1')^{k_1} \\
&=p_1(r) f_1(u_1,u_2) .
\end{aligned} \label{fc}
\end{equation}
Using  inequality \eqref{conv} in \eqref{fc} we obtain
\begin{align*}
p_1(r) f_1(u_1,u_2)
&\leq C_{N-1}^{k_1-1}u_1''(\frac{u_1'}{r})^{k_1-1}
 +C_{N-1}^{k_1-1}\frac{N-k_1}{r^{N}k_1}f_1(
u_1,u_2) \int_{0}^{r}\frac{s^{N-1}}{C_{0}}p_1(s) ds
\\
&\quad +\frac{b_1(r) f_1(u_1,u_2) }{r^{N-k_1}}
\int_{0}^{r}\frac{s^{N-1}}{C_{0}}p_1(s) ds,
\end{align*}
from which we have
\begin{align*}
&f_1(u_1,u_2) [p_1(r) -(C_{N-1}^{k_1-1}
\frac{N-k_1}{k_1r^{N}}-\frac{b_1(r) }{r^{N-k_1}}
)\int_{0}^{r}\frac{s^{N-1}}{C_{0}}p_1(s) ds] \\
&\leq C_{N-1}^{k_1-1}u_1''(\frac{u_1'}{r})^{k_1-1},
\end{align*}
which completes the proof of $u_1''(r) \geq 0 $.
A similar argument produces $u_2''(r)\geq 0$.
 We also remark that, in the simple case $b_1=b_2=0$,
 $s^{N-1}p_1(s) $ and $s^{N-1}p_2(s) $ are
increasing then \eqref{ing1} and \eqref{ing2} hold.
\end{proof}

\subsection*{Acknowledgements}
The author would like to thank to the editors and
reviewers for valuable comments and suggestions which contributed to improve
this article.

\begin{thebibliography}{99}

\bibitem{CA} C. O. Alves. A. R. .F de Holanda;
\emph{Existence of blow-up solutions for a class of elliptic systems}, 
Differential Integral Equations,
\textbf{26} (1/2) (2013), 105-118.

\bibitem{BA} C. Bandle, E. Giarrusso;
\emph{Boundary blow-up for semilinear elliptic equations with nonlinear 
gradient terms}, Advances in Differential Equations, \textbf{1} (1996), 133--150.

\bibitem{BAOII} J. Bao, X. Ji, H. Li;
\emph{Existence and nonexistence theorem for entire subsolutions of k-Yamabe 
type equations}, J. Differential Equations \textbf{253} (2012), 2140--2160.

\bibitem{CL} P. Cl\'{e}ment, R. Man\'{a}sevich, E. Mitidieri;
 \emph{ Positive solutions for a quasilinear system via blow up}, Communications in
Partial Differential Equations, \textbf{18 }(12) (1993), 2071--2106.

\bibitem{CD} D.-P. Covei;
\emph{Existence and non-existence of solutions for an elliptic system},
 Applied Mathematics Letters, \textbf{37} (2014), 118-123.

\bibitem{FI} D. G. De Figueiredo, Y. Jianfu;
\emph{Decay, symmetry and existence of solutions of semilinear elliptic systems}, 
Nonlinear Analysis: Theory, Methods \& Applications, \textbf{33} (1998), 211--234.

\bibitem{GA} V. Galaktionov, J.-L. V\'{a}zquez;
\emph{The problem of blow-up in nonlinear parabolic equations}, 
Discrete and Continuous Dynamical Systems - Series A, \textbf{8} (2002), 399--433.

\bibitem{GI} E. Giarrusso;
\emph{On blow up solutions of a quasilinear elliptic equation}, 
Mathematische Nachrichten, \textbf{213} (2000), 89--104.

\bibitem{GR} X. Jiang, X. Lv;
\emph{Existence of entire positive solutions for semilinear elliptic systems 
with gradient term}, Archiv der Mathematik, \textbf{99} (2) (2012), 169-178.

\bibitem{K} J. B. Keller;
\emph{On solution of $\Delta u=f(u)$},
Communications on Pure and Applied Mathematics, \textbf{10} (1957), 503-510.

\bibitem{J} H. Jian;
 \emph{Hessian equations with infinite Dirichlet boundary value}, 
Indiana University Mathematics Journal, \textbf{55} (2006), 1045--1062.

\bibitem{BAO} X. Ji, J. Bao;
\emph{Necessary and sufficient conditions on solvability for Hessian inequalities}, 
Proceedings of the American Mathematical Society, \textbf{138} (1), 175--188.

\bibitem{LY} Q. Li, Z. Yang;
\emph{Large solutions to non-monotone quasilinear elliptic systems}, 
Applied Mathematics E-Notes, \textbf{15} (2015), 1-13.

\bibitem{Z} H. Li, P. Zhang, Z. Zhang;
\emph{A remark on the existence of entire positive solutions for a class of 
semilinear elliptic systems}, Journal of Mathematical Analysis and Applications, 
\textbf{365} (2010), 338--341.

\bibitem{O} R. Osserman;
\emph{On the inequality $\Delta u\geq f(u)$},
Pacific Journal of Mathematics, \textbf{7} (1957), 1641-1647.

\bibitem{PW} J. Peterson, A. W. Wood;
\emph{Large solutions to non-monotone semilinear elliptic systems}, 
Journal of Mathematical Analysis and Applications, \textbf{384} (2011), 284-292.

\bibitem{PRI} C. Pripoae;
\emph{Non-Holonomic Economical Systems},
Geometry Balkan Press, \textbf{10} (2004), 142-149.

\bibitem{Q} P. Quittner;
 \emph{Blow-up for semilinear parabolic equations with a gradient term}, 
Mathematical Methods in the Applied Sciences, \textbf{14} (1991), 413--417.

\bibitem{R} H. Rademacher;
\emph{Finige besondere probleme partieller Differentialgleichungen}, in:
 Die Differential und Integralgleichungen der
Mechanick und Physik I, 2nd ed., Rosenberg, New York, Pages 838-845, 1943.

\bibitem{SA} P. Salani;
\emph{Boundary blow-up problems for Hessian equations}, 
Manuscripta Mathematica, \textbf{96} (1998), 281--294.

\bibitem{V} J. A. Viaclovsky;
\emph{Conformal geometry, contact geometry, and the calculus of variations}, 
Duke Mathematical Journal, \textbf{101} (2000), 283--316.

\bibitem{VI} J. A. Viaclovsky;
\emph{Estimates and existence results for
some fully nonlinear elliptic equations on Riemannian manifolds},
Communications in Analysis and Geometry, \textbf{10} (2002), 815--846.

\bibitem{ZZ} Z. Zhang, S. Zhou;
\emph{Existence of entire positive k-convex radial solutions to Hessian 
equations and systems with weights},
Applied Mathematics Letters, \textbf{50} (2015), 48--55.

\bibitem{ZEJ} Z. Zhang;
\emph{Existence of entire positive solutions for a class of semilinear elliptic 
systems},  Electronic Journal of Differential
Equations, \textbf{2010} (16) (2010), Pages 1-5.

\end{thebibliography}

\end{document}
