Accurate numerical solution of the five-body Schrodinger equation is effected via variational Monte Carlo in a correlated quark model. The spectrum is assumed to exhibit a narrow resonance with strangeness S = +1. A fully antisymmetrized and pair-correlated five-quark wave function is obtained for the assumed non-relativistic Hamiltonian which has spin, isospin, and color dependent pair interactions and many-body confining terms which are fixed by the non-exotic spectra. Gauge field dynamics are modeled via flux tube exchange factors. The energy determined for the ground states with spin-parity = 1/2 -(1/2+) is 2.22 GeV (2.50 GeV). A lower energy negative parity state is consistent with recent lattice results.