  We use a recently proposed diagonalization/Monte Carlo computational scheme
to study relativistic two-body and three-body bound states in 1+1 dimensional
phi^6-phi^4 theory. We find that the approach is well-suited for calculating
bound state energies and wavefunctions.
