code: def MatrixStudentT(
      symbol, nu, location_matrix, scale_matrix_1, scale_matrix_2
  ):
      """
      Creates a random variable with Matrix Gamma Distribution.
  
      The density of the said distribution can be found at [1].
  
      Parameters
      ==========
  
      nu: Positive Real number
          degrees of freedom
      location_matrix: Positive definite real square matrix
          Location Matrix of shape ``n x p``
      scale_matrix_1: Positive definite real square matrix
          Scale Matrix of shape ``p x p``
      scale_matrix_2: Positive definite real square matrix
          Scale Matrix of shape ``n x n``
  
      Returns
      =======
  
      RandomSymbol
  
      Examples
      ========
  
      >>> from sympy import MatrixSymbol,symbols
      >>> from sympy.stats import density, MatrixStudentT
      >>> v = symbols('v',positive=True)
      >>> M = MatrixStudentT('M', v, [[1, 2]], [[1, 0], [0, 1]], [1])
      >>> X = MatrixSymbol('X', 1, 2)
      >>> density(M)(X)
      gamma(v/2 + 1)*Determinant((Matrix([[-1, -2]]) + X)*(Matrix([
      [-1],
      [-2]]) + X.T) + Matrix([[1]]))**(-v/2 - 1)/(pi**1.0*gamma(v/2)*Determinant(Matrix([[1]]))**1.0*Determinant(Matrix([
      [1, 0],
      [0, 1]]))**0.5)
  
      References
      ==========
  
      .. [1] https://en.wikipedia.org/wiki/Matrix_t-distribution
  
      """
      if isinstance(location_matrix, list):
          location_matrix = ImmutableMatrix(location_matrix)
      if isinstance(scale_matrix_1, list):
          scale_matrix_1 = ImmutableMatrix(scale_matrix_1)
      if isinstance(scale_matrix_2, list):
          scale_matrix_2 = ImmutableMatrix(scale_matrix_2)
      args = (nu, location_matrix, scale_matrix_1, scale_matrix_2)
      return rv(symbol, MatrixStudentTDistribution, args)
  
imports: from math import prod
  from sympy.core.basic import Basic
  from sympy.core.numbers import pi
  from sympy.core.singleton import S
  from sympy.functions.elementary.exponential import exp
  from sympy.functions.special.gamma_functions import multigamma
  from sympy.core.sympify import sympify, _sympify
  from sympy.matrices import ImmutableMatrix, Inverse, Trace, Determinant, MatrixSymbol, MatrixBase, Transpose, MatrixSet, matrix2numpy
  from sympy.stats.rv import _value_check, RandomMatrixSymbol, NamedArgsMixin, PSpace, _symbol_converter, MatrixDomain, Distribution
  from sympy.external import import_module
label_desc: Use image node_id 5 for calling a global function with example usage: MatrixStudentT(symbol, nu, location_matrix, scale_matrix_1, scale_matrix_2) and returns: rv