Mathematica 9 is now available
Student Support Forum
Student Support Forum: 'Bug in Integration' topicStudent Support Forum > General > Archives > "Bug in Integration"

< Previous Comment | Next Comment >Help | Reply To Comment | Reply To Topic
Author Comment/Response
Peter Pein
08/26/12 5:17pm

Hi Amir,

after losing patience I used a simpler approach. y is just a phase shift and has no influence to the value of the integral. 2 Pi is a multiple of the period w.r.t. x and due to the even powers there is a bit more symmetry.

So defining
f[x_, y_, a_] = (1 + a^4 Tan[2 x - 2 y]^2)/(1 + a^2 Tan[2 x - 2 y]^2);

the value at x==k*Pi/4, k integer, y==0 has to be determined by the use of limits. So I added an intermediate point of integration at Pi/4.

the value of the integral is

value = Simplify[
4*Integrate[f[x, 0, a], {x, 0, Pi/4, Pi/2},
Assumptions -> Element[a, Reals]]

this evaluates to
Pi*(2 - a + 2*a^2 - Abs[a])
in a few seconds.

For a in {-1, 0, 1} this value is 2*Pi.

maybe the wrong result you got is due to the removable singularities at x = k*Pi/4 + Mod[y, Pi/4]?


URL: ,

Subject (listing for 'Bug in Integration')
Author Date Posted
Bug in Integration Amir 08/21/12 09:32am
Re: Bug in Integration Nasser M. Ab... 08/21/12 11:41am
Re: Bug in Integration jf 08/21/12 1:13pm
Re: Re: Bug in Integration Amir 08/24/12 10:31am
Re: Bug in Integration Peter Pein 08/26/12 5:17pm
Re: Re: Bug in Integration Amir 08/31/12 09:41am
< Previous Comment | Next Comment >Help | Reply To Comment | Reply To Topic