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

wangjun98989898

如题，已知概率密度函数，拟依据其生成随机变量，想法是先对其积分求累计概率分布，再求反函数，然后由（0，1）均匀分布随机数来生成随机变量，但是，其累计概率密度函数是一个含有hypergeom（超几何分布函数）的函数，怎么办呢？？？
具体问题（符合贝塔分布的一个概率密度函数）：
$ M+ ]1 v3 @( K2 x. p/ \/ F  kf=int((1/216815.2)*((x/196000)^-0.05)*(1-(x/196000))^-0.05)；MATLAB运行结果是：
f =
1 i5 m2 a- f3 [, Q6 c; \! s
                      2374999999999999993
                      -------------------
                  -5  2500000000000000000
4 Q- c8 p6 ]! L$ E; ]5 {5 V3 p3 n   0.8928931123 10   x                    hypergeom(
4 P; Y" x( Q( f: a  U* u5 N                                                                      -5
- [2 G0 n) v% a' n- L- }         [0.05000000000, 0.9500000000], [1.950000000], 0.5102040816 10   x)
当我继续运行
( K3 J" t7 b/ b. m3 M+ B' l& _7 yf=rand;
/ ^- M' y5 D- S1 U, B  J: D x=
       x
4 i: J5 M/ s+ V5 W, H. c& O7 {这可怎么办呢？？？急死了，弄了一天了谢谢各位啊！
3 S5 h8 z; q, c附上函数：负荷贝塔分布的概率密度函数，其中P为待求随机变量，R为已知数

qzl19820902

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