code: def create_traj_ranges(start, stop, N):
      
      steps = (1.0/(N-1)) * (stop - start)
      if np.isscalar(steps):
          return steps*np.arange(N) + start
      else:
          return steps[:, None]*np.arange(N) + start[:, None]
docs: Fills in the trajectory range.
  
      # Adapted from https://stackoverflow.com/a/40624614