code: def test_beta():
      a, b, x = symbols("a b x", positive=True)
      e = x ** (a - 1) * (-x + 1) ** (b - 1) / beta(a, b)
      Q = QQ[a, b].get_field()
      h1 = expr_to_holonomic(e, x, domain=Q)
  
      _, Dx = DifferentialOperators(Q.old_poly_ring(x), "Dx")
      h2 = HolonomicFunction(
          (a + x * (-a - b + 2) - 1) + (x**2 - x) * Dx, x
      )
  
      assert h1 == h2
  
imports: from sympy.holonomic import DifferentialOperator, HolonomicFunction, DifferentialOperators, from_hyper, from_meijerg, expr_to_holonomic
  from sympy.holonomic.recurrence import RecurrenceOperators, HolonomicSequence
  from sympy.core import EulerGamma
  from sympy.core.numbers import I, Rational, pi
  from sympy.core.singleton import S
  from sympy.core.symbol import Symbol, symbols
  from sympy.functions.elementary.exponential import exp, log
  from sympy.functions.elementary.hyperbolic import asinh, cosh
  from sympy.functions.elementary.miscellaneous import sqrt
  from sympy.functions.elementary.trigonometric import cos, sin
  from sympy.functions.special.bessel import besselj
  from sympy.functions.special.beta_functions import beta
  from sympy.functions.special.error_functions import Ci, Si, erf, erfc
  from sympy.functions.special.gamma_functions import gamma
  from sympy.functions.special.hyper import hyper, meijerg
  from sympy.printing.str import sstr
  from sympy.series.order import O
  from sympy.simplify.hyperexpand import hyperexpand
  from sympy.polys.domains.integerring import ZZ
  from sympy.polys.domains.rationalfield import QQ
  from sympy.polys.domains.realfield import RR
label_desc: Use image node_id 22 for calling a global function with example usage: test_beta() without return types