code: def test_zipfian_asymptotic(self):
      # test limiting case that zipfian(a, n) -> zipf(a) as n-> oo
      a = 6.5
      N = 10000000
      k = np.arange(1, 21)
      assert_allclose(zipfian.pmf(k, a, N), zipf.pmf(k, a))
      assert_allclose(zipfian.cdf(k, a, N), zipf.cdf(k, a))
      assert_allclose(zipfian.sf(k, a, N), zipf.sf(k, a))
      assert_allclose(
          zipfian.stats(a, N, moments="msvk"),
          zipf.stats(a, moments="msvk"),
      )
  
imports: import pytest
  import itertools
  from scipy.stats import betabinom, betanbinom, hypergeom, nhypergeom, bernoulli, boltzmann, skellam, zipf, zipfian, binom, nbinom, nchypergeom_fisher, nchypergeom_wallenius, randint
  import numpy
  from numpy.testing import assert_almost_equal, assert_equal, assert_allclose, suppress_warnings
  from scipy.special import binom
  from scipy.optimize import root_scalar
  from scipy.integrate import quad
label_desc: Use image node_id 1 for calling the TestZipfian obj's underlying member method code with example usage: obj.test_zipfian_asymptotic() without return types