code: def quantile(arg, quantile, interpolation='linear'):
      
      if isinstance(quantile, collections.abc.Sequence):
          op = ops.MultiQuantile(arg, quantile, interpolation)
      else:
          op = ops.Quantile(arg, quantile, interpolation)
      return op.to_expr()
docs: Return value at the given quantile, a la numpy.percentile.
  
      Parameters
      ----------
      quantile : float/int or array-like
          0 <= quantile <= 1, the quantile(s) to compute
      interpolation : {'linear', 'lower', 'higher', 'midpoint', 'nearest'}
  
          This optional parameter specifies the interpolation method to use,
          when the desired quantile lies between two data points `i` and `j`:
  
          * linear: `i + (j - i) * fraction`, where `fraction` is the
            fractional part of the index surrounded by `i` and `j`.
          * lower: `i`.
          * higher: `j`.
          * nearest: `i` or `j` whichever is nearest.
          * midpoint: (`i` + `j`) / 2.
  
      Returns
      -------
      quantile
          if scalar input, scalar type, same as input
          if array input, list of scalar type