code: def max_frequency(signal, fs):
      """Computes maximum frequency of the signal.
  
      Feature computational cost: 2
  
      Parameters
      ----------
      signal : nd-array
          Input from which maximum frequency is computed
      fs : float
          Sampling frequency
  
      Returns
      -------
      float
          0.95 of maximum frequency using cumsum
      """
      f, fmag = calc_fft(signal, fs)
      cum_fmag = np.cumsum(fmag)
  
      try:
          ind_mag = np.where(cum_fmag > cum_fmag[-1] * 0.95)[0][0]
      except IndexError:
          ind_mag = np.argmax(cum_fmag)
  
      return f[ind_mag]
  
imports: import scipy.signal
  from tsfel.feature_extraction.features_utils import *
label_desc: Use image node_id 39 for calling a global function with example usage: max_frequency(signal, fs) and returns: f