xiefuzuo
e2<-function(a) integrate(function(x) x,lower=0,upper=a)$integral
integrate(e2,lower=0, upper=0.5)$integral
陈放
你编的完全看不懂，还不如你把具体数学式，别人帮你编。
yihui
晕，这么简单个积分需要用integrate函数求么，如果我眼神儿没问题的话，第一行是1/2*a^2，第二行是1/6*.5^3……
关于出错，为啥就是不看帮助，就是不看帮助，不看到底integrate返回的是什么，有没有$integral这个子对象，不要想象。
yihui
原来是重复发帖……就算是重复发帖也把软件说清楚嘛，我不用S-Plus，不知道integrate返回的是什么，上楼说的是R。
多重积分参见?integrate的See Also部分。
xiefuzuo
[attachment=230]
我就是想用软件把g求出来。 R和S-Plus都行。 qnorm是求分位数函数
xiefuzuo
这个积分比较简单， 但是，我想参考这个积分的代码去写一个更复杂点的积分
xiefuzuo
别人给了个对的代码如下，呵呵，R上运行是正确的。奇怪，在S-Plus上有点问题，可能是我的S-Plus盗版的吧
e2 <- function(x) {
out <- 0*x
for(i in 1:length(x))
out <-integrate(function(y) qnorm(y),lower=0,upper=x)$value
out }
integrate(e2,lower=0, upper=a)$value
