code: def Wishart(symbol, n, scale_matrix):
      """
      Creates a random variable with Wishart Distribution.
  
      The density of the said distribution can be found at [1].
  
      Parameters
      ==========
  
      n: Positive Real number
          Represents degrees of freedom
      scale_matrix: Positive definite real square matrix
          Scale Matrix
  
      Returns
      =======
  
      RandomSymbol
  
      Examples
      ========
  
      >>> from sympy.stats import density, Wishart
      >>> from sympy import MatrixSymbol, symbols
      >>> n = symbols('n', positive=True)
      >>> W = Wishart('W', n, [[2, 1], [1, 2]])
      >>> X = MatrixSymbol('X', 2, 2)
      >>> density(W)(X).doit()
      exp(Trace(Matrix([
      [-1/3,  1/6],
      [ 1/6, -1/3]])*X))*Determinant(X)**(n/2 - 3/2)/(2**n*3**(n/2)*sqrt(pi)*gamma(n/2)*gamma(n/2 - 1/2))
      >>> density(W)([[1, 0], [0, 1]]).doit()
      exp(-2/3)/(2**n*3**(n/2)*sqrt(pi)*gamma(n/2)*gamma(n/2 - 1/2))
  
      References
      ==========
  
      .. [1] https://en.wikipedia.org/wiki/Wishart_distribution
  
      """
      if isinstance(scale_matrix, list):
          scale_matrix = ImmutableMatrix(scale_matrix)
      return rv(symbol, WishartDistribution, (n, scale_matrix))
  
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 3 for calling a global function with example usage: Wishart(symbol, n, scale_matrix) and returns: rv