我正在尝试求解 R 中的积分。但是,当我尝试求解该积分时出现错误。
我试图解决的方程如下:
$$ C_m = \frac{{abs{x}}e^{2x}}{\pi^{1/2}}\int_0^t t^{-3/2}e^{-x^2/t-t}dt $$
我正在使用的代码如下:
a <- seq(from=-10, by=0.5,length=100)
## Create a function to compute integration
Cfun <- function(XX, upper){
integrand <- function(x)x^(-1.5)*exp((-XX^2/x)-x)
integrated <- integrate(integrand, lower=0, upper=upper)$value
(final <- abs(XX)*pi^(-0.5)*exp(2*XX)*integrated) }
b<- sapply(a, Cfun, upper=1)
我得到的错误如下:
Error in integrate(integrand, lower = 0, upper = upper) :
the integral is probably divergent
这是否意味着我无法求解积分?
任何可能解决此问题的方法都将受到高度赞赏。
谢谢。