code: def elimination_matrix(n):
      """
      Create the elimination matrix L_n which satisfies vech(M) = L_n vec(M) for
      any matrix M
  
      Parameters
      ----------
  
      Returns
      -------
      """
      n = int_like(n, "n")
      vech_indices = vec(np.tril(np.ones((n, n))))
      return np.eye(n * n)[vech_indices != 0]
  
imports: from __future__ import annotations
  from statsmodels.compat.python import Literal, lrange
  import warnings
  import numpy
  import pandas
  from pandas import DataFrame
  from pandas.tseries import offsets
  from pandas.tseries.frequencies import to_offset
  from statsmodels.tools.data import _is_recarray, _is_using_pandas
  from statsmodels.tools.sm_exceptions import ValueWarning
  from statsmodels.tools.typing import NDArray
  from statsmodels.tools.validation import array_like, bool_like, int_like, string_like
label_desc: Use image node_id 14 for calling a global function with example usage: elimination_matrix(n) and returns: np