Я пытаюсь интегрировать функцию (от 0 до бесконечности), которая, как я знаю, должна сходиться, потому что она имеет форму, подобную t, умножающую отрицательную экспоненту, но я получаю следующую ошибку:
Ошибка в интегрировать (c11t, нижний = 0, верхний = Inf, подразделы = 2000): значение конечной функции
Этот код показывает ошибку:
c11<-function(t){
if(Tr^2>4*D){
exp(-Tr*abs(t)/2)*sigma^2*(R*cosh(R*abs(t)/2)*(D+gamma22^2)+Tr*(-D+gamma22^2)*sinh(R*abs(t)/2))/(2*D*R*Tr)
}else if(Tr^2<4*D){
exp(-Tr*abs(t)/2)*sigma^2*(R*cos(R*abs(t)/2)*(D+gamma22^2)+Tr*(-D+gamma22^2)*sin(R*abs(t)/2))/(2*D*R*Tr)
}else{
exp(-Tr*abs(t)/2)*sigma^2*(Tr^2*(2-Tr*abs(t))+4*(2+Tr*abs(t))*gamma22^2)/(4*Tr^3)
}
}
c11abs<-function(t){
abs(c11(t))
}
c11t<-function(t){
t*c11abs(t)
}
gamma22<--10
Tr<-1.020408
D<-0.2040816
sigma<-0.2
R<-0.4742429
t11<-integrate(c11t,lower=0,upper=Inf,subdivisions = 2000)$value/integrate(c11abs,lower=0,upper=Inf,subdivisions = 2000)$value
Это выглядело очень странно мне, и я попытался построить функцию c11t, получая кривую, которая начинается с нуля при t = 0, немного растет, а затем быстро уменьшается экспоненциально, вероятно, до нуля, в то время как c11 затухает как экспонента (и всегда больше, чем ноль) предполагая, что интеграция должна быть возможной. Может кто-нибудь, пожалуйста, помогите мне с этой проблемой? Заранее большое спасибо