q1, dexp是指数分布的概率密度;
r*是生成*分布的随机数,p是*分布的概率分布,d是概率密度,还有q是d的反函数,貌似
R可向量操作,数乘向量是用数乘向量内的每个元素
q2, pi在R里是一个常量,可直接用,3.141593,没有你要求的精度高,你可pi<-3.1415926535...
q3, 不会
q4, integrate 函数能算数值积分;蒙特卡洛积分我不会,甚至都没听过,不是数学专业的,正常正常...
以 ∫ e^(-x) dx 上限: 无限, 下限:0为例啊;
f<-function(x) {exp(-x)}integrate(f, lower=0, upper=Inf)
若是 ∫ e^(-x^3) dx的话,
f<-function(x) {exp(-x^3)}就行了
“面积”的意思是积分吗?还是说x轴之下的部分也算正的面积?x必须是正的吗?
我就当x必须是正的了哈。
只考虑积分的话
F<-function(x) integrate(sin,0,x)$value-1.5uniroot(F,c(0,pi))$root
不然就是
F<-function(x) integrate(function(t) abs(sin(t)),0,x)$value-1.5uniroot(F,c(0,20))$root #20应该够大了
(虽然结果都一样o(╯□╰)o)