The paper deals with the existence and multiplicity of positive solutions for the system of nonlinear Hammerstein integral equations u(x) = ∫ 1 0 k(x, y)f1(y, u(y), v(y), w(y))dy, v(x) = ∫ 1 0 k(x, y)f2(y, u(y), v(y), w(y))dy, w(x) = ∫ 1 0 k(x, y)f3(y, u(y), v(y), w(y))dy. We use concave functions to characterize growing and interacting behaviors of our nonlinearities so that f1, f2, f3 cover three cases: the first with all superlinear, the second with all sublinear, and the last with two superlinear and the other sublinear. Based on a priori estimates achieved by using Jensen’s integral inequality, we use fixed point index theory to establish our main results. As an application, we use our main results to establish the existence and multiplicity of positive solutions for a system of nth order boundary value problems for nonlinear ordinary differential equations. MSC : 45G15, 45M20, 47H07, 47H11, 34A34, 34B18