The theoretical underpinnings of the proposed algorithm design were developed extensively, with complete implementations in the diatomic case and initial implementation in the more complex triatomic case. The treatment of a finite region Schrodinger equation proved much more difficult than anticipated; in particular, the boundary conditions required was unexpected, but yielded great insight into the heavily utilised R-matrix technique. We found a solid choice of basis set for this problem, implemented it within the DVR3D program and achieved promising results for model systems and HOCl. Two publications have so far resulted, with additional publications in preparation. Further, this work was presented at international conferences.