code: def triple_odf_maxima(vertices, odf, width):
      indmax1 = np.argmax([odf[i] for i, v in enumerate(vertices)])
      odfmax1 = odf[indmax1]
      pole = vertices[indmax1]
      eqvert = equatorial_zone_vertices(vertices, pole, width)
      indmax2, odfmax2 = equatorial_maximum(vertices, odf, pole, width)
      indmax3 = eqvert[
          np.argmin(
              [
                  np.abs(np.dot(vertices[indmax2], vertices[p]))
                  for p in eqvert
              ]
          )
      ]
      odfmax3 = odf[indmax3]
      """
      cross12 = np.cross(vertices[indmax1],vertices[indmax2])
      cross12 = cross12/np.sqrt(np.sum(cross12**2))
      indmax3, odfmax3 = patch_maximum(vertices, odf, cross12, 2*width)
      """
      return [
          (indmax1, odfmax1),
          (indmax2, odfmax2),
          (indmax3, odfmax3),
      ]
  
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 12 for calling a global function with example usage: triple_odf_maxima(vertices, odf, width) and returns: list