code: def test_nhypergeom_accuracy():
      # Check that nhypergeom.rvs post-gh-13431 gives the same values as
      # inverse transform sampling
      np.random.seed(0)
      x = nhypergeom.rvs(22, 7, 11, size=100)
      np.random.seed(0)
      p = np.random.uniform(size=100)
      y = nhypergeom.ppf(p, 22, 7, 11)
      assert_equal(x, y)
  
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 8 for calling a global function with example usage: test_nhypergeom_accuracy() without return types