How to avoid Singularity in numirical integration
显示 更早的评论
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.
采纳的回答
Raynier Suresh
2020-3-24
0 个投票
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,
Split integration domain: https://www.mathworks.com/help/matlab/math/singularity-on-interior-of-integration-domain.html
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
