now using spherical HO with Laguerre from GSL (previous function errorneously gave zeros for small r and large n)