def Wood_1966(Re, eD):
    """Calculates Darcy friction factor using the method in Wood (1966) [2]_
    as shown in [1]_.

    .. math::
        f_d = 0.094(\\frac{\\epsilon}{D})^{0.225} + 0.53(\\frac{\\epsilon}{D})
        + 88(\\frac{\\epsilon}{D})^{0.4}Re^{-A_1}

    .. math::
        A_1 = 1.62(\\frac{\\epsilon}{D})^{0.134}

    Parameters
    ----------
    Re : float
        Reynolds number, [-]
    eD : float
        Relative roughness, [-]

    Returns
    -------
    fd : float
        Darcy friction factor [-]

    Notes
    -----
    Range is 4E3 <= Re <= 5E7;  1E-5 <= eD <= 4E-2.

    Examples
    --------
    >>> Wood_1966(1E5, 1E-4)
    0.021587570560090762

    References
    ----------
    .. [1] Winning, H. and T. Coole. "Explicit Friction Factor Accuracy and
       Computational Efficiency for Turbulent Flow in Pipes." Flow, Turbulence
       and Combustion 90, no. 1 (January 1, 2013): 1-27.
       doi:10.1007/s10494-012-9419-7
    .. [2] 	Wood, D.J.: An Explicit Friction Factor Relationship, vol. 60.
       Civil Engineering American Society of Civil Engineers (1966)
    """
    A1 = 1.62 * eD ** 0.134
    return 0.094 * eD ** 0.225 + 0.53 * eD + 88 * eD ** 0.4 * Re ** -A1
