How to avoid Singularity in numirical integration

49 views (last 30 days)
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.

Accepted Answer

Raynier Suresh
Raynier Suresh on 24 Mar 2020
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,

More Answers (0)

Products


Release

R2018b

Community Treasure Hunt

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

Start Hunting!