code: def test_sample_seed():
      X = MatrixNormal(
          "M", [[5, 6], [3, 4]], [[1, 0], [0, 1]], [[2, 1], [1, 2]]
      )
  
      libraries = ["scipy", "numpy", "pymc"]
      for lib in libraries:
          try:
              imported_lib = import_module(lib)
              if imported_lib:
                  s0, s1, s2 = [], [], []
                  s0 = sample(X, size=10, library=lib, seed=0)
                  s1 = sample(X, size=10, library=lib, seed=0)
                  s2 = sample(X, size=10, library=lib, seed=1)
                  for i in range(10):
                      assert (s0[i] == s1[i]).all()
                      assert (s1[i] != s2[i]).all()
  
          except NotImplementedError:
              continue
  
file_path: sympy/sympy/stats/tests/test_matrix_distributions.py
imports: from sympy.concrete.products import Product
  from sympy.core.numbers import pi
  from sympy.core.singleton import S
  from sympy.core.symbol import Dummy, symbols
  from sympy.functions.elementary.exponential import exp
  from sympy.functions.elementary.miscellaneous import sqrt
  from sympy.functions.special.gamma_functions import gamma
  from sympy.matrices import Determinant, Matrix, Trace, MatrixSymbol, MatrixSet
  from sympy.stats import density, sample
  from sympy.stats.matrix_distributions import MatrixGammaDistribution, MatrixGamma, MatrixPSpace, Wishart, MatrixNormal, MatrixStudentT
  from sympy.testing.pytest import raises, skip
  from sympy.external import import_module
label_desc: Use image node_id 8 for calling a global function with example usage: test_sample_seed() without return types