Use uniroot() to find roots and compare with bisection outcome
uniroot.beta(x1, x2, alpha, par.list, R, p.mat, k)