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

wangjun98989898

如题，已知概率密度函数，拟依据其生成随机变量，想法是先对其积分求累计概率分布，再求反函数，然后由（0，1）均匀分布随机数来生成随机变量，但是，其累计概率密度函数是一个含有hypergeom（超几何分布函数）的函数，怎么办呢？？？
. I, q7 A# B; d: v2 h9 A; J4 `) Q% s* } 具体问题（符合贝塔分布的一个概率密度函数）：
" x3 Q. \( p  A, n" n9 L" lf=int((1/216815.2)*((x/196000)^-0.05)*(1-(x/196000))^-0.05)；MATLAB运行结果是：
: a. u) r8 z. |5 z4 ^& af =

                      2374999999999999993
2 u7 F; N1 x! A: o/ r7 m" m                      -------------------
* U8 Q* G$ V$ D4 F4 _- K4 d                  -5  2500000000000000000
: A+ W" o" H- R* h   0.8928931123 10   x                    hypergeom(
" @9 N# F* w% ^                                                                      -5
         [0.05000000000, 0.9500000000], [1.950000000], 0.5102040816 10   x)
" |6 w  s4 n. z; B8 F当我继续运行
9 P  K' Y' P5 Y: wf=rand;
: `: ?) f, v  b x=
       x
这可怎么办呢？？？急死了，弄了一天了谢谢各位啊！
0 q  L( U# J4 n) g' G附上函数：负荷贝塔分布的概率密度函数，其中P为待求随机变量，R为已知数
% p& o, _* M7 D% l0 g, @

qzl19820902

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