已知概率密度函数，求随机变量问题

wangjun98989898

如题，已知概率密度函数，拟依据其生成随机变量，想法是先对其积分求累计概率分布，再求反函数，然后由（0，1）均匀分布随机数来生成随机变量，但是，其累计概率密度函数是一个含有hypergeom（超几何分布函数）的函数，怎么办呢？？？
4 ^: g1 T! g. w9 J1 A 具体问题（符合贝塔分布的一个概率密度函数）：
% d. F8 L% s3 G) K+ d  b* Ef=int((1/216815.2)*((x/196000)^-0.05)*(1-(x/196000))^-0.05)；MATLAB运行结果是：
7 B; Z% k# ]6 `% Y2 k* h; q$ Xf =

                      2374999999999999993
                      -------------------
" J/ L# x! O8 B7 v                  -5  2500000000000000000
   0.8928931123 10   x                    hypergeom(
/ f, w) S) i. C9 j& P+ L7 ^  A$ `                                                                      -5
         [0.05000000000, 0.9500000000], [1.950000000], 0.5102040816 10   x)
6 p) k5 w2 d# Z; d5 u- H* Q- ?当我继续运行
' O: K7 V8 U7 l% ^( F- af=rand;
x=
       x
) T  X3 w3 @( y2 q* S% c4 X这可怎么办呢？？？急死了，弄了一天了谢谢各位啊！
4 X8 X( f9 l8 R附上函数：负荷贝塔分布的概率密度函数，其中P为待求随机变量，R为已知数
" n! `* e; h1 p- I

qzl19820902

贝塔分布的随机变量的生成，不能用逆函数的方法，可以用gamma函数的办法