Abstract. We prove the existence of nonconstant positive solutions for a system of the form −ε 2 ∆u + u = g(v), −ε 2 ∆v + v = f (u) in Ω, with Neumann boundary conditions on ∂Ω, where Ω is a smooth bounded domain and f , g are power-type nonlinearities having superlinear and subcritical growth at infinity. For small values of ε, the corresponding solutions uε and vε admit a unique maximum point which is located at the boundary of Ω.