code: def _hc(k, cs, rho, omega):
      return (
          cs
          / cupy.sin(omega)
          * (rho**k)
          * cupy.sin(omega * (k + 1))
          * cupy.greater(k, -1)
      )
  
imports: import cupy
  import cupyx.scipy.ndimage
  from cupyx.scipy.signal._iir_utils import apply_iir_sos
  from cupyx.scipy.signal._splines import _symiirorder1_nd, _symiirorder2_nd
  from cupyx.scipy.interpolate._bspline import BSpline
  import numpy
label_desc: Use image node_id 5 for calling a global function with example usage: _hc(k, cs, rho, omega) and returns: unknown