q1,dexp是指数分布的概率密度;
r*是生成*分布的随机数,p是*分布的概率分布,d是概率密度,还有q是d的反函数,貌似
R可向量操作,数乘向量是用数乘向量内的每个元素
q2,pi在R里是一个常量,可直接用,3.141593,没有你要求的精度高,你可piq4,integrate 函数能算数值积分;蒙特卡洛积分我不会,甚至都没听过,不是数学专业的,正常正常...
以 ∫ e^(-x) dx 上限:无限,下限:0为例啊;
fintegrate(f, lower=0, upper=Inf)
若是 ∫ e^(-x^3) dx的话,
f 就行了