def r_to_rq(r):
    """
    reduced coordinate inverse:                                   
                                                        
               (1+rq)
           r = ------
               (1-rq)
    
    Hernquist & Ostriker 1992 eq. 2.16
    """
    return (r-1.0)/(r+1.0);