Make Fluxion with the distance between body i and j

def flux_r(q_vars: List[fl.Var], i: int, j: int):
    """Make Fluxion with the distance between body i and j"""
    return fl.sqrt(flux_r2(q_vars, i, j))