A variational R-matrix finite element procedure for solving ultra-cold collision problems