docs: def timefreq(x,fs=200):
      """
      TIMEFREQ
      
      This function takes the time series and the sampling rate and calculates the
      total number of points, the maximum frequency, the minimum (or change in)
      frequency, and the vector of frequency points F.
      
      Version: 2011may04
      """
      from numpy import size, shape, arange, append
      
      maxfreq=float(fs)/2.0 # Maximum frequency
      minfreq=float(fs)/float(size(x,0)) # Minimum and delta frequency -- simply the inverse of the length of the recording in seconds
      F=arange(minfreq,maxfreq+minfreq,minfreq) # Create frequencies evenly spaced from 0:minfreq:maxfreq
      F=append(0,F) # Add zero-frequency component
      
      return F