Abstract. We consider a system of the form −ε 2 ∆u + u = g(v), −ε 2 ∆v + v = f (u) in Ω with Dirichlet boundary condition on ∂Ω, where Ω is a smooth bounded domain in R N , N 3 and f, g are power-type nonlinearities having superlinear and subcritical growth at infinity. We prove that the least energy solutions to such a system concentrate, as ε goes to zero, at a point of Ω which maximizes the distance to the boundary of Ω.2000 Mathematics Subject Classification: 35J50, 35J55, 38E05.