  In present work, we study an numerical approach to one dimensional finite
volume three-body interaction, the method is demonstrated by considering a toy
model of three spinless particles interacting with pair-wise $\delta$-function
potentials. The numerical results are compared with the exact solutions of
three spinless bosons interaction when strength of short-range interactions are
set equal for all pairs.
