\documentstyle[twoside]{article}
\input amssym.def     % used for R in Real numbers
\pagestyle{myheadings}
\markboth{\hfil Nonexistence of Positive Singular Solutions 
\hfil EJDE--1996/08}%
{EJDE--1996/08\hfil Cecilia S. Yarur \hfil}
\begin{document}
\title{\vspace{-1in}\parbox{\linewidth}{\footnotesize\noindent
{\sc  Electronic Journal of Differential Equations},
Vol.\ {\bf 1996}(1996), No.\ 08, pp. 1--22. \newline
ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu
\newline ftp (login: ftp) 147.26.103.110 or 129.120.3.113}
 \vspace{\bigskipamount} \\
Nonexistence  of Positive Singular Solutions for
a  Class of Semilinear Elliptic Systems
\thanks{ {\em 1991 Mathematics Subject Classification:}
35J60, 31A35.\newline\indent
{\em Key words and phrases:} Elliptic systems, Removable singularity, Biharmonic equation.
\newline\indent
\copyright 1996 Southwest Texas State University  and University of
North Texas.\newline\indent
Submitted June 6, 1996. Published September 6, 1996.} }
\date{}
\author{Cecilia S. Yarur}
\maketitle
\begin{abstract}
We study nonexistence and removability results for nonnegative  
sub-solutions to
$$\left. \begin{array}{rcl}
\Delta u &=& a(x) v^p \\
\Delta v &=& b(x) u^q \end{array}\right\}
\mbox{ in } \Omega \subset \Bbb R^N,\quad N\ge 3\,, 
$$
where  $p\geq 1$, $q\geq 1$, $pq>1$, and  $a$ and $b$ are nonnegative 
functions. As a consequence of this work, we obtain new results for 
biharmonic equations.
\end{abstract}

\renewcommand{\theequation}{\arabic{section}.\arabic{equation}}
\newtheorem{teorema}{Theorem}[section]
\newtheorem{lema}{Lemma}[section]
\newtheorem{coro}{Corollary}[section]

\section{ Introduction}
The aim of this paper  is to study  nonexistence and removability results 
for nonnegative solutions of the inequality system
\begin{equation}\left. \begin{array}{rcl}
\Delta u &\ge& a(x) v^p \\
\Delta v &\ge& b(x) u^q \end{array}\right\}
\mbox{ in } \Omega \subset \Bbb R^N,\quad N\ge 3 
\end{equation}
where   $p\ge 1$ , $q\ge 1$ and $pq>1$.
We assume that  the functions $a$ and $b$ are nonnegative functions 
defined in $L_{\rm loc}^{\infty}\left(\Omega\right)$.

We will give a unified treatment for the cases  $\Omega = \Bbb{R}^N$, 
$ \Omega = B_1(0) \backslash \{0\} \;$  and $\Omega = \Bbb{R}^N \backslash \{0\}$ 
in (1.1). For this purpose we will base our arguments essentially on 
a priori bounds results for (1.1) in the one-dimensional case in exterior domains 
(Theorem 2.1, Theorem 2.2 and Corollary 2.1 below). 

One reason for tackling this type of problem is the study of 
nonnegative solutions for the semilinear biharmonic equation
\begin{equation}
\Delta^2 u =  u^q 
 \;\;\; \mbox {in \quad } \Bbb{R}^N,\  N\ge 3\,. 
\end{equation}
As a consequence of our results for system (1.1) we will prove   
that all the nonnegative nontrivial  solutions  of (1.2) are 
super-harmonic functions in $\Bbb{R}^N$ (Corollary~3.1). Then, for 
instance, nonexistence results of positive super-harmonic functions for 
(1.2) proved by Mitidieri in \cite{mi1,mi2} are now nonexistence 
results of positive solutions for the biharmonic equation. 

Moreover,  the system 
\begin{equation}
\left. \begin{array}{rcl}
-\Delta u &=&  |v|^{p-1}v \\
-\Delta v &=&  u^{q-1}u \end{array}\right\}\quad
\mbox{ in }\Omega \subset \Bbb R^N,\quad N\ge 3 
\end{equation}
with $u$ positive and $v$ negative can be treated as a particular 
case of (1.1). For  the system (1.3) we refer to \cite{sp,  vd} and the references therein.

In the case that  $\Omega = \Bbb{R}^N$, we will assume that
$a$ and $b$ in $(1.1)$ satisfy the following condition at infinity:
\begin{equation}
\begin{array}{lll}
 a_p(|x|):= \left(\frac{1}{|S_{N-1}|}\int_{S_{N-1}}
a(|x|\sigma)^{-1/(p-1)} \; d\sigma \right)^{1-p}  & \ge & c |x|^{-\alpha} \cr
 b_q(|x|):= \left(\frac{1}{|S_{N-1}|}\int_{S_{N-1}}
b(|x|\sigma)^{-1/(q-1)} \; d\sigma \right)^{1-q}  & \ge & c |x|^{-\beta}, \cr
\end{array}  
\end{equation}
for some  positive constant $c$. Let us define 
\begin{equation}
\gamma_1(\alpha, \beta) = {{\alpha-2 + (\beta -2)p }\over{pq-1}}
\mbox{ and }\gamma_2(\alpha, \beta) = 
{{\beta-2 + (\alpha -2)q }\over{pq-1}}\,.
\end{equation}
 Our main result for the system (1.1)  in $\Bbb{R}^N$ reads as follows  

\paragraph{Theorem 3.4} {\it Let $(u,v) \in \left(C\left({\Bbb{R}}^N\right)
\right)^2$ be a positive solution of} (1.1).
{\it Let} $p \ge 1, q \ge 1$ {\it and }$pq>1$. {\it Assume  $a$ and $b$ 
are nonnegative functions  defined in $\Bbb{R}^N$ satisfying} (1.4)   
{\it for $ |x|$ near infinity with 
$\alpha, \beta$   such that }
$$
\min \left\{\gamma_1(\alpha,\; \beta),\gamma_2(\alpha,\; \beta)  \right\} 
\le 0.
$$
{\it Then $u\equiv 0$  and $v\equiv  0$.}
\bigskip

Ni \cite{ni} has proven that, for $\alpha < 2$, the  equation
\begin{equation}
\Delta u = a(x)u^q \qquad  \mbox{in} \quad \Bbb{R}^N
\end{equation}
does not have any positive solution. This result was improved by 
 F.H. Lin \cite{flin} for  $\alpha \le 2$. On the other hand
 for   $\alpha > 2$, Ni \cite{ni}, and Naito \cite{na}, among others, 
have proven existence results. In this case, there is no sign  restriction
 for the function $a$, but now $|a(x)|\le c|x|^{-\alpha}.$
Thus  $\alpha= 2$ is a {\it critical exponent} for the equation (1.6) 
in $\Bbb{R}^N$. We point out that for  equation (1.6) we have  
 $ \alpha= \beta$, $p= q$,  and
$\gamma_1= \gamma_2= \frac{\alpha - 2}{p-1}$. Thus, the critical 
exponent
$\alpha= 2$ is represented  now by $\min\{\gamma_1, \gamma_2\}= 0$. 
Therefore, Theorem 3.4 generalizes the early works  
\cite{ni} and \cite{flin} to the nonlinear system (1.1).
In exterior domains  the behavior near infinity  of any  solution $u$ of
\begin{equation}
\Delta u = |u|^{q-1}u,
\end{equation}
has been given by  V\'eron  \cite{ve}. 

In the case the 
$\Omega = B_1(0) \backslash \{0\}$,  we are interested in removability 
results for system (1.1), that is, when   all nonnegative solutions of 
(1.1) are bounded at zero and satisfy (1.1) in the sense of 
distributions in ${\cal D}'(B_1(0))$. The main  result that we  will 
prove in this direction is  the following.

\paragraph {Theorem 4.3} {\it Let $(u,v) \in \left(C\left(B_1(0) \backslash \{0\}\right)\right)^2$ 
be a positive solution
of} (1.1)   {\it  in  $B_1(0) \backslash \{0\}$}.  {\it Let 
 $p \ge 1, q \ge 1$  and $pq>1$.  Assume  $a$ and $b$ are nonnegative
functions  defined in $B_1(0) \backslash \{0\}$ satisfying} (1.4) {\it 
 for $ |x|$ near $0$, 
with
$\alpha, \beta$  such that, either }
 \begin{enumerate}
\item[ {\rm(i)}] $\min \left\{\gamma_1(\alpha,\; \beta), \; 
\gamma_2(\alpha,\; \beta)\right \} \ge 2-N$, {\it or}
 \item[ {\rm (ii)}] $\max \left\{\gamma_1(\alpha,\; \beta), \; 
\gamma_2(\alpha,\; \beta)\right \} \ge 2-N $, 
$p\ge (2-\alpha)/(N-2)$, and \newline $q \ge (2-\beta)/(N-2)$.
\end{enumerate}
{\it Then $u$ and $v$ are bounded near zero, and  
$(u, v)$ satisfies (1.1)  in ${\cal{D}}'( B_1(0))$.
\bigskip

Loewner and Nirenberg  \cite{ln} proved removability results 
for (1.7)  
with  $p=(N+2)/(N-2)$.  Later,  Br\`esis and V\'eron  \cite{bv}  
 improved the Loewner-Nirenberg result for $p \ge N/(N-2)$.
 If $1<p<N/(N-2)$, there are solutions of (1.7) with isolated 
singularities.
Therefore,  for   equation (1.7), the critical exponent for 
removability results in  a ball 
is $p=N/(N-2)$, which  is exactly the  condition (i) ( or (ii)) in  
Theorem 4.3.

Finally, in the case that $\Omega=\Bbb{R}^N \backslash \{0\}$, we prove 
nonexistence  of nonnegative solutions  (singular or not) for 
the system (1.1).
 We remark that for the equation
\begin{equation}
\Delta u -V(|x|)u= a(x)u^{p},
\end{equation}
nonexistence of nonnegative  sub-solutions was proven in 
\cite{bly} under decay conditions on
$a(x )$ for $x$ near zero and infinity. For  existence results for (1.8) 
see also \cite{sy} .

The rest of the  paper is organized as follows:
In Section 2, we give some preliminary results for the one dimensional 
 case in (1.1) in exterior domains. Section 3 is devoted 
to the cases where $\Omega$ in (1.1) is 
 either the whole space or   an exterior domain and in 
Section 4 we study removability results for (1.1). Finally, in  Section 5
we prove nonexistence results in $\Bbb{R}^N \backslash \{0\}$.

\section{ Preliminary results}  \setcounter{equation}{0}
In this section we prove some results that are needed later in the 
proof of our main theorems. The first two lemmas
are proven in  \cite{bly} (see also \cite{ni} for the second one).
We also need the spherical average of a function $f$, which is defined  by
$$
\bar{f}(r)=\frac{1}{|S_{N-1}|}\int_{S_{N-1}}f(r\sigma)d\sigma,
$$
where $d\sigma$ denotes the invariant measure on the sphere
$$
S_{N-1}=\left\{x\in\Bbb{R}^N\;:\;\sum_{i=1}^{N}x_{i}^{2}=1\right\}.
$$
Here, $|S_{N-1}|$ denotes the volume of the unit sphere. We denote by 
$\Bbb{R}_{0}^N$ the set $\Bbb{R}^{N}\backslash\{0\}$.
We will say that $(u, v) \in\left( C\left(\Omega\right)\right)^2$ is a 
nonnegative  solution of (1.1)  if $u$ and $v$ are nonnegative in 
$\Omega$ and $(u, v)$   satisfies (1.1) in $ {\cal{D}}' (\Omega)$.

The following lemma is  a nonexistence result for 
positive sub-harmonic functions with prescribed behavior at zero and at 
infinity (see \cite{bly})

\begin{lema} Let $u \ge 0 \in L_{\rm loc}^{1}(\Bbb{R}_0^N)$ such that $\Delta u \ge 0$ 
and assume
\begin{equation}
\lim_{r \rightarrow 0}r^{N-2}\bar u(r) = 0
\end{equation}
and
\begin{equation}
\lim_{r \rightarrow \infty}\bar u(r) = 0.
\end{equation}
Then $u \equiv 0$ in $\Bbb{R}^N$.
\end{lema}

The next lemma is used to reduce the study of a partial differential 
problem to the  study of an ordinary differential one (see \cite{bly} 
and \cite{ni})

\begin{lema} Let $ f(x,t) = a(x)t^p$,$ \; a(x) \ge 0$,  $p \ge 1$
and let $v $ be a nonnegative function. Then 

\begin{equation}
\overline{av^p}(|x|) \ge a_p(|x|){\bar v}^p(|x|)\,,
\end{equation}
where 
$$
a_p(r) = \left(\frac{1}{|S_{N-1}|}\int_{S_{N-1}}
a(r\sigma)^{-1/(p-1)} \; d\sigma \right)^{1-p} \qquad \mbox {for}\quad p>1
$$
and $a_1(r) =\min_{ \sigma \in  S_{N-1}} a(r\sigma)$ for $p=1$.
If $\int_{S_{N-1}} a(r\sigma)^{-1/(p-1)} \, d{\sigma}  = \infty$,
we put $a_1(r) =0$.  
\end{lema}

Having reduced the partial differential problem to an 
ordinary differential one, 
we need some previous results for solutions of system (1.1) in 
 one dimension. To begin with, we give some power  solutions for the 
system
\begin{eqnarray}
(r^{N-1}u')' &=& ar^{N-1-\alpha}v^p \nonumber\\
(r^{N-1}v')' &=& br^{N-1-\beta}u^q\,, 
\end{eqnarray}  
with $a$ and $b$ positive constants, 
which will play an important role in determining the regions of nonexistence 
as well as bounds for the solutions of (1.1). This is not surprising, 
 since for the equation 
  
\begin{equation}
(r^{N-1}u')' = ar^{N-1-\alpha}u^q
\end{equation}
those solutions have an outstanding role, too.
 If we try to get solutions to (2.4) of power type, that is 
\begin{equation}
\begin{array}{rcl}
u(r)&=&l_1 r^{\gamma_1(\alpha, \beta)} \\
v(r)&=&l_2 r^{\gamma_2(\alpha, \beta)} \,.
\end{array}
\end{equation}
We then find that  $l_1, l_2, \gamma_1$ and $\gamma_2$ must satisfy
\begin{equation}
 \begin{array}{rcl}
l_1 \gamma_1 (\gamma_1 +N-2)&=& a {l_2}^p  \\
l_2 \gamma_2 (\gamma_2 +N-2)&=& b {l_1}^q  
\end{array} 
\end{equation}
and
\begin{equation}
 \gamma_1(\alpha, \beta) = {{\alpha-2 + (\beta -2)p }\over{pq-1}}\,,
\quad
 \gamma_2(\alpha, \beta) = {{\beta-2 + (\alpha -2)q }\over{pq-1}}\,.
\end{equation}
We write at our convenience $\gamma_1(\alpha, \beta)$ and 
$ \gamma_2(\alpha, \beta)$, but $\gamma_1,\gamma_2$
 certainly depend also on $p$ and  $q$. 

The existence of positive constants $l_1, l_2$ which satisfy (2.7) 
is equivalent to 
$$\gamma_i (\gamma_i+N-2)>0, \quad  \mbox{ for} \quad i=1,2\,.$$
We observe that for $N \ge 3$ and $ \;\min\{\gamma_1, \gamma_2\} > 0$,
we get the existence   of power  solutions
 for the system (2.4) in the whole space. This fact is very relevant
 in view of Theorem 3.4. Moreover, for some values of $\alpha, \beta, 
 p$ and $q$ we have existence of a solution  of (2.4) satisfying (2.6)
in $\Bbb{R}^N \backslash \{0\}$, with $u$ bounded near zero and $v$ 
going to infinity and vice versa.  

Now, we state the main results of this section, that belong to the case  
 $N= 1$ for the system (1.1).  Theorem 2.1 and
Theorem 2.2, or their  equivalents in higher dimensions
 (Theorem 3.1 and Theorem 3.2), 
will be the key to demonstrate nonexistence results for the coming sections.
 The proof will be shown at the end of this section because some 
preliminary lemmas are required. 

\begin{teorema} Let $(w_1, w_2)$ be a nonnegative solution of
\begin{equation} \left.
\begin{array}{rcl}
\ddot{w}_1(s) & \geq & c_1 s^{-\delta_1} w_{2}^{p}\\
\ddot{w}_2(s) & \geq & c_2 s^{-\delta_2} w_{1}^{q}
\end{array}\right\}    \qquad\mbox{for all}\qquad s\ge s_0,
\end{equation}
for some  $s_0$ positive.  Assume that  $p,q >0$ and $pq>1$. Moreover,
 we assume that either 
\begin{enumerate}
\item[{\rm (i)}] $\gamma_1(\delta_1, \delta_2) \le 1$, or 
\item[{\rm(ii)}]$\gamma_2(\delta_1, \delta_2) \le 1$ and ${\delta}_2 \le q+2$.
\end{enumerate}

\noindent Then $w_1$ is bounded. 
\end{teorema} 

Similarly, we have the following 

\begin{teorema} Let $(w_1, w_2)$ be a nonnegative solution of (2.9) with 
  $p,q >0$ and $pq>1$. Moreover, we assume that either 

\begin{enumerate}
\item[ {\rm (i)}] $\gamma_2(\delta_1, \delta_2) \le 1$, or  
\item[{\rm (ii)}]$\gamma_1(\delta_1, \delta_2) \le 1$ and ${\delta}_1 \le p+2$.
\end{enumerate}
\noindent Then $w_2$ is bounded. 
\end{teorema}

\begin{coro} Let $(w_1, w_2)$ be a nonnegative solution of (2.9) 
 with  $p,q >0$ and $pq>1$. Moreover, we assume that
\begin{equation}
\min{ \left\{\gamma_1(\delta_1, \delta_2), \;\gamma_2(\delta_1, \delta_2) \right\}} \le 1.
\end{equation}
\noindent Then $w_1$ or $w_2$ is bounded.
\end{coro}

\begin{coro}Let $(w_1, w_2)$ be a nonnegative solution of (2.9) 
 with    $p,q >0$ and $pq>1$. Moreover, assume that either

\begin{enumerate}
\item[{\rm (i)}] $\max \left\{ \gamma_1(\delta_1, \delta_2),
\gamma_2(\delta_1,\; \delta_2) \right\}\le 1$, or
\item[{\rm (ii)}] ${\delta}_1\le p+2$, ${\delta}_2\le q+2 $ and 
$ \min \left\{ \gamma_1(\delta_1, \delta_2),\gamma_2(\delta_1,\; \delta_2) \right\}\le 1.$
\end{enumerate}
Then $w_1$ and $w_2$ are bounded.
\end{coro}

The next lemma is a generalization of Lemma 2.4 in \cite{bly}
for a systems.

\begin{lema}  Let $p$ and $q$ be two positive real numbers
 such that $pq > 1$, and let  $(w_1, w_2)$ be a nonnegative solution of
\begin{equation}
\begin{array}{lll}
\ddot{w}_1(s) & \geq & X_1(s)w_{2}^{p}\\
\ddot{w}_2(s) & \geq & X_2(s)w_{1}^{q},\\
\end{array}  
\end{equation}
for all $ s \ge s_0$, for some $s_0 > 0$. Here $X_1(s) \ge 0$,  
$ X_2(s) \ge 0$ are continuous and non-increasing functions on $s \ge s_0.$ 
Moreover, we assume the following hypotheses:
\begin{itemize} 
\item[{\rm (H1)}] $\int^{\infty}X_1(s)\;s^{p}\;ds = \infty$ and 
$\int^{\infty}X_2(s)\;s^{q}\;ds = \infty$,
\item[{\rm (H2)}] There exist three positive constants  $\alpha_1 > 1$, 
$\alpha_2 > 1$, and $c$ such that
$$
\frac{\alpha_1}{p+1} + \frac{\alpha_2}{q+1} = 1,$$
and, for all $s $ large enough 
$$
\max\{s^{{-\alpha_1 + 1}}, s^{{-\alpha_2 + 1}}\} \le 
c \int_{s}^{\infty}\, X_1(s)^{\alpha_1/(2(p+1))} 
X_2(s)^{\alpha_2/(2(q+1))}\,ds\,.
$$
\end{itemize}
Then $w_1$ and $w_2$ are bounded.
\end{lema}

\paragraph{\bf Remarks}  In the above lemma  we have    that
$\alpha_1= \alpha_2 = (p+1)/2$ for the equation (2.5).      
Lemma 2.3 can be generalized  for more general functions than $t^p$ 
and $t^q$. 

\paragraph{Proof of  Lemma 2.3.} First, we will show that 
it is enough to consider the case in which $w_1$ and 
$w_2$  are both unbounded near infinity. This fact will
be fundamentally a consequence of the hypothesis  (H1).

Since $X_2$ is nonnegative, the function $w_2$ is convex and we have  
 the following two possibilities.  Either:
\begin{itemize}
\item[{\rm (a)}]\quad $\dot{w}_2(s)\leq 0$, for all $s$, or
\item[{\rm (b)}]\quad there is an $s_1$, such that $\dot{w}_2(s)>0$ for 
all $s\geq s_1.$
\end{itemize}
If (a) holds then $w_2$ is bounded. If we assume that $w_1$ is not bounded,
then $\dot{w}_1(s)\ge 0$ for all
large $s$, then  since $ w_1$ is convex we get,
 $w_1(s)\ge cs$ for some constant $c$ positive.
By integrating (2.11), it follows that 
$$
\begin{array}{lll}
\dot{w}_2(s) & \geq & \dot{w}_2(s_1)+\int_{s_1}^{s}X_2\;w_{1}^{q}(t)\;dt\\
\\
             & \geq & \dot{w}_2(s_1)+c\int_{s_1}^{s}X_2\;t^{q}\;dt.\\  
\end{array}
$$
Hence from (H1), $\dot{w}_2(s)$ goes to infinity as $s\to\infty$, 
which contradicts (a).
Thus we conclude that  $w_1$ is bounded if  $w_2$ is bounded.

Now, if (b) holds, arguing as in case (a) we also have $\dot{w}_1(s)>0$ for large $s$, 
and $\dot{w_i}(s)$ goes to infinity as $s\to\infty$, for $i= 1,2$. 
 Therefore, we can assume that  $w_1$ and $w_2$ are both unbounded.

Now, multiplying the first inequality in (2.11) by $\dot{w}_2$ and the
second one  by $\dot{w}_1$ and  then adding both expressions,  we get
\begin{equation}
\frac{d}{ds}(\dot{w}_1\dot{w}_2)\geq X_1\frac{d}{ds}
\left(\frac{w_{2}^{p+1}}{p+1}\right)+X_2\frac{d}{ds}\left(\frac{w_{1}^{q+1}}{q+1}\right)
\end{equation}
for all $s\ge \tilde{s}$, for some $\tilde{s}$. Integrating (2.12 ) 
from $\tilde{s}$ to $s$ we have 
\begin{equation}  
\dot{w}_1\dot{w}_2(s)\geq\int_{\tilde{s}}^{s}X_1\frac{d}{ds}
\left(\frac{w_{2}^{p+1}}{p+1}\right)+\int_{\tilde{s}}^{s}X_2\frac{d}{ds}
\left(\frac{w_{1}^{q+1}}{q+1}\right).
\end{equation}

Moreover, since  $X_1$ and $X_2$ are non-increasing functions for large $s$,
 from (2.13) we get   
\begin{equation}
\dot{w}_1\dot{w}_2(s)\geq X_1(s)\left( \frac{w_{2}^{p+1}}{p+1}(s)  - 
\frac{w_{2}^{p+1}}{p+1}(\tilde{s}) \right)
 + X_2(s) \left(\frac{w_{1}^{q+1}}{q+1}(s)- 
\frac{w_{1}^{q+1}}{q+1}(\tilde{s})\right).
\end{equation}

If $s$ is large enough,   $s\geq s_2$ for some $s_2$, we can take
$w_i(s)\geq\frac{1}{2}\;w_i(\tilde{s})$, for $i= 1,2$, and we obtain 
\begin{equation}
\dot{w}_1(s)\dot{w}_2(s)\geq c\left(X_1(s)\;w_2^{p+1}(s)+
X_2(s)\;w_1^{q+1}(s)\right),
\end{equation}
for all $s\geq s_2$. Here $c$ is a positive constant.

Now, we use the following relation between the geometric and arithmetic 
means
\begin{equation}
 {a_1}^{p_1} {a_2}^{p_2} \le \left( {{p_1 a_1 + p_2 a_2}
\over{p_1 + p_2}}\right) ^{p_1+p_2} 
\end{equation}
where $a_1, a_2, p_1,$ and $ p_2 $ are positive  numbers.
 We can choose $p_1$ and $p_2$ as follows
$$
\frac{p_1}{p_1+p_2} = \frac{ \alpha_1}{p+1}, \quad 
\frac{p_2}{p_1+p_2} = \frac{ \alpha_2}{q+1}.
$$
Then if we apply (2.16) into (2.15) with $a_1$ and $a_2$ defined by 
$${{p_1 a_1}\over{p_1+p_2}}=X_1 {w_2}^{p+1}, 
\quad {{p_2 a_2}\over{p_1+p_2}}=X_2 {w_1}^{q+1}$$
we get
$$
{{\dot{w}_1 \dot{w}_2 (s)}\over{{w}_1^{\alpha_2} {w_2}^{\alpha_1}}} 
\ge c X_1^{\alpha_1/(p + 1)} X_2^{\alpha_2/(q + 1)}.
$$
Hence,
$$
{{\left(\dot{w}_1 \dot{w}_2 \right)^{1/2}}
\over{{w_1}^{\alpha_2/2} {w_2}^{\alpha_1/2}}} 
\ge c X_1^{\alpha_1/(2(p+1))} X_2^{\alpha_2/(2(q+1))}
$$
which in turn implies
$$
{{\dot{w}_1 \over{w_1}^{{\alpha_2}}}} +
{{\dot{w}_2 \over{w_2}^{{\alpha_1}}}}
\ge c X_1^{\alpha_1/(2(p+1))} X_2^{\alpha_2/(2(q+1))}
$$
for all $s \ge s_2$. Then integrating from $s \ge s_2$ to $\infty$ we 
get
$$
\int^{\infty}_{w_1(s)} {{dt}\over {t^{{\alpha_2}}}} + 
\int^{\infty}_{w_2(s)} {{dt}\over {t^{{\alpha_1}}}} \ge
c \int^{\infty}_{s} X_1^{\alpha_1/(2(p+1))} X_2^{\alpha_2/(2(q+1))}\;dt,
$$
which because of (H2), and since
 $\lim_{s \rightarrow \infty}w_i(s)/s = + \infty$,
 for $i=1,2$, gives us a contradiction.

The next result is a particular case of the above lemma, and is the 
key for proving the main results of this section.   

\begin{lema} Let $(w_1, w_2)$ be a nonnegative solution of (2.9).
Assume that  $p,q >0$ and $pq>1$. Moreover, assume that  $\delta_1 \le p+1$ 
 and $\delta_2 \le q+1$. Then  $w_1$ and $w_2$ are bounded near infinity.
\end{lema}
\smallskip
{\bf Proof.} \quad Let us call ${\delta_i}^{+}= \max \{ \delta_i, 0 \}$
 for $ i=1,2$. We can take on the above lemma,
 $X_i= c_i s^{-{\delta_i}^{+}},$  for $i=1,2.$ 
Then $X_1$ and $X_2$ are non-increasing functions, and $(w_1, w_2)$ is a 
nonnegative solution of 
$$
\begin{array}{lll}
\ddot{w}_1& \ge& X_1 w_2^p \cr
\ddot{w}_2& \ge& X_2  w_2^q,\cr
\end{array}
$$
for all $s$ large.  We have to prove 
 the validity of the conditions (H1) and (H2) given on the above result.
\begin{enumerate}
\item[(H1):]  $\int^{\infty} X_1 s^p=\infty$, is equivalent with
${\delta_1}^{+} \le p+1$, which is satisfied since  ${\delta_1}\le p+1$. 
In the same way 
$\int^{\infty} X_2 s^q=\infty$,  since  ${\delta_2}\le q+1$.
\item[(H2):]  We have to find $\alpha_1$ and $\alpha_2$ satisfying 
 condition (H2) on Lemma 2.3. Let us denote $x = \alpha_1/(p+1)$ and 
$y =  \alpha_2/(q+1)$. 
The problem of finding  $ \alpha_1$ and 
$\alpha_2 $  is reduced to find $x, y$ which verify the  following 
conditions
\begin{eqnarray*}
&x+y=1, \quad  x>{1\over{p+1}},\quad y>{1\over {q+1}},& \\ 
&\left( 2(q+1) - {\delta_2}^{+}\right) y \ge {\delta_1}^{+} x,
\mbox{ and } 
\left( 2(p+1) - {\delta_1}^{+}\right) x \ge {\delta_2}^{+} y\,.&
\end{eqnarray*}

Let 
$$ a={{{\delta_2}^{+}}\over {2(p+1) -{\delta_1}^{+} + 
{\delta_{2}}^{+}}}
\mbox{ and } 
b={{2(q+1)-{\delta_2}^{+}}\over {{\delta_1}^{+}+
2(q+1) - {\delta_2}^{+}}}\,.$$
Then $a$ and $b$ are well defined and $(a, 1-a)$ is the intersection of the 
lines $x+y=1$, 
 $\left( 2(p+1) - {\delta_1}^{+}\right) x = {\delta_2}^{+} y$
and $(b, 1-b)$ is the intersection of $x+y = 1$ with
 $\left( 2(q+1) - {\delta_2}^{+}\right) y = {\delta_1}^{+} x.$

Now, since $pq> 1$ and $\delta_1^+ \le p+1,\; \delta_2^+ \le q+1$, we always have 
$$
a< \frac{q}{q+1} \quad \mbox{and} \quad    \frac{1}{p+1}< b.
$$
Also $a \le b$, so that

$$ A \equiv \max \left\{ {1\over{p+1}},a \right\} \le  
\min \left\{ {q\over{q+1}},b \right\} \equiv B.$$

If $A \ne B$, we can choose any $x$ such 
that $A < x < B$. On the contrary, if $A = B$, it can be 
proved that 
 $A = a = b$. In this case, we choose $x = a$. 
\end{enumerate} 


 The above  systems  can have only  one component bounded 
but not the other. This is enough for some of our purposes, as we will see on 
section 3. 
The following two lemmas are concerned with the boundedness of  at least 
one of the 
components of the pair $(w_1,w_2)$.

\smallskip
 
\begin{lema} Let $(w_1,w_2)$ be a  positive solutions of (2.9) for some 
 $p,q >0$ and $pq>1$. Let us call $\bar{\delta}_1\equiv
\delta_1-p-1, \quad \bar{\delta}_2\equiv \delta_2-q-1$.
 Assume that 
$$
\bar{\delta}_1 \le 0 \qquad \mbox{and} \qquad \gamma_2(\delta_1, \delta_2)
\le 1.
$$
Then $w_2$ is bounded. 
\end{lema}

\paragraph {Proof.} \quad The proof is divided into three cases, 
depending on the values of $\delta_1$ and $\delta_2$. 

\paragraph{Case 1:} $\overline{\delta}_1 \le 0$ and $\overline{\delta}_2\le 0$.
 We are in the previous lemma. 

\paragraph{Case 2:} Assume next that \quad $\overline{\delta}_1 < 0$ and
$\gamma_2({\delta}_1, {\delta}_2) < 1.$ The condition $\gamma_2<1$ 
is equivalent to $ \bar{\delta}_1q +\bar {\delta}_2 < 0$.
We proceed by contradiction. If $w_2$ is not bounded, then there exists an $s_0$ such that 
$\dot{w}_2
(s_0)>0$. Now, since $w_2$ is convex we get 
$w_2(s) \ge cs$ for all large $s$  and for some nonnegative constant $c$. 
Going back to (2.9) we get
\begin{equation}
\ddot{w}_1 \ge cs^{-\overline{\delta}_1-1},
\end{equation}
for all $s$ large enough.
Integrating twice from $s_0$ to $s$ in the above inequality and using the fact
 that $ \bar{\delta}_1 < 0$,
 we obtain
\begin{equation} 
w_1(s) \ge cs ^{-\overline{\delta}_1+1},
\end{equation}
for all $s$  large.
 Applying the estimate (2.18) into  (2.9), we have 
the following for $w_2$:
$$w_2(s) \ge cs^{-\overline{\delta}_2 -\overline{\delta}_1 q+1}$$
for all large $s$.
Iterating the above  process, as in \cite{gmmy}, we get for $n \in \Bbb{N}$
$$w_1(s) \ge cs^{p_n}$$
$$w_2(s) \ge cs^{q_n}$$
for $s$ large, where
$$
\begin{array}{rcl}
p_n &=&-{\delta}_1+2 +pq_n \\
q_{n+1}&=&-{\delta}_2+2 +qp_n\\
 q_1&=&1. \end{array}
$$
(The constant $c$ represents any positive value). Due to the condition 
 $\overline{\delta}_1q + \overline{\delta}_2<0,$
  we deduce that  the sequences $\left\{ p_n \right\} $ and 
$\left\{ q_n \right\}$  are strictly increasing.
  Let us call
$$P=\lim_{n\rightarrow \infty} p_n \quad \mbox{and} \quad 
Q= \lim_{n\rightarrow \infty} q_n.$$
Then either  $P=Q=\infty$   or
\begin{equation}
P=-{\delta}_1+2 +pQ \mbox{ and }
Q=-{\delta}_2+2 +qP\,.
\end{equation}
Thus, multiplying the first equation on (2.19) by  $q$ and adding the 
 second one, we get
$$0=-\overline{\delta}_1q-\overline{\delta}_2+(Q-1)(pq-1),$$
which  is a contradiction to   $\; Q > q_1= 1$ and $-\overline{\delta}_1q-
\overline{\delta}_2>0$.

Now, if $P=Q=\infty$, then for all $p'$ and $q'$ with,
$ p'<p$ and $ q'<q$, we have 
\begin{eqnarray*}
&\ddot{w}_1 \ge cs^{-{\delta_1}} {w_2}^p \ge w_2^{p'}&\\  
&\ddot{w}_2 \ge cs^{-{\delta_2}} {w_1}^q \ge w_1^{q'}\,.
\end{eqnarray*}
Moreover, choosing  $p'$ and $q'$ such that $p'q'>1$,  from Lemma 2.3, we 
deduce  that 
$w_1$ and  $w_2$ are bounded which is a contradiction.

\paragraph{Case 3:} $\gamma_2({\delta}_1,{\delta}_2)=1$ and 
$\overline{\delta}_1<0$.
As in the previous case, we proceed by contradiction. If $w_2$ is not 
bounded, we claim that  for all $k>0$ 
$$
\lim_{s\to \infty}{\frac{w_1(s)}{ s^k}} = \infty \mbox{ and }
\lim_{s\to \infty}{\frac{w_2(s)}{ s^k}} = \infty\,.  
$$
If the claim is true, then arguing as we did at the end of Case 2, we will 
get a contradiction.
Next we will prove the claim.
Since we are assuming that $w_2$ is not bounded, one can prove the 
following  estimate for $w_2$ near infinity:
$$
w_2(s) \ge cs\log s\,,
$$
so that
\begin{equation}
\lim_{s\to \infty}\frac{w_2(s)}{ s} = \infty\,.
\end{equation}
Also, $ w_1$ and $w_2$ are increasing functions for large $s$. 
Integrating the first inequality on (2.9) from $s$ to $2s$, we get
\begin{equation}
\dot{w}_1(2s) \ge \dot{w}_1(2s)-\dot{w}_1(s)\ge
c \int^{2s}_s t^{-\delta_1} {w_2}^p(t)\;dt.
\end{equation}
Hence,
\begin{equation}
\dot{w}_1(2s) \ge  c \int^{2s}_s t^{-\delta_1} 
 w_2^p(t)\; dt 
\ge  c w_2^p (s) s^{-{\delta}_1+1}
\end{equation}
Integrating (2.22) from $s$ to $2s$, and arguing as above, we get
\begin{equation}
w_1(4s) \ge c {w_2}^p(s) s^{-\delta_1 + 2}
\end{equation}
In the same way, but now starting  with the second inequality on (2.9), 
we get 
\begin{equation}
 w_2(4s) \ge c {w}_1^q(s) s^{-\delta_2 + 2}
\end{equation}
If we use (2.23) in (2.24) we obtain
$$w_2 (16s) \ge c {w_2}^{pq}(s) s^{-\overline{\delta}_1q -
\overline{\delta}_2+1-pq}.$$
From the hypothesis $\overline{\delta}_1q+\overline{\delta}_2=0$,
we then  have 
\begin{equation}
w_2(16s) \ge c {w}_2^{pq}(s) s^{1-pq}.
\end{equation}
We rewrite (2.25) in the form
\begin{equation}
 {{w_2(16s)}\over{16s}} \ge c \left({{{w_2}(s)}
\over s}\right)^{pq}.
\end{equation}
For $n\in \Bbb{N}$, choose  $s=2^{4n}$ in (2.26), and 
$\displaystyle x_n= c^{1/(pq-1)}w_2 (2^{4n})/2^{4n}$. Then
\begin{equation}
x_{n+1} \ge  {x_n}^{pq},
\end{equation} 
for all $n$ large,  $n \ge n_0$, for some $n_0$.
A repeated iteration on (2.27) leads to the estimate
$$ 
x_{n+1} \ge  {x_{n_0}}^{(pq)^{n+1-n_0}}, 
$$
for all $n \ge n_0$. From  (2.20),  $x_n \rightarrow \infty$ as 
$n\rightarrow\infty$, then  we can take $n_0 \in \Bbb{N}$   such that
$$ x_{n_0} > 1\,.$$
Therefore, for all $\beta >0$ we obtain 
$$\lim_{n\rightarrow \infty} {{x_{n+1}}\over
{\left( 2^{4(n+1)}\right)^{\beta}}} =\infty.$$
Going back  to the definition of $x_n$, we deduce 
$$\lim_{n\rightarrow \infty} {{w_2\left( 2^{4n} \right)}\over
{\left( 2^{4n} \right)^{\beta+1}}} =\infty.$$

Next, we prove that 
$\displaystyle \lim_{s\rightarrow \infty} {{w_2(s)}\over{s^{\beta+1}}}
=\infty$. Let $s$ be sufficiently large and $n\in \Bbb{N}$ be such that
$s\in \left[ 2^{4n}, {2}^{4(n+1)}\right]$. Since  $w_2(s)$ is 
nondecreasing, then
$$ {{w_2(s)}\over{s^{\beta+1}}} \ge {{w_2\left( 2^{4n} \right)}\over
{ 2^{4(n+1)(\beta+1)}}},$$
which implies 
$\displaystyle \lim_{s \rightarrow\infty}{{w_2(s)}\over{s^{\beta+1}}}
=\infty$, for all $\beta>0$ and the claim follows from (2.23).

In analogous form, we obtain
\begin{lema} Let $(w_1,w_2)$ be a  positive solution of (2.9) with 
 $p,q >0$ and $pq>1$. Let us call $\bar{\delta}_1\equiv
\delta_1-p-1, \quad \bar{\delta}_2\equiv \delta_2-q-1$.
 Assume that 
$$
\bar{\delta}_2 \le 0 \qquad \mbox{and} \qquad \gamma_1(\delta_1, \delta_2)
 \le 1.
$$
Then $w_1$ is bounded. 
\end{lema}

In the following lemmas, we prove  that for certain values of 
$\delta_1, \delta_2, p $ and $q$ in (2.9), if one component of the pair
$(w_1, w_2)$  is bounded, then the other is bounded, too. This  allows
extending the regions of boundedness of $w_1 $ and $w_2$ obtained in 
previous lemmas.

\begin{lema} Let $(w_1, w_2)$ be a  positive solution of (2.9) 
with $p,q>0$ and $pq>1$. Assume that $w_2$ is bounded and 
$$\min\left\{ \gamma_1(\delta_1, \delta_2), \;\bar{\delta}_2 \right\} \le 1.$$
Then $w_1$ is bounded.
\end{lema}
With respect to the boundedness of  $w_2$, assuming boundedness of $w_1$, we have 

\begin{lema} Let $(w_1, w_2)$ be a  positive solution of (2.9)  
with $p,q>0$ and $pq>1$. Assume that $w_1$ is bounded and 
$$\min\left\{ \gamma_2(\delta_1, \delta_2), \;\bar{\delta}_1 \right\} \le 1.$$
Then $w_2$ is bounded.
\end{lema}
 
\paragraph{Proof of Lemma 2.7.} We distinguish two cases, according to
whether $\gamma_1 \le 1$ or 
$\bar{\delta}_2 \le 1$. We assume first 

\paragraph{Case 1:}  $\gamma_1 \le 1$. This case is equivalent to
 $\bar{\delta}_1 +p\bar{\delta}_2 \le 0$.
Now,  since $w_2$ is bounded at infinity it must be a non-increasing function for 
all $s$ large. Suppose by contradiction that $w_1$ is 
not bounded near infinity. Then $w_1$
is increasing for $s$ large enough. Integrating the first inequality 
on (2.9) from $s/2$ to $s$ it follows that
\begin{equation}
\dot{w}_1(s) \ge c \left( \int^s_{s/2} t^{-\delta_1}\right)
{w_2}^p (s)
\ge c s^{-\delta_1+1} {w_2}^p(s)\,. 
\end{equation}
Integrating once again from $s/2$ to $s$ in (2.28), we obtain
\begin{equation}
 w_1(s) \ge c s^{-\delta_1+2}{w_2}^p(s).
\end{equation}
Similarly, but now integrating from  $s$ to $2s$ in the second inequality of
 (2.9), we get 

\begin{equation}
w_2(s) \ge c s^{-\delta_2+2}{w_1}^q(s).
\end{equation}
Therefore, by using  (2.29) and (2.30), 
in the first inequality of (2.9) we have  the following for $w_1$
\begin{equation}
\begin{array}{rcl}
\ddot{w}_1 &\ge& c s^{-\delta_1+p(-\delta_2+2)} {w_1}^{pq} \\
&\equiv& c s^{-\gamma} {w_1}^{pq}, 
\end{array}
\end{equation}
where $\gamma =\delta_1 - p(-\delta_2 +2)$. By the assumption
$\overline{\delta}_1 +p\overline{\delta}_2\le 0$, it follows that
$\gamma \le pq+1$. Thus, by Lemma 2.4 (see also \cite{bly}) $w_1$ must be bounded.

\paragraph{Case 2:} $ \bar{\delta}_2  \le 1$. Assume that $w_1$ is 
not bounded, then  $w_1 \ge cs$\quad for $s$ large. As before, 
from (2.9) it follows that
$$w_2(s) \ge c s^{-\overline{\delta}_2+1}\,,$$
which in turn implies that 
 $w_2 \rightarrow \infty $ as $ s \rightarrow \infty$
if $\overline{\delta}_2 <1$.
Now, if $\overline{\delta}_2=1$\quad we get  the same conclusion by
integrating in
$$\ddot{w}_2(s) \ge c s^{-\delta_2} {w_1}^q \ge c s^{-1}.$$


\paragraph{Remark} Theorem 2.1 is a consequence of 
Lemma2 2.5, 2.6, and 2.7.  
Similarly, Theorem 2.2 is a consequence of the Lemmas 2.5, 
2.6 and 2.8.

\section{Nonexistence in $\Bbb{R}^N$} \setcounter{equation}{0}
In  this section we consider $\Omega$ in (1.1) to be either an exterior 
domain, for instance  $\Omega = \{x : |x| \ge 1 \}$, or 
$\Omega= \Bbb{R}^N$. For exterior 
domains, 
we will give  bounds near infinity for one or both of the 
components of the pair $(u, v)$, where $(u, v)$ is a  nonnegative 
solution of (1.1) (Theorem 3.1,  
Theorem 3.2 and Theorem 3.3).
  In the whole space we will prove
a nonexistence result, Theorem 3.4, for nonnegative  
nontrivial solutions of (1.1).  
We remark that Theorem 3.4 is optimal for the system (2.4). 

 Throughout this section we will assume that $a$ and $b$ are nonnegative
functions 
in $L_{\rm loc}^{\infty}\left(\Omega\right)$. Moreover, there exist three 
constants $\alpha$, 
$\beta$ and $c$, with $c$ positive, such that  
\begin{equation}\left.
\begin{array}{rcl}
 a_p(|x|) & \ge & c |x|^{-\alpha} \cr
 b_q(|x|) & \ge & c |x|^{-\beta} \cr
\end{array}\right\} \qquad \mbox{at infinity,} 
\end{equation}
 where $a_p$ and $b_q$ are defined in Lemma 2.2. 

\begin{teorema} Let $(u,v) \in \left( C\left(|x| \ge 1 \right) \right)^2$  be 
a positive solution of
\begin{equation}\left.
\begin{array}{rcl}
\Delta u &\ge& a(x) v^p \cr
\Delta v &\ge& b(x) u^q \cr
\end{array}\right\} \qquad \mbox{in} \qquad |x| \ge 1\,,
\end{equation}
where $p\ge 1, q\ge 1$ and  $ \; pq>1$. Assume  $a$ and $b$ are nonnegative
functions  defined in $|x| \ge 1$  and satisfying (3.1) with 
$\alpha, \beta$  such that either
\begin{enumerate}
\item[ {\rm (i)}]\quad  $\gamma_1(\alpha,\; \beta) \le 0$, or
\item[ {\rm (ii)}]\quad  $\gamma_2(\alpha,\; \beta) \le 0$ and $\beta \le N.$
\end{enumerate}
Then $|x|^{N-2}u$ is bounded.
\end{teorema}

For the equation (1.6) the  conditions ${\rm (i)}$ and ${\rm (ii)}$ in Theorem 3.1  
are equivalent with $\alpha\le 2$. 

Before  proving  Theorem 3.1 let us enunciate the boundedness for $v$.

\begin{teorema} Let $(u,v) \in \left( C\left(|x| \ge 1 \right)\right)^2$  be a positive solution
of (3.2). Let  
  $p\ge 1, q\ge 1$ and $ \; pq>1$. Assume  $a$ and $b$ are nonnegative
functions  defined in $|x| \ge 1$  and satisfying (3.1) with
$\alpha, \beta$  such that either

\begin{enumerate}
\item[ {\rm (i)}]\quad  $\gamma_2(\alpha,\; \beta) \le 0$, or 
\item[ {\rm (ii)}] \quad  $\gamma_1(\alpha,\; \beta) \le 0$ and $\alpha \le N.$
\end{enumerate}
Then $|x|^{N-2}v$ is bounded.
\end{teorema}

\paragraph {Proof of Theorem 3.1.} From (3.1), (3.2)  and Lemma 2.2, 
we have 
\begin{equation}
\begin{array}{lll}
\bar{u}'' + \frac{N-1}{r}\bar{u}' & \ge & c r ^{-\alpha}\;\bar{v}^p\\

\bar{v}'' + \frac{N-1}{r}\bar{v}' & \ge & c r^{-\beta}\;\bar{u}^q,\\
\end{array}
\end{equation}
for all $r$ large enough. Let  $s=r^{N-2}$ and let
$$
\begin{array}{rcl}
w_1(s) & = & s\bar{u}(r)\\
w_2(s) & = & s\bar{v}(r)\,.
\end{array}
$$
Then $w_1$ and $w_2$ satisfy
\begin{equation}
\begin{array}{rcl}
\ddot{w}_1(s) & \geq & c s^{-\delta_1} w_{2}^{p}\\
\ddot{w}_2(s) & \geq & c s^{-\delta_2}w_{1}^{q}
\end{array}
\end{equation}
where 
$$
\delta_1 =  \frac{\alpha-2}{N-2}+p+1  
\quad\mbox{and}\quad \delta_2 =  \frac{\beta-2}{N-2}+q+1.
$$
It follows from the hypothesis on $\alpha, \beta, p, q$ and  Theorem 2.1 that 
$w_1$ is bounded. Thus, from the definition of $w_1$, we get that $r^{N-2}\bar u$ is bounded. To prove that  
 $|x|^{N-2}u$ is also  bounded we use the following mean value inequality
for sub-harmonic functions (see \cite{gt})
$$
u(x) \le \frac1{|B_{|x|/2}(x)|}\int_{B_{|x|/2}(x)}u(y) \;dy,
$$
then 
\begin{equation}
u(x) \le c |x|^{-N}\int_{|x|/2}^{3|x|/2} r^{N-1}\bar{u}(r) \;dr.
\end{equation}
Since $r^{N-2}\bar u$ is bounded  for $r$ large enough and $u$ satisfies (3.5), then 
  the conclusion of the  theorem follows. 

Next we apply the previous results to the biharmonic.

 \begin{coro} Let  $q > 1$ and $u \in  C^2\left({\Bbb{R}}^N\right)$ be 
a positive solution of 
\begin{equation}
\Delta^2 u = b(x) u^q 
 \qquad \mbox{in} \qquad  \Bbb{R}^N
\end{equation}
 Assume that $b$ is a  nonnegative
function  defined in $\Bbb{R}^N$  and satisfying 
$$
b_q(x) \ge c |x|^{-\beta}, \qquad \mbox{for all} \quad |x| \quad \mbox{large},
$$
with $\beta \le 2(q+1)$. Then $u$ is a super-harmonic function in $\Bbb{R}^N$.
\end{coro}

\paragraph{Proof.} Let us define $v:= \Delta u$. Then, the pair 
$(u, v)$ is a solution for
 \begin{equation}\left.
\begin{array}{rcl}
\Delta u &=&  v \\
\Delta v &=& b(x) u^q 
\end{array}\right\} \qquad \mbox{in} \qquad  \Bbb{R}^N\,.
\end{equation}
Since $v$ is a sub-harmonic function in $\Bbb{R}^N$ we get the following
two possibilities for $\bar v$, either 

\begin{enumerate}
\item[{\rm (1)}] \quad There is a positive $r_0$ so that $\bar v(r)\ge 0$, for all 
$r$ larger than $r_0$. Moreover, 
$\lim_{r \rightarrow\infty}{r^{N-2} \bar v(r)}= \infty$, or 
\item[{\rm (2)}] \quad  $\bar v(r) \le 0$, for all $r > 0$.
\end{enumerate}

Theorem 3.2  and the hypothesis on $\beta$ imply that case 1 is 
impossible and then $\bar v \le 0$. Repeating the above argument for 
the functions $v_{y}(x):= v(x+y)$ with $y \in \Bbb{R}^N$, we obtain that 
$\overline{v_{y}} \le 0$ for all $y$. Then the conclusion follows.
As a consequence of the two previous theorems we obtain the following, 
which gives us the boundedness of  $u$  and  $v$ at the same time.

\begin{coro} Let $(u,v) \in \left( C\left( |x| \ge 1 \right)\right)^2$  
be a positive solution of  (3.2). Let  
  $p\ge 1, q\ge 1$ and $ \; pq>1$. Assume  $a$ and $b$ are nonnegative
functions  defined in $|x| \ge 1$ satisfying (3.1) with
$\alpha, \beta$  such that either

\begin{enumerate}
\item[{\rm (i)}] \quad $\max \left\{ \gamma_1(\alpha, 
\beta),\gamma_2(\alpha,\; \beta) \right\}\le 0$, or
\item[{\rm (ii)}] \quad $\alpha\le N, \; \beta \le N$  \quad and \quad 
$ \min \left\{ \gamma_1(\alpha,\; \beta), \gamma_2(\alpha,\; \beta)\right\} \le 0$.
\end{enumerate}
Then $|x|^{N-2}u$ and $|x|^{N-2}v$ are bounded.
\end{coro}

Our main result of this section, in a way, extends those of \cite{ni} and 
\cite{flin}. 
 
\begin{teorema} Let $(u,v) \in  \left(C\left({\Bbb{R}}^N\right)\right)^2$ be 
a positive solution of 

\begin{equation} \left.
\begin{array}{rcl}
\Delta u &\ge& a(x) v^p \\
\Delta v &\ge& b(x) u^q 
\end{array}\right\} \qquad \mbox{in} \qquad  \Bbb{R}^N\,,
\end{equation}
Let $p \ge 1, q \ge 1$ and $pq>1$. Assume  $a$ and $b$ are nonnegative
functions  defined in $\Bbb{R}^N$  and satisfying (3.1)  with
$\alpha, \beta$  such that 
\begin{equation}
\min \left\{\gamma_1(\alpha,\; \beta),\gamma_2(\alpha,\; \beta)  \right\} \le 0
\end{equation}
Then $u\equiv 0$ and $v\equiv  0$.
\end{teorema}

\paragraph{Proof.} The proof follows from Lemma 2.1, 
Theorem 3.1, and Theorem 3.2.

\paragraph{Remark.} For the equation  (1.6), condition  (3.9) 
 in the above theorem is 
the well known condition $\alpha \le 2$ (see \cite{ni} and  \cite{flin}).
 If  (3.9) in the above theorem is not satisfied, then $\gamma_1$ and 
$\gamma_2$ are both positive. Therefore, we can get a positive radial 
solution  $(u, v)$
for the system {\rm (2.4)} in $\Bbb{R}^N$, 
with $u(r) = l_1 r^{\gamma_1}$ and $ v(r) = l_2r^{\gamma_2}.$  

\section{Removable singularities} \setcounter{equation}{0}

Br\`esis and V\'eron (\cite{bv}) have proven removable singularities for
nonlinear elliptic equations in a ball. In the sequel we give the same 
type of result but now for a system. To obtain the behavior of solutions
to (1.1) at zero, we use the Kelvin transform together with the results in 
section~3. 
Let $ B_1(0)$ be the open unit ball centered at zero of 
$\Bbb{R}^N$, with $N\ge 3$. 
Throughout  this section the functions $a$ and $b$ are nonnegative functions in 
$L_{\mbox{loc}}^{\infty}\left( B_1(0) \backslash \{0\}\right)$ such that 
\begin{equation}\left.
\begin{array}{rcl}
a_p(|x|)&\ge& c |x|^{-\alpha} \\
b_q(|x|)&\ge& c |x|^{-\beta} 
\end{array}\right\} \qquad \mbox{for all} \; x \; \mbox{small,}
\end{equation}
for some positive constant $c$, and $a_p$ and $b_q$ defined in 
Lemma 2.2.

\begin{teorema} Let $(u,v) \in \left(C\left(B_1(0) \backslash \{0\}\right)\right)^2$ 
be a positive solution of
\begin{equation}\left.
\begin{array}{rcl}
\Delta u &\ge& a(x) v^p \\
\Delta v &\ge& b(x)u^q 
\end{array}\right\} \qquad \mbox{in}  \quad  B_1(0) \backslash \{0\}
\end{equation}
where   $p\ge 1, q\ge1$, and  $pq>1 $. Assume that $a$ and $b$ are 
nonnegative functions satisfying (4.1) with  
$\alpha, \beta $ such that either
\begin{enumerate}
\item[ {\rm (i)}] \quad $\gamma_1(\alpha,\; \beta)\ge 2-N $, or 
\item[ {\rm (ii)}] \quad $\gamma_2(\alpha,\; \beta)\ge 2-N$ and 
$q \ge (2-\beta)/(N-2)$. 
\end{enumerate}
Then $u$ is bounded near zero.
\end{teorema} 


\paragraph{Proof.} This result is a  consequence of those of section 3;
we transform our problem near zero to  a problem  near infinity.
Let $u_1$ and $v_1$ be the Kelvin transform of $u$ and $v$, that is
$$
\left. \begin{array}{rcl}
u_1(x) & = & |x|^{2-N}\;u(x/{|x|^2})\\
v_1(x) & = & |x|^{2-N}\; v(x/{|x|^2})                        
\end{array}\right\} \qquad \mbox{ for} \quad |x|\ge 1\,,
$$
then, $(u_1, v_1)$ satisfies (\cite{gt}) 
\begin{equation}\left.
\begin{array}{rcl}
\Delta{u_1} & \ge & a_1(x)v_1^p \\
\Delta{v_1} & \ge & b_1(x)u_1^q 
\end{array}\right\} \qquad \mbox{ for} \quad  |x|\ge 1\,,
\end{equation}
where $a_1$ and $b_1$ satisfy
$$
\begin{array}{lcccl}
a_1(x) &=&|x|^{(N-2)p-(N+2)} a(x/{|x|^2}) &\ge& c|x|^{-\alpha_1} \cr
b_1(x)&=& |x|^{(N-2)q-(N+2)}b(x/{|x|^2})&\ge& c|x|^{-\beta_1}, 
\end{array}
$$
and  $\alpha_1$, $\beta_1$ are defined by
$$
\begin{array}{rcl}
\alpha_1 &=& N+2-(N-2)p- \alpha \cr
\beta_1 &=& N+2-(N-2)q- \beta. \cr
\end{array}
$$
Then we obtain 
$$
\gamma_i(\alpha_1, \beta_1)= -\gamma_i(\alpha, \beta) -(N-2), \quad \mbox{for i=1,2}.
$$
 From here, we easily get that  $\alpha_1, \beta_1, p$ and $q$ satisfy the 
hypotheses of Theorem 3.1, thus $|x|^{N-2}u_1$ is bounded at infinity. 
Therefore $u$ is bounded near zero.


\paragraph{Remark.} If in the previous theorem,  $p= q$, 
$\alpha = 0 = \beta$ and $u= v$, then we obtain   
 Theorem 1 of \cite{bv}. In this case conditions 
{\rm (i)} and {\rm (ii)} on   Theorem  4.1  are equivalent to  
$ p \ge N/(N-2)$. Analogously, we get for $v$ the following theorem:

\begin{teorema} Let $(u,v) \in \left(C\left(B_1(0) \backslash \{0\}\right)\right)^2$ 
 be a positive solution of (4.2),
where   $p\ge 1, q\ge1$, and  $pq>1 $. Assume that $a$ and $b$ are nonnegative 
functions satisfying (4.1) with  
$\alpha, \; \beta $   such that either
\begin{enumerate}
\item[\rm{ (i)}] \quad $\gamma_2(\alpha,\; \beta)\ge 2-N$, or 
\item[ \rm{(ii)}] \quad $\gamma_1(\alpha,\; \beta)\ge 2-N$ and 
$p \ge (2-\alpha)/(N-2)$. 
\end{enumerate}
Then $v$ is bounded near zero.
\end{teorema} 

The intersection of the region of  $(\alpha, \beta)$ where 
$u$ is bounded  with the region where $v$ is bounded gives us the main 
result of the section.

\begin{teorema} Let $(u,v) \in \left(C\left(B_1(0) \backslash \{0\}\right)\right)^2$ 
be a positive solution of {\rm (4.2)},
where   $p\ge 1, q\ge1$, and  $pq>1 $. Assume that $a$ and $b$ are nonnegative 
functions satisfying {\rm (4.1)} with  
$\alpha, \; \beta $   such that either
\begin{enumerate}
\item[ {\rm (i)}] \quad $\min \left\{\gamma_1(\alpha,\; \beta), \; 
\gamma_2(\alpha,\; \beta)\right \} \ge 2-N$, or 
\item[ {\rm (ii)}] \quad $\max \left\{\gamma_1(\alpha,\; \beta), \; \gamma_2(\alpha,\; \beta)\right \} \ge 2-N$ 
  and 
$p \ge (2-\alpha)/(N-2)$, $q \ge (2-\beta)/(N-2)$. 
\end{enumerate} 
Then $u$ and $v$ are bounded near zero, and
$(u, v)$ satisfies (4.2) in ${\cal{D}}'\left( B_1(0)\right)$.
\end{teorema} 

As a consequence of the above result we can state the following 
for the biharmonic case:

\begin{coro} Let $u \in C^2\left(B_1(0) \backslash \{0\}\right)$ 
be a positive sub-harmonic solution of 
\begin{equation}
\Delta^2 u = |x|^{-\beta}u^q,
\end{equation}
where $q > 1$. Assume that either 
\begin{enumerate}
\item[ {\rm (i)}] \quad $\beta \ge 4$, or
\item[ {\rm (ii)}] \quad $ N > 4$,  $\beta < 4$ and 
$q \ge (N+2-\beta)/(N-2)$.
\end{enumerate}
Then $ u$ is bounded near zero.
Moreover, $u$ satisfies (4.4) in ${\cal{D}}'\left( B_1(0)\right)$.
\end{coro} 

Soranzo \cite{sor} has proven removability  results for nonnegative  
super-harmonic solutions of (4.4). We remark that for a radially symmetric 
nonnegative solution $u$ of (4.4), we get that $u$  is either sub-harmonic 
or super-harmonic near zero.

\section{ Nonexistence of singular solutions in 
$\Bbb{R}^N \backslash \{0\}.$}  \setcounter{equation}{0}
This section is devoted to nonexistence results of nonnegative solutions 
 (singular or not) for (1.1) in $ \Bbb{R}^N \backslash \{0\}$. 
These results  
can be obtained as a consequence of those of the previous sections. 
We give them without proof.

 Throughout  this section the functions $a$ and $b$ are nonnegative 
functions in $L_{\mbox{loc}}^{\infty}\left( \Bbb{R}^N \backslash \{0\}\right)$. 
In some of the next results  we need also the following properties for 
$a$ and $b$
 

\begin{equation}\left.
\begin{array}{rcl}
a_p(|x|)&\ge& c |x|^{-\alpha_{0}} \\
b_q(|x|)&\ge& c |x|^{-\beta_{0}} 
\end{array}\right\} \qquad \mbox{for all} \; x \; \mbox{small}
\end{equation}
and
\begin{equation}\left.
\begin{array}{lll}
a_p(|x|)&\ge& c |x|^{-\alpha_{\infty}} \\
b_q(|x|)&\ge& c |x|^{-\beta_{\infty}} 
\end{array}\right\} \qquad \mbox{for all}  \; x \; \mbox{large enough}
\end{equation}
where $a_p$ and $b_q$ are defined in Lemma 2.2, and $c$ is some positive 
constant.

\begin{teorema} Let $(u,v) \in \left(C\left(\Bbb{R}^N \backslash \{0\}\right)\right)^2$ 
be a positive solution
of
\begin{equation}\left.
\begin{array}{rcl}
\Delta u &\ge& a(x)  v^p \\
\Delta v &\ge& b(x) u^q 
\end{array}\right\}\qquad in \quad  \Bbb{R}^N \backslash \{0\} 
\end{equation}
where  $p\ge 1, q\ge1$, with  $pq>1 $. Moreover, we assume that 
$a$ and $b$ satisfy  (5.1),
 with $\alpha_{0}, \beta_{0}$ satisfying either 
 \begin{enumerate}
\item[{\rm (i)}]$\gamma_1(\alpha_{0},\; \beta_{0})\ge 2-N $, or 
\item[{\rm (ii)}]$\gamma_2(\alpha_{0},\; \beta_{0})\ge 2-N $ and 
$q \ge (2-\beta_0)/(N-2)$ .
\end{enumerate}
Then  the system (5.3) does not possess any positive solution $(u, v)$ 
with  $u$ going to $0$ at infinity.
\end{teorema}
Likewise, we get the following

\begin{teorema} Let $(u,v) \in \left(C\left(\Bbb{R}^N \backslash \{0\}\right)\right)^2$ 
be a positive solution
of (5.3). Let   $p\ge 1, q\ge1$, and  $pq>1 $. Moreover, we assume that 
$a$ and $b$ satisfy (5.1),
 with $\alpha_{0}, \beta_{0}$ satisfying either 
 \begin{enumerate}
\item[{\rm (i)}]$\gamma_2(\alpha_{0},\; \beta_{0})\ge 2-N$, or
\item[{\rm (ii)}]$\gamma_1(\alpha_{0},\; \beta_{0})\ge 2-N $ and 
$p \ge (2-\alpha_0)/(N-2)$. 
\end{enumerate}
Then  the system (5.3) does not posses any positive solution $(u, v)$ 
with $v$ going to $0$ at infinity.
\end{teorema}

In \cite{bly} Benguria, Lorca, and Yarur prove, among others, the nonexistence of nonnegative 
 singular solutions for the equation (1.6), with decay conditions on 
$a(x)$ for $x$ near zero and infinity. Our next two results extend those 
of   \cite{bly} to the system (5.3).

\begin{teorema} Let $(u,v) \in \left(C\left({\Bbb{R}}_{0}^N 
\right)\right)^2$  be 
a positive solution of (5.3). Let  
 $p\ge 1,q \ge 1$ and $pq>1.$ Moreover,  we assume that 
$ a(x),  \; b(x)$  satisfies (5.1) and (5.2). Suppose that  
$ \alpha_{\infty}$ and  $\beta_{\infty}$
 are such that either 
\begin{enumerate}
\item[ {\rm (i)}]\quad  $\gamma_1(\alpha_{\infty},\; \beta_{\infty}) 
\le 0$,  or
\item[ {\rm (ii)}]\quad  $\gamma_2(\alpha_{\infty},\; \beta_{\infty}) 
\le 0$ 
and $\beta_{\infty} \le N.$
\end{enumerate}
For  $\alpha_{0}$ and  $ \beta_{0}$ we assume that either
\begin{enumerate}
\item[ {\rm (i)}$_0$]\quad  $\gamma_1(\alpha_{0},\; \beta_{0}) \ge 2-N$,
or 
\item[ {\rm (ii)}$_0$]\quad  $\gamma_2(\alpha_{0},\; \beta_{0}) \ge 2-N$, and 
$ q \ge (2-\beta_0)/(N-2)$.
\end{enumerate}
Then $u\equiv 0$ and $v\equiv  0$.
\end{teorema}

\begin{teorema} Let $(u,v) \in \left(C\left({\Bbb{R}}_{0}^N \right)\right)^2$  be 
a positive solution of (5.3). Let  
 $p\ge 1$, $q \ge 1$ and $pq>1$. Moreover,  we assume that 
$ a(x),  \; b(x)$  satisfy (5.1) and (5.2). Suppose that  
$ \alpha_{\infty}$ and  $\beta_{\infty}$ are such that either 
\begin{enumerate}
\item[ {\rm (i)}]\quad  $\gamma_2(\alpha_{\infty},\; \beta_{\infty}) 
\le 0$, or 
\item[ {\rm (ii)}]\quad  $\gamma_1(\alpha_{\infty},\; \beta_{\infty}) 
\le 0$ and $\alpha_{\infty} \le N$.
\end{enumerate}
For  $\alpha_{0}$ and  $ \beta_{0}$ we assume that either
\begin{enumerate}
\item[ {\rm (i)$_0$}]\quad  $\gamma_2(\alpha_{0},\; \beta_{0}) \ge 2-N$, or
\item[ {\rm (ii)$_0$}]\quad  $\gamma_1(\alpha_{0},\; \beta_{0}) \ge 2-N$ and 
$p \ge (2-\alpha_0)/(N-2)$.
\end{enumerate}
Then $u\equiv 0$ and $v\equiv  0$.
\end{teorema}

\begin{thebibliography}{99}


\bibitem{bly}  R. Benguria, S. Lorca and C. Yarur, Nonexistence of 
positive solution of semilinear elliptic equations, {\em Duke Math. J.} 
{\bf 74} (1994).

\bibitem{br} H. Br\`esis, Semilinear equations in $\Bbb{R}^N$ without 
condition at 
infinity, {\em Appl. Math. Optim.} {\bf 12} (1984). 

\bibitem{bv} H. Br\`esis and L. V\'eron, Removable singularities of some 
nonlinear elliptic equations, {\em Arch. Rational Mech. Anal.} {\bf 75}
(1980), 1-6. 

\bibitem{gmmy} M. Garc\'{\i}a-Huidobro , R. Man\'asevich R., E. Mitidieri 
and C. Yarur, 
Existence and Nonexistence of Positive Singular Solutions for
a  Class of Semilinear Elliptic Systems, to appear in {\em 
Arch. Rat. Mech \& Anal.}

\bibitem{gt} D. Gilbarg and N.S. Trudinger, Elliptic Partial Differential 
Equations of Second Order, Springer-Verlag, 1977.

\bibitem{flin} F.H. Lin, On the elliptic equation  $D_i[a_{ij}(x)D_jU] - 
k(x)U + K(x)U^p = 0$, {\em Proceedings of the AMS} {\bf 95} (1985), 219-226.  

\bibitem{ll} Ch. Lin and S. Lin, Positive radial solutions for 
$\Delta u + K(r)u^{(N+2)/(N-2)} = 0$ in $\Bbb{R}^N$ and related topics, 
{\em Appl. Anal.} {\bf 38} (1990), 121-159.  

\bibitem{ln} C. Loewner and L. Nirenberg, Partial differential equations under 
conformal or projective transformations, {\em Contributions to Analysis} 
Academic Press (1974), 245-272.

\bibitem{mi1} E. Mitidieri, 
Nonexistence of positive solutions of semilinear elliptic systems 
in $\Bbb{R}^N$, {\em Differential and Integral Equations,} 
in press.


\bibitem{mi2} E. Mitidieri, A Rellich type identity and 
applications, {\em Commun. in Partial Differential Equations} {\bf 18} (1993), 125-151.

\bibitem{na} M. Naito,  A note on bounded positive entire solutions 
of semilinear elliptic equations,{\em Hiroshima Math. J.} {\bf 14} (1984), 211-214.

\bibitem{ni} W.M. Ni,  On the elliptic equation 
$\Delta u + K(x)u^{(n+2)/(n-2)}=0,$ its generalizations, and
 applications in geometry, {\em Indiana Univ. Math. J. }{\bf 31} (1982),493-529.

\bibitem{sp} P. Pucci and J. Serrin,  A general variational identity, {\em Indiana University Mathematics Journal.}
{\bf 35} (1986), 681-703. 

\bibitem{sor} R. Soranzo, 
Isolated singularities of positive solutions of a superlinear 
biharmonic equation, {\em J.  Potential Theory,} in press.

\bibitem{sy} H. Soto and C. Yarur, Some existence results of semilinear elliptic 
equations, {\em Rendiconti di Matematica} {\bf 15} (1995),109-123.

\bibitem{vd} R.C.A.M. van der Vorst,  {\em Variational problems with a strongly indefinite structure.}  Doctoral thesis (1994). 

\bibitem{ve} L. V\'eron,  Comportement asymptotique des solutions 
d'\'equations elliptiques semi-lin\'eaire dans $\Bbb{R}^N$, {\em Ann.
Math. Pura Appl.}
{\bf 127} (1981), 25-50. 
\end{thebibliography}

{\sc Ceclia S. Yarur\newline
Departamento de Matem\'aticas \newline
Universidad de Santiago de Chile 
Casilla 307, Correo 2\newline
Santiago, Chile.\newline}
E-mail: cyarur@usach.cl
  




\end{document}

