code: def test_visualise_gradient_table_G_Delta_rainbow():
      gtab_4d = generate_gtab4D()
      qtdmri.visualise_gradient_table_G_Delta_rainbow(gtab_4d)
  
      gtab_4d.small_delta[
          4
      ] += 0.001  # so now the gtab has multiple small_delta
      assert_raises(
          ValueError,
          qtdmri.visualise_gradient_table_G_Delta_rainbow,
          gtab_4d,
      )
  
imports: import warnings
  import numpy
  import scipy.integrate
  from numpy.testing import assert_, assert_almost_equal, assert_array_almost_equal, assert_equal, assert_raises
  import pytest
  from dipy.core.gradients import gradient_table_from_qvals_bvecs
  from dipy.data import get_gtab_taiwan_dsi, get_sphere
  from dipy.reconst import qtdmri, mapmri
  from dipy.reconst.shm import descoteaux07_legacy_msg
  from dipy.sims.voxel import multi_tensor, add_noise
  from dipy.testing.decorators import set_random_number_generator
label_desc: Use image node_id 21 for calling a global function with example usage: test_visualise_gradient_table_G_Delta_rainbow() without return types