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