code: def test_jacobi_precise():
      x, w = gauss_jacobi(3, Rational(-1, 2), S.Half, 40)
      assert [str(r) for r in x] == [
          "-0.6234898018587335305250048840042398106323",
          "0.2225209339563144042889025644967947594664",
          "0.9009688679024191262361023195074450511659",
      ]
      assert [str(r) for r in w] == [
          "0.3379547635663543330553835737094171534907",
          "1.097332224279111467485302294320899710461",
          "1.706305665744327437921957515249186020246",
      ]
  
      x, w = gauss_jacobi(3, 2, 3, 40)
      assert [str(r) for r in x] == [
          "-0.4611587037808976179121958105554375981274",
          "0.1043853303832390210914918407615869143233",
          "0.6295006461249313240934312425211234110769",
      ]
      assert [str(r) for r in w] == [
          "0.1793761350221326596137764371503859752628",
          "0.6159564099114715430909548532229749439714",
          "0.2713341217330624639619353762933057474325",
      ]
  
file_path: sympy/sympy/integrals/tests/test_quadrature.py
imports: from sympy.core import S, Rational
  from sympy.integrals.quadrature import gauss_legendre, gauss_laguerre, gauss_hermite, gauss_gen_laguerre, gauss_chebyshev_t, gauss_chebyshev_u, gauss_jacobi, gauss_lobatto
label_desc: Use image node_id 14 for calling a global function with example usage: test_jacobi_precise() without return types