docs: def Morsi_Alexander(Re):
      r'''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
          Reynolds number of the sphere, [-]
  
      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:
          Cd = 24./Re
      elif Re < 1:
          Cd = 22.73/Re + 0.0903/Re**2 + 3.69
      elif Re < 10:
          Cd = 29.1667/Re - 3.8889/Re**2 + 1.222
      elif Re < 100:
          Cd = 46.5/Re - 116.67/Re**2 + 0.6167
      elif Re < 1000:
          Cd = 98.33/Re - 2778./Re**2 + 0.3644
      elif Re < 5000:
          Cd = 148.62/Re - 4.75E4/Re**2 + 0.357
      elif Re < 10000:
          Cd = -490.546/Re + 57.87E4/Re**2 + 0.46
      else:
          Cd = -1662.5/Re + 5.4167E6/Re**2 + 0.5191
      return Cd