Abstract: | In this paper we consider a system of heat equations ut = Δu, vt = Δv in an unbounded domain Ω??N coupled through the Neumann boundary conditions uv = vp, vv = uq, where p>0, q>0, pq>1 and ν is the exterior unit normal on ?Ω. It is shown that for several types of domain there exists a critical exponent such that all of positive solutions blow up in a finite time in subcritical case (including the critical case) while there exist positive global solutions in the supercritical case if initial data are small. |