Mathematica Stack Exchange is a question and answer site for users of Mathematica. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

In Mathematica how can I compute this integral:$$ \iiint_{D}\sqrt{(1-9z^2)(1-4y^2-9z^2)}\,dx\,dy\,dz$$ where D is the domain:

$$D: x^2 +4y^2+9z^2\le1$$

Please I need help!!!

share|improve this question

closed as off-topic by Sjoerd C. de Vries, Jens, Michael E2, ciao, m_goldberg Jun 4 '14 at 1:34

This question appears to be off-topic. The users who voted to close gave this specific reason:

  • "This question arises due to a simple mistake such as a trivial syntax error, incorrect capitalization, spelling mistake, or other typographical error and is unlikely to help any future visitors, or else it is easily found in the documentation." – Sjoerd C. de Vries, Jens, Michael E2, ciao, m_goldberg
If this question can be reworded to fit the rules in the help center, please edit the question.

I asked how to in fact... – user14738 Jun 3 '14 at 15:26
Where.., there are Integrate Sqrt and all.. :) – Öskå Jun 3 '14 at 15:27
I know...but which is the syntax I've to use? – user14738 Jun 3 '14 at 15:28
Yeah...but no one talks about triple integrals – user14738 Jun 3 '14 at 15:29
Please, check the documentation of Integrate. Specifically, look into the Scope section, subsection Integrals over Regions – Sjoerd C. de Vries Jun 3 '14 at 15:31

I just did this:

Integrate[Boole[x^2 + 4*y^2 + 9*z^2 <= 1]*
 Sqrt[(1 - 9*z^2)*(1 - 4*y^2 - 9*z^2)], {z, -Infinity, Plus[Infinity]},{y, -Infinity, Plus[Infinity]}, {x, -Infinity, Plus[Infinity]}]

MMA quickly returned

share|improve this answer
The order in which you did the integrations is important — doing x then y then z gives the quickest route through the nested integral. Other orderings seem not to evaluate. – Stephen Luttrell Jun 3 '14 at 19:42
Heh. Sometimes it's better to be lucky than to be smart! – Cassini Jun 3 '14 at 19:56

This is perhaps too "creative". Some health checks needed for the series behavior:

region = (x x + 4 y y + 9 z z <= 1);
paregion = Region`ParametricRegion[{{x, y, z}, region}];
k = FullSimplify@
     Normal@Series[Sqrt[(1-9 z^2) (1-9 z^2-4 y^2)], {z,0, #},{y,0, #}] &/@ Range[1, 10, 2];
res = N@Integrate[#, {x, y, z} ∈ paregion] & /@ k

(* {0.698132, 0.488692, 0.480154, 0.477423, 0.476201} *)

So the result is near to 0.476


Mathematica graphics

This is where I've read first about this way for using Integrate[]

share|improve this answer
I don't think this is the most practical way to integrate this. I'm posting it more as a curiosity. – Dr. belisarius Jun 3 '14 at 19:10
+1 for the undocumented integrate use :) – RunnyKine Jun 3 '14 at 19:10
@RunnyKine Thanks for remembering me that! Added a link. – Dr. belisarius Jun 3 '14 at 19:51

Not the answer you're looking for? Browse other questions tagged or ask your own question.