Mathematica - what is wrong with ImplicitRegion[]

  • Context: Mathematica 
  • Thread starter Thread starter PrecPoint
  • Start date Start date
  • Tags Tags
    Mathematica
Click For Summary
SUMMARY

The discussion centers on the issue of integrating the lateral surface area of a cylinder using Mathematica's ImplicitRegion and ParametricRegion functions. The integral $$\int_S z\,dS$$ is computed manually to yield -8π, while Mathematica returns 8π for both regions defined. The user identifies this discrepancy as a bug in Mathematica, confirming that the manual calculation is accurate and the software's output is incorrect.

PREREQUISITES
  • Understanding of surface integrals in multivariable calculus
  • Familiarity with Mathematica version 12.3 or later
  • Knowledge of cylindrical coordinates and parametrization techniques
  • Experience with Mathematica's ImplicitRegion and ParametricRegion functions
NEXT STEPS
  • Investigate the latest updates and bug fixes in Mathematica related to ImplicitRegion and ParametricRegion
  • Explore advanced surface integral techniques in multivariable calculus
  • Learn about alternative methods for surface area calculations in Mathematica
  • Review documentation on Mathematica's integration functions for potential workarounds
USEFUL FOR

Mathematics students, educators, and researchers using Mathematica for surface integrals, as well as developers seeking to understand integration discrepancies in computational software.

PrecPoint
Messages
15
Reaction score
3
TL;DR
Trying to understand how to use an implicit region in Integrate.
I'm working the integral
$$\displaystyle \int_S z\,dS$$

Where S is the lateral (surface) area of the cylinder y^2+z^2=4 cut off by two planes, x=y-3 and x=6-z.

Using the parametrization r(x,\theta)=(x,2\cos(\theta), 2\sin(\theta)) this is pretty straight forward

$$\int 2\sin(\theta)|\frac{\partial r}{\partial x}\times\frac{\partial r}{\partial \theta}|\,dxd\theta=-8\pi$$

However, when I use Mathematica I get the following:

Code:
myReg2 = ImplicitRegion[y^2 + z^2 == 4, {{x, y - 3, 6 - z}, y, z}];
Integrate[z, {x, y, z} \[Element] myReg2]

>8 \[Pi]

Even using my exact parametrization:

Code:
myRegion = 
  ParametricRegion[{x, 2*Cos[\[Theta]], 
    2*Sin[\[Theta]]}, {{x, 2*Cos[\[Theta]] - 3, 
     6 - 2*Sin[\[Theta]]}, {\[Theta], 0, 2 \[Pi]}}];
Integrate[z, {x, y, z} \[Element] myRegion]

>8 \[Pi]

Obviously I'm missing something, but I can't figure it out and it drives me nuts. Any ideas would be welcome!
 
Physics news on Phys.org
Quick update*; so I was not missing something. Apparently, this is a bug in Mathematica.
 
  • Like
Likes   Reactions: jim mcnamara and DrClaude

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
5
Views
2K
  • · Replies 13 ·
Replies
13
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K