The Neumann boundary value problem is examined for systems of elliptic equations of the form $\Delta u + g(u) = f(x), x \in \omega.$ It is assumed that $g \in 2 C(\mathbb(R)^N,\mathbb(R)^N)$ is a bounded function which may vanish at infinity. Leray-Schauder degree methods are used.