Computed as exp(lbeta(a, b)).
beta(a, b)
# S4 method for class 'numeric,numeric'
beta(a, b)
# S4 method for class 'dualr,dualr'
beta(a, b)
# S4 method for class 'dualr,numeric'
beta(a, b)
# S4 method for class 'numeric,dualr'
beta(a, b)beta(a, b) with derivative.
beta(2, 3) # 1/12
#> [1] 0.08333333
a <- dual_variable(2)
value(beta(a, 3))
#> [1] 0.08333333