Re: Trouble calculating a volume

• To: mathgroup at smc.vnet.net
• Subject: [mg33840] Re: [mg33806] Trouble calculating a volume
• From: Rob Pratt <rpratt at email.unc.edu>
• Date: Fri, 19 Apr 2002 02:28:56 -0400 (EDT)
• Sender: owner-wri-mathgroup at wolfram.com

```On Tue, 16 Apr 2002, Jay D. Martin wrote:

> Respected Colleagues:
>
> I am attempting to calculate the volume of the frustum of a cone.
> I have been able to calculate the volume of a cone using a triple
> integral as follows:
>
> Integrate[1, {x, 0, h}, {y, (-r)*(x/h), r*(x/h)},
>    {z, -Sqrt[(r*(x/h))^2 - y^2], Sqrt[(r*(x/h))^2 - y^2]}]
>
> This integral does evaluate to the correct answer
>
> (1/3)*h*Pi*r^2
>
> where r is the base diameter and h is the height.
>
> In order to calculate the volume of the frustum of a cone, I have used the
> following equation :
>
> Simplify[Integrate[1, {x, 0, h},
>     {y, -(r1 + ((r2 - r1)*x)/h), r1 + ((r2 - r1)*x)/h},
>     {z, -Sqrt[(r1 + ((r2 - r1)*x)/h)^2 - y^2],
>      Sqrt[(r1 + ((r2 - r1)*x)/h)^2 - y^2]}]]
>
>
> where r1 and r2 are the top and bottom radii of the section and h
> is the height. This integral evaluates to include Logarithms and
> imaginary numbers
>
> (I*Sqrt[h^2*r2]*(r1^2 + r1*r2 + r2^2)*
>     (Log[-2*I*h*r2] - Log[2*I*h*r2]))/(3*Sqrt[r2])
>
> when the solution should be
>
> (1/3)*h*Pi*(r1^2+r1*r2+r2^2)
>
> Which, if you substitute r1=0 and r2=r you get the volume of
> a cone.
>
> I am using Mathematica Version 4.1.2.0
> running on Windows 2000
>
> Thank you for any help
>
> Jay Martin
>
> PS: what I am ultimately after are the moments of inertia of the section
> That is why I am using the triple integral in terms of x, y, and z

Using Simplify with assumptions {h >= 0, r2 > 0} yields your desired

Rob Pratt
Department of Operations Research
The University of North Carolina at Chapel Hill

rpratt at email.unc.edu

http://www.unc.edu/~rpratt/

```

• Prev by Date: MathML Conference: First Announcement
• Next by Date: Re: Passing arguments and pattern matching
• Previous by thread: Re: Trouble calculating a volume
• Next by thread: Redefing a Numberedfigure-Style in a MyReport Style Sheet