def _f(x): return mp.gamma(g+1)* \ np.power(x/mpf(2), -g) * \ mp.exp(1j*x) * mp.besselj(g, x)
def test_J6(): assert mpmath.besselj(2, 1 + 1j).ae(mpc('0.04157988694396212', '0.24739764151330632'))
def test_J6(): assert mpmath.besselj(2, 1 + 1j).ae( mpc('0.04157988694396212', '0.24739764151330632'))