code: def test_hist_color_multi_band(image_array_2bands):
      """Test that multiple colors work for multiband images."""
      f, ax = ep.hist(image_array_2bands, colors=["red", "blue"])
      colors = [
          a.patches[0].__dict__.get("_original_facecolor") for a in ax
      ]
      expected_colors = [
          np.array([1.0, 0.0, 0.0, 1.0]),
          np.array([0.0, 0.0, 1.0, 1.0]),
      ]
      for i in range(2):
          assert np.array_equal(colors[i], expected_colors[i])
      plt.close(f)
  
imports: import numpy
  import pytest
  import matplotlib.pyplot
  import earthpy.plot
label_desc: Use image node_id 10 for calling a global function with example usage: test_hist_color_multi_band(image_array_2bands) without return types