dbeta(x, shape1, shape2) pbeta(q, shape1, shape2, ncp=0) qbeta(p, shape1, shape2) rbeta(n, shape1, shape2)
The algorithm for computing noncentrality in pbeta is described in Posten (1993).
Posten, H. O. (1993). An effective algorithm for the noncentral beta distribution function. The American Statistician 47 129-131.
rbeta(20,2,3) #sample of 20 with shape parameters 2 and 3