Warning: Explicit integral could not be found.

Hi
I'm trying to solve this integral but it gives me this error: Warning: Explicit integral could not be found.
syms a b
Q = (1/(2*pi*lams)) * exp( - (a^2 + b^2)/(2*lams));
q = int(int(Q,b,0,sqrt((lams^2) - (a^2))),a,0,lams);

 Réponse acceptée

Roger Stafford
Roger Stafford le 26 Juin 2013
Apparently 'int' is not quite smart enough to make a change of variables. If you change the variables to polar coordinates:
a = r*cos(theta)
b = r*sin(theta)
the double integral is easily found and has the value (if my calculus isn't too rusty)
q = 1/4*(1-exp(-lams/2))
(The area covered by the integration in the a,b coordinates is a quarter circle.)

Plus de réponses (1)

Walter Roberson
Walter Roberson le 26 Juin 2013

0 votes

It does not have any obvious closed form integral. Do you have reason to expect that there is a closed form for it?

3 commentaires

Roger Stafford
Roger Stafford le 26 Juin 2013
But it does, Walter, as I have just pointed out, but 'int' is not able to find it.
Walter Roberson
Walter Roberson le 26 Juin 2013
My answer was being prepared before yours was posted ;-)
Roger Stafford
Roger Stafford le 26 Juin 2013
I'm sorry. I should have made allowance for that before complaining.

Connectez-vous pour commenter.

Community Treasure Hunt

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

Start Hunting!

Translated by