发散函数的解析或数值积分

计算科学 收敛 积分方程
2021-12-04 15:21:21

我正在尝试以数字方式执行以下积分

rinroutΣ(r)rr(1rrE(L)+1r+rK(L))dr,

在哪里

L=4rr(r+r)2,

是第一类和第二类完全椭圆积分都是是一些积极的,表现良好函数(没有奇点或不连续性,...)。KErrinrout>0rinrroutΣ(r)

最后,我需要为给定的对结果进行数值评估,但是由于两个奇点,我常用的 python 积分器或无法收敛或给出不可靠的结果:Σ(r)scipy.integrate.quadscipy.integrate.romberg

  • 1/(rr)第一个加法中的
  • K(L) as在第二个加法(rrL1

任何人都可以帮助我以数字方式评估它或重写它以使通常的方法不会失败吗?

0个回答
没有发现任何回复~