code: def patch_maximum(vertices, odf, pole, width):
      eqvert = patch_vertices(vertices, pole, width)
      # need to test for whether eqvert is empty or not
      if len(eqvert) == 0:
          print(
              "empty cone around pole %s with with width %f"
              % (np.array_str(pole), width)
          )
          return np.Null, np.Null
      eqvals = [odf[i] for i in eqvert]
      eqargmax = np.argmax(eqvals)
      eqvertmax = eqvert[eqargmax]
      eqvalmax = eqvals[eqargmax]
      return eqvertmax, eqvalmax
  
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 9 for calling a global function with example usage: patch_maximum(vertices, odf, pole, width) and returns: eqvertmax,  eqvalmax, np,  np