我用R计算函数的积分,代码如下:
TT <- c(0.5120,0.5990,0.6860,0.7300,0.8170,0.9040,0.9920,1.0790,1.1660,1.2530)
C1_B <- function(T){integrate(function(tau)f(tau),lower=0.01*T,upper=0.99*T)$value}
f <- function(tau) {integrate(function(tau1)sqrt(1/(tau-tau1)),lower=0.01*tau,upper=0.99*tau)$value}
C <- matrix(rep(0,length(TT)*2),nrow=length(TT),ncol=2)
for (j in 1:length(TT)){
T=TT[j]
C[j,] <- c(T,C1_B(T))
}
但是R报错,显示:
Error in integrate(function(tau1) sqrt(1/(tau - tau1)), lower = 0.01 * :
non-finite function value In addition: Warning message:
In sqrt(1/(tau - tau1)) : NaNs produced
请教大家,这个问题如何解决呢?谢谢了!