code: def test_sample_pymc():
      distribs_pymc = [
          MatrixNormal(
              "M", [[5, 6], [3, 4]], [[1, 0], [0, 1]], [[2, 1], [1, 2]]
          ),
          Wishart("W", 7, [[2, 1], [1, 2]]),
      ]
      size = 3
      pymc = import_module("pymc")
      if not pymc:
          skip("PyMC is not installed. Abort tests for _sample_pymc.")
      else:
          for X in distribs_pymc:
              samps = sample(X, size=size, library="pymc")
              for sam in samps:
                  assert Matrix(sam) in X.pspace.distribution.set
          M = MatrixGamma("M", 1, 2, [[1, 0], [0, 1]])
          raises(NotImplementedError, lambda: sample(M, size=3))
  
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 7 for calling a global function with example usage: test_sample_pymc() without return types