Message Boards Message Boards

7 Replies
2 Total Likes
View groups...
Share this post:

Double Integration runs a long time

Posted 10 years ago


I am trying to solve following expression

int(int((1-exp(-5.5/cos(x)))*sin(x), x = 0 .. arctan(300*cos(y)+sqrt(12.25-90000*sin(y)^2))), y = 0 .. Pi)

But I am not getting results.

Please help me!!!


7 Replies

An additional reply that was here should have been a new discussion. It was moved to

POSTED BY: Moderation Team

Thanks Bruce. I'll make sure to do this in future posts.

POSTED BY: Michael McCain

I think this is more a mathematical question, but anyway...

In this case one is lucky, because the inner integral can be done analytically! What remains is then just a single integral which can be calculated numerically at least. Here is my result:

In[1]:= ClearAll["Global`*"];

In[2]:= fInt[x_] = Integrate[(1 - Exp[-5.5/Cos[x]]) Sin[x], x]

Out[2]= -1. Cos[x] + E^(-5.5 Sec[x]) Cos[x] + 
 5.5 ExpIntegralEi[-5.5 Sec[x]]

In[3]:= NIntegrate[
 fInt[ArcTan[300 Cos[y] + Sqrt[12.25 - 90000 Sin[y]^2]]] - 
  fInt[0], {y, 0, Pi}]

Out[3]= 3.13316 + 0. I

In[4]:= Chop[%]

Out[4]= 3.13316

Cheers Henrik

POSTED BY: Henrik Schachner

Thank you for including code and not just a picture of code.

POSTED BY: Bruce Miller

I tried what you tried... and did not get a result. The kernel just keeps running without evaluating. I then tried a numeric integration and got a result. I could be wrong here. Anyone else have any ideas?

POSTED BY: Michael McCain

Same thing is happening with me. But I have solved analytically the inner integral.Now the problem is with second integral, i.e. integration with respect to y, from 0 to pi. I got following result by first integralenter image description here


I apologize. I missed a parenthesis in my code. I fixed it. I am getting the same result as Henrik. Here is it...

I believe Henrik is correct. The inner integral is solvable analytically. But the outer integral might not be (or it's too large for Mathematica to calculate), hence using numerical integration on the outer integral. I'm leave it running for an hour or so just in case and get back to you if Mathematica is able to do the outer integral analytically.

POSTED BY: Michael McCain

Group Abstract Group Abstract