The existence of positive solutions for a nonlocal boundary-value problem with vector-valued response is investigated. We develop duality and variational principles for this problem. Our variational approach enables us to approximate solutions and give a measure of a duality gap between the primal and dual functional for minimizing sequences.