code: def test_expm_multiply_interval_matrix(self):
      np.random.seed(1234)
      start = 0.1
      stop = 3.2
      endpoint = True
      for num in (14, 13, 2):
          for n in (1, 2, 5, 20, 40):
              for k in (1, 2):
                  A = scipy.linalg.inv(np.random.randn(n, n))
                  B = np.random.randn(n, k)
                  X = expm_multiply(
                      A,
                      B,
                      start=start,
                      stop=stop,
                      num=num,
                      endpoint=endpoint,
                  )
                  samples = np.linspace(
                      start=start, stop=stop, num=num, endpoint=endpoint
                  )
                  for solution, t in zip(X, samples):
                      assert_allclose(
                          solution, scipy.linalg.expm(t * A).dot(B)
                      )
  
imports: import numpy
  from numpy.testing import assert_allclose, assert_, assert_equal, suppress_warnings
  from scipy.sparse import SparseEfficiencyWarning
  import scipy.linalg
  from scipy.sparse.linalg._expm_multiply import _theta, _compute_p_max, _onenormest_matrix_power, expm_multiply, _expm_multiply_simple, _expm_multiply_interval
label_desc: Use image node_id 3 for calling the TestExpmActionInterval obj's underlying member method code with example usage: obj.test_expm_multiply_interval_matrix() without return types