求解扔骰子概率问题，R语言

Question

用R语言模拟扔3个骰子，扔出后3个骰子的相乘的值小于90的概率，至少模拟10万次求概率，保留6位数字。

还有一道题，一种灯使用寿命是指数分布，平均寿命3000小时，求这一个这种灯使用寿命在2000-3000小时之间的概率。

在线求解，求高手解答。最好把步骤写下来。

Answer 1

第一个取名为dice

dice<-function(N){ #建一个名字叫dice的方程
count<-0 #计数器：初始化为0
for(n in 1:N){ #运行N次
result<-1 #这个是色子点数的乘积，初始化为1
for(i in 1:3){ #投三次
result<-result*sample(1:6,1) #sample(1:6,1)的意思是从1到6中随机选一个
}
if(result<90) #乘积小于90的话
count<-count+1 #就自增1
}
prob<-round(count/N,6) #概率保留6位小数
return(prob) #函数返回这个概率
}

编好了在控制台source一下你这个R文件，然后输入dice(100000)
> dice(100000)
[1] 0.84911
> dice(100000)
[1] 0.85275
我跑了2次

第二题都不用编函数了，直接算吧。平均寿命3000小时，那么指数分布的参数lambda就是1/3000,如果X~Exp(1/3000)，题目求P(2000<X<3000)的值
R里有函数叫pexp(q,rate=r)是指数分布的分布函数P(X<q)，其中 q是寿命，r是lambda
所以P(2000<X<3000)等于
>pexp(3000,rate=1/3000)- pexp(2000,rate=1/3000)
[1] 0.1455377