def METHOD_NAME(
    R_0: float, A: float, B_0: float, I_p: float, kappa: float, delta: float
) -> float:
    """
    Estimate safety factor at the 95th percentile flux surface based on an empirical fit.
    Parameters
    ----------
    R_0:
        Plasma major radius [m]
    A:
        Plasma aspect ratio
    B_0:
        Toroidal field at major radius [T]
    I_p:
        Plasma current [A]
    kappa:
        Plasma elongation
    delta:
        Plasma triangularity
    Notes
    -----
    Uckan et al., ITER Physics Design Guidelines, 1989, sec. 2.3
    https://inis.iaea.org/search/search.aspx?orig_q=RN:21068960
    Ref [11] in the above does not appear to include the geometry factor
    """
    eps = 1 / A
    geometry_factor = (1.17 - 0.65 * eps) / (1 - eps**2) ** 2
    q_star = calc_qstar_uckan(R_0, A, B_0, I_p, kappa, delta)
    return q_star * geometry_factor