docs: def stdev(data, xbar=None):
      """Return the square root of the sample variance.
  
      See ``variance`` for arguments and other details.
  
      """
      var = variance(data, xbar)
      try:
          return var.sqrt()
      except AttributeError:
          return math.sqrt(var)