Integration of some function containing a singularity

11 次查看（过去 30 天）
Xianjie
Xianjie 2012-3-20
Dear All,
These days I am working on some integration integrate(cos(x)*sin(x)/x^t),x(0,a),a is a value. the value fot t is 1, 2, 3. The function also have the form cosx*cosx and sinx*sinx. The value x = 0 is the singularity. How to overcome this singularity. The function quadgk can integrate some funtion with endpoint singularity. x = 0 is the end point singularity. But the integral is not correct when use quadgk. How to solve this problem.
Thank you for your concern.

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

采纳的回答

Mike Hosea
Mike Hosea 2012-3-20
QUADGK has no trouble when t=1, but for t=2,3,4,... the correct answer, if I am not mistaken, is Inf. You were expecting QUADGK to return Inf for you?
  3 个评论
显示 1更早的评论
Xianjie
Xianjie 2012-3-21
But when the x=0 is not selected, choose x=0.001, the integration is ok. How about this phenomenon. Is there some numerical evaluate the results? I calculate the strain energy of circular plate. x=0 is the center of circular plate.
Thanks
Mike Hosea
Mike Hosea 2012-3-21
I believe QUADGK handles that case properly, though you may well run into precision issues.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 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