R语言生成随机变量

x是（0，1）上的随机变量，服从参数为（a1，a2）的贝塔分布，R语言如何用程序实现

rbeta(n, shape1, shape2, ncp = 0)返回服从beta分布的随机数
n：产生随机数的个数
shape1, shape2：beta分布的两个参数
ncp：non-centrality parameter非中心化参数，对于标准beta分布取ncp=0

标准beta分布的取值范围为[0,1]
参数shape1=a、shape2=b的Beta分布的密度函数为：Γ(a+b)/(Γ(a)*Γ(b))*x^(a-1)*(1-x)^(b-1)

rbeta(n, shape1, shape2, ncp = 0)

rbeta(n, a1, a2)