code: def quadratic_interpolate_peak(left, middle, right):
      
      bin_offset = 1.0/2.0 * (left - right) / (left - 2 * middle + right)
      peak_value = middle + 0.25 * (left - right) * bin_offset
      return bin_offset, peak_value
docs: Interpolate the peak and offset using a quadratic approximation
  
      Parameters
      ----------
      left : numpy array
          Values at a relative bin value of [-1]
      middle : numpy array
          Values at a relative bin value of [0]
      right : numpy array
          Values at a relative bin value of [1]
  
      Returns
      -------
      bin_offset : numpy array
          Array of bins offsets, each in the range [-1/2, 1/2]
      peak_values : numpy array
          Array of the estimated peak values at the interpolated offset