def Morsi_Alexander(Re):
    """Calculates drag coefficient of a smooth sphere using the method in
    [1]_ as described in [2]_.

    .. math::
        C_D = \\left\\{ \\begin{array}{ll}
        \\frac{24}{Re} & \\mbox{if $Re < 0.1$}\\\\
        \\frac{22.73}{Re}+\\frac{0.0903}{Re^2} + 3.69 & \\mbox{if $0.1 < Re < 1$}\\\\
        \\frac{29.1667}{Re}-\\frac{3.8889}{Re^2} + 1.2220 & \\mbox{if $1 < Re < 10$}\\\\
        \\frac{46.5}{Re}-\\frac{116.67}{Re^2} + 0.6167 & \\mbox{if $10 < Re < 100$}\\\\
        \\frac{98.33}{Re}-\\frac{2778}{Re^2} + 0.3644 & \\mbox{if $100 < Re < 1000$}\\\\
        \\frac{148.62}{Re}-\\frac{4.75\\times10^4}{Re^2} + 0.3570 & \\mbox{if $1000 < Re < 5000$}\\\\
        \\frac{-490.5460}{Re}+\\frac{57.87\\times10^4}{Re^2} + 0.46 & \\mbox{if $5000 < Re < 10000$}\\\\
        \\frac{-1662.5}{Re}+\\frac{5.4167\\times10^6}{Re^2} + 0.5191 & \\mbox{if $10000 < Re < 50000$}\\end{array} \\right.

    Parameters
    ----------
    Re : float
        Particle Reynolds number of the sphere using the surrounding fluid
        density and viscosity, [-]

    Returns
    -------
    Cd : float
        Drag coefficient [-]

    Notes
    -----
    Range is Re <= 2E5.
    Original was reviewed, and confirmed to contain the cited equations.

    Examples
    --------
    >>> Morsi_Alexander(200)
    0.7866

    References
    ----------
    .. [1] Morsi, S. A., and A. J. Alexander. "An Investigation of Particle
       Trajectories in Two-Phase Flow Systems." Journal of Fluid Mechanics
       55, no. 02 (September 1972): 193-208. doi:10.1017/S0022112072001806.
    .. [2] Barati, Reza, Seyed Ali Akbar Salehi Neyshabouri, and Goodarz
       Ahmadi. "Development of Empirical Models with High Accuracy for
       Estimation of Drag Coefficient of Flow around a Smooth Sphere: An
       Evolutionary Approach." Powder Technology 257 (May 2014): 11-19.
       doi:10.1016/j.powtec.2014.02.045.
    """
    if Re < 0.1:
        return 24.0 / Re
    elif Re < 1:
        return 22.73 / Re + 0.0903 / Re ** 2 + 3.69
    elif Re < 10:
        return 29.1667 / Re - 3.8889 / Re ** 2 + 1.222
    elif Re < 100:
        return 46.5 / Re - 116.67 / Re ** 2 + 0.6167
    elif Re < 1000:
        return 98.33 / Re - 2778.0 / Re ** 2 + 0.3644
    elif Re < 5000:
        return 148.62 / Re - 47500.0 / Re ** 2 + 0.357
    elif Re < 10000:
        return -490.546 / Re + 578700.0 / Re ** 2 + 0.46
    else:
        return -1662.5 / Re + 5416700.0 / Re ** 2 + 0.5191
