Message Boards Message Boards

0
|
3247 Views
|
3 Replies
|
2 Total Likes
View groups...
Share
Share this post:

[?] Calculate the volume of an exponentially distorted ellipsoid?

Consider the following code:

a = 2; b = 1.5; c = 1; tx = 0.175; ty = 0.2; tz = 0.25;
    AbsoluteTiming[
     Volume[ImplicitRegion[
       Exp[tx*x] (x/a)^2 + Exp[ty*y] (y/b)^2 + Exp[tz*z] (z/c)^2 <= 1, {x,
         y, z}], Method -> "NIntegrate"]]

 {39.1684, 
 RegionMeasure[
  ImplicitRegion[
   1/4 E^(0.175 x) x^2 + 0.444444 E^(0.2 y) y^2 + E^(0.25 z) z^2 <= 
    1, {x, y, z}], 3, WorkingPrecision -> MachinePrecision, 
  Method -> "NIntegrate"]}
POSTED BY: Frank Kampas
3 Replies

Alternative way:

 $Version
 (* "11.3.0 for Microsoft Windows (64-bit) (March 7, 2018)" *)

 a = 2; b = 1.5; c = 1; tx = 0.175; ty = 0.2; tz = 0.25;
 R = ImplicitRegion[Exp[tx*x] (x/a)^2 + Exp[ty*y] (y/b)^2 + Exp[tz*z] (z/c)^2 <= 1, {x,y, z}];
 V = DiscretizeRegion[R, Method -> "RegionPlot3D"] // Quiet;
 Volume[V]

(* 6.924361674 *)

or:

a = 2; b = 1.5; c = 1; tx = 0.175; ty = 0.2; tz = 0.25;
AbsoluteTiming[N[Volume[ImplicitRegion[Exp[tx*x]*(x/a)^2 + Exp[ty*y]*(y/b)^2 + Exp[tz*z]*(z/c)^2 <= 1, {x, y, z}]]]]

(*  {202.8073994, 6.924361674} *)
POSTED BY: Mariusz Iwaniuk

Thanks

POSTED BY: Frank Kampas

I had to correct the answer because I had a small mistake.

Regards,MI

POSTED BY: Mariusz Iwaniuk
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract