code: def test_polynomial1():
      # Polynomial part of Gaussian densities are Hermite polynomials.
      for dim in range(1, 10):
          q = rft.Gaussian().quasi(dim)
          h = hermitenorm(dim - 1)
          assert_almost_equal(q.c, h.c)
  
imports: import numpy
  import scipy.stats
  from scipy.special import gammaln, hermitenorm
  import pytest
  from numpy.testing import assert_almost_equal
  from ..None import rft
label_desc: Use image node_id 6 for calling a global function with example usage: test_polynomial1() without return types