第二个问题 最终的寿命即表示此时它死亡。 所以d(x) 即为你要的函数
比方说 d(20)=0.01 它有1%的概率寿命正好是20岁，即他有1%的概率刚好这个时候死掉。
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)}就行了
