利用概率密度函数生成随机变量过程中出现超几何分布函数问题
 

[IMG]http://emuch.net/bbs/viewthread.php?tid=4339838&pid=2&page=1#pid2[/IMG]
如题，已知概率密度函数，拟依据其生成随机变量，想法是先对其积分求累计概率分布，然后由（0，1）均匀分布随机数来生成随机变量，但是，其累计概率密度函数是一个含有hypergeom（超几何分布函数）的函数，怎么办呢？？？
具体问题（符合贝塔分布的一个概率密度函数）：
f=int((1/216815.2)*((x/196000)^-0.05)*(1-(x/196000))^-0.05)；MATLAB运行结果是：
f =

2374999999999999993
-------------------
-5 2500000000000000000
0.8928931123 10 x hypergeom(
-5
[0.05000000000, 0.9500000000], [1.950000000], 0.5102040816 10 x)
当我继续运行
f=rand;
x=
x
这可怎么办呢？？？急死了，弄了一天了谢谢各位啊！
附上函数：负荷贝塔分布的概率密度函数，其中P为待求随机变量，R为已知数