求大佬看看计算积分误差
mathematica吧
全部回复
仅看楼主
level 5
熵与焓 楼主
NIntegrate[
423091.1143*((2*(1 + 1.46*(Csch[1344750858*x/2])^2)^0.5) - (1.1*
Exp[-1344750858*x])/(1 + 1.46*(Csch[1344750858*x/2])^2) -
2) + 4.518*10^-22*(1.596*10^8*x + 2)/(12*Pi*
x^3*(1 + 5.32*10^7*x)^2) + 1.5*10^10*Exp[-x*10^11/5], {x,
1.6*10^-9, \[Infinity]}]
NIntegrate::ncvb: NIntegrate failed to converge to prescribed accuracy after 9 recursive bisections in x near {x} = {2.01188*10^-8}. NIntegrate obtained 0.00019636995474246995` and 1.3201614457152002`*^-6 for the integral and error estimates.
2020年08月03日 08点08分 1
level 5
熵与焓 楼主
出现这种问题需要咋调整一下?不胜感激各位大佬。
2020年08月03日 08点08分 2
吧务
level 15
不要一看到警告信息就慌。冷静地看一下就会发现误差估计其实并不大。真看它不顺眼的话,解决方法也已经在警告里有提示了:
NIntegrate[423091.1143*((2*(1 + 1.46*(Csch[1344750858*x/2])^2)^0.5) - (1.1*
Exp[-1344750858*x])/(1 + 1.46*(Csch[1344750858*x/2])^2) - 2) +
4.518*10^-22*(1.596*10^8*x + 2)/(12*Pi*x^3*(1 + 5.32*10^7*x)^2) +
1.5*10^10*Exp[-x*10^11/5], {x, 1.6*10^-9, \[Infinity]}, MaxRecursion -> 20]
2020年09月05日 07点09分 3
1