code: def eigval_decomp(sym_array):
      """
      Returns
      -------
      W: array of eigenvectors
      eigva: list of eigenvalues
      k: largest eigenvector
      """
      # check if symmetric, do not include shock period
      eigva, W = linalg.eig(sym_array, left=True, right=False)
      k = np.argmax(eigva)
      return W, eigva, k
  
file_path: statsmodels/statsmodels/tsa/vector_ar/util.py
imports: from statsmodels.compat.pandas import frequencies
  from statsmodels.compat.python import asbytes
  from statsmodels.tools.validation import array_like, int_like
  import numpy
  import pandas
  from scipy import stats, linalg
  import statsmodels.tsa.tsatools
label_desc: Use image node_id 10 for calling a global function with example usage: eigval_decomp(sym_array) and returns: W,  eigva,  k