1
$\begingroup$

I want to visualize in 3d intersecting bodies and compute the volume of their intersection at the same time. It would be very nice if I can rotate 3d view and switch between bodies. If it is too hard or takes long time, it is ok for me to have 3 different 3d view in a document.

My bodies are given by inequalities:

  1. $x^2+y^2+z^2 \ge 10z;$ $x^2+y^2-2z \le 12$.
  2. $x^2+y^2+z^2 -24\le -2z;$ $2 \sqrt{x^2+y^2} \le z+3$.
  3. $x^2+y^2+z^2 +6z\le 16;$ $x^2+y^2+z^2 -6z\le 4$.

Can you please help me?

Update

My current code:

R = 
  ImplicitRegion[
    x^2 + y^2 + z^2 + 6 z <= 16 && x^2 + y^2 + z^2 - 6 z <= 4, {x, y, z}]
RegionBounds[R]
RegionMeasure[R]
RegionBoundary[R]
Show[{R, DiscretizeRegion[ℛ, RegionBounds[R]]}]

Positive things:

$\quad$I can compute the volume of the intersection.

Negative things:

  1. My 3D view show shows something wrong. In my case it should be intersection of two spheres and it shows something close to it, but not what it must be. - SOLVED.
  2. I want to extract the equation of the boundaries to compute multiple integrals, but I do not see how to do it. If possible I want to get cylindrical and spherical boundaries as well.
  3. I do not see the coordinate axes in my 3D view.

Important question on borders:

With my tiny brain I can deduce borders of the region:

Integrate[1, 
  {x, -3, 3}, 
  {y, -Sqrt[9 - x^2], Sqrt[9 - x^2]}, 
  {z, 3 - Sqrt[13 - x^2 - y^2], -3 + Sqrt[25 - x^2 - y^2]}]

How can I extract this information from region $R$? I am interested in Cartesian, cylindrical and spherical coordinates at the same time. $R$ must store at least one these representations internally.

$\endgroup$
5
  • 2
    $\begingroup$ What have you tried so far? Have you converted the formulas to Mathematica syntax? Have you looked at RegionPlot3D, ContourPlot3D, ImplicitRegion, RegionIntersection, Volume? $\endgroup$
    – Szabolcs
    May 21, 2017 at 17:14
  • $\begingroup$ I've updated my answer and posted code, please check it. $\endgroup$
    – Hedgehog
    May 21, 2017 at 17:30
  • $\begingroup$ Ok, it is getting better, can you please check my question about borders of region? Thanks! $\endgroup$
    – Hedgehog
    May 21, 2017 at 17:49
  • 1
    $\begingroup$ For the borders: look up CylindricalDecomposition[]. For the axes: Show[(* stuff *), Axes -> True] $\endgroup$ May 21, 2017 at 18:31
  • $\begingroup$ @J. M. Wow, thanks! One last (but very important thing). CylindricalDecomposition gives me Cartesian coordinates. How can I get cylindrical and spherical coordinates from them with Mathematica? Thanks again! $\endgroup$
    – Hedgehog
    May 21, 2017 at 18:57

1 Answer 1

1
$\begingroup$

I think the intersection is right and what you are seeing is how it should be. see the two spheres below

enter image description here

$\endgroup$
1
  • $\begingroup$ It is correct with BoundaryDiscretizieRegion, but it failed with DiscretizieRegion. Please check my question re borders. $\endgroup$
    – Hedgehog
    May 21, 2017 at 17:46

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

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