How to avoid Singularity in numirical integration

39 次查看（过去 30 天）
ABDEL MOUMEN
ABDEL MOUMEN 2019-11-27
评论： ABDEL MOUMEN 2020-3-26
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
Raynier Suresh 2020-3-24
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

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by