code: def test_T1():
      # O-dim EC density should be tail probality.
      x = np.linspace(0.1, 10, 100)
      for dfd in [40, 50]:
          t = rft.TStat(dfd=dfd)
          assert_almost_equal(t(x), scipy.stats.t.sf(x, dfd))
      t = rft.TStat(dfd=np.inf)
      assert_almost_equal(t(x), scipy.stats.norm.sf(x))
  
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 12 for calling a global function with example usage: test_T1() without return types