code: def fhti(rmin, rmax, n, q, mu):
      r
  
      # Central point log10(r_c) of periodic interval
      logrc = (rmin + rmax)/2
  
      # Central index (1/2 integral if n is even)
      nc = (n + 1)/2.
  
      # Log spacing of points
      dlogr = (rmax - rmin)/n
      dlnr = dlogr*np.log(10.)
  
      # Get low-ringing kr
      y = 1j*np.pi/(2.0*dlnr)
      zp = special.loggamma((mu + 1.0 + q)/2.0 + y)
      zm = special.loggamma((mu + 1.0 - q)/2.0 + y)
      arg = np.log(2.0)/dlnr + (zp.imag + zm.imag)/np.pi
      kr = np.exp((arg - np.round(arg))*dlnr)
  
      # Calculate required input x-values (freq); angular freq -> freq
      freq = 10**(logrc + (np.arange(1, n+1) - nc)*dlogr)/(2*np.pi)
  
      # Calculate tcalc with adjusted kr
      logkc = np.log10(kr) - logrc
      tcalc = 10**(logkc + (np.arange(1, n+1) - nc)*dlogr)
  
      # rk = r_c/k_r; adjust for Fourier transform scaling
      rk = 10**(logrc - logkc)*np.pi/2
  
      return freq, tcalc, dlnr, kr, rk
docs: r"""Return parameters required for FFTLog.