code: def test_normal(self):
      # with two cumulants, it's just a gaussian
      ne2 = ExpandedNormal([3, 4])
      x = np.linspace(-2.0, 2.0, 100)
      assert_allclose(ne2.pdf(x), stats.norm.pdf(x, loc=3, scale=2))
  
imports: import warnings
  import numpy
  from numpy.testing import assert_equal, assert_raises, assert_allclose
  import numpy.testing
  from scipy.special import gamma, factorial, factorial2
  import scipy.stats
  from statsmodels.distributions.edgeworth import _faa_di_bruno_partitions, cumulant_from_moments, ExpandedNormal
label_desc: Use image node_id 3 for calling the TestExpandedNormal obj's underlying member method code with example usage: obj.test_normal() without return types