code: def squared_radial_component(x, tol=0.01):
      """Part of the GQI2 integral
  
      Eq.8 in the referenced paper by Yeh et al. 2010
      """
      with warnings.catch_warnings():
          warnings.simplefilter("ignore")
          result = (2 * x * np.cos(x) + (x * x - 2) * np.sin(x)) / (
              x**3
          )
      x_near_zero = (x < tol) & (x > -tol)
      return np.where(x_near_zero, 1.0 / 3, result)
  
imports: import numpy
  from dipy.reconst.odf import OdfModel, OdfFit
  from dipy.reconst.cache import Cache
  import warnings
  from dipy.reconst.multi_voxel import multi_voxel_fit
label_desc: Use image node_id 2 for calling a global function with example usage: squared_radial_component(x, tol) and returns: np