How to avoid Singularity in numirical integration

Question

ABDEL MOUMEN 2019-11-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/493441-how-to-avoid-singularity-in-numirical-integration

评论： ABDEL MOUMEN 2020-3-26

采纳的回答： Raynier Suresh

I want to numerically integrate the following with Matlab:

fun=@(t) (((1 - (-1)^n*cos(n*pi*t/R)))./t);

hn = integral(fun,0,(sqrt(R^2-x^2)));

For now, I just avoided to evaluate the integral in 0, are they any way to encounter the singularity at this point.

Thank you.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Raynier Suresh 2020-3-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/493441-how-to-avoid-singularity-in-numirical-integration#answer_421637

The “quadgk” function can handle singularity if the singularity is present at the boundary. In case if your singularity is not at the boundary you can split the integration domain to place the singularity at the boundary. Refer to the below links for more information,

Quadgk: https://www.mathworks.com/help/matlab/ref/quadgk.html#References

Split integration domain: https://www.mathworks.com/help/matlab/math/singularity-on-interior-of-integration-domain.html