code: def robust_std(x, debug=False):
      
  
      x = numpy.asarray(x)
  
      # compute percentiles and robust estimator
      q25 = numpy.percentile(x, 25)
      q75 = numpy.percentile(x, 75)
      sigmag = 0.7413 * (q75 - q25)
  
      if debug:
          print('debug|sigmag -> q25......................:', q25)
          print('debug|sigmag -> q75......................:', q75)
          print('debug|sigmag -> Robust standard deviation:', sigmag)
  
      return sigmag
docs: Compute a robust estimator of the standard deviation
  
      See Eq. 3.36 (page 84) in Statistics, Data Mining, and Machine
      in Astronomy, by Ivezic, Connolly, VanderPlas & Gray
  
      Parameters
      ----------
      x : 1d numpy array, float
          Array of input values which standard deviation is requested.
      debug : bool
          If True prints computed values
  
      Returns
      -------
      sigmag : float
          Robust estimator of the standar deviation