If C is the temperature in degrees Celsius and F is the temperature in degrees Fahrenheit,
then the following exact relation holds:
F + 40 = 1.8 ( C + 40 )
This allows an immediate answer to the trivia question:
"What temperature is the same in both scales, Fahrenheit and Celsius?"
(Answer: "40 below").
More importantly, it's useful to remember the conversion in the form given above,
because the numbers involved are easy to manipulate mentally:_{ }
To multiply by 1.8 (from C to F), you multiply by 2 (easy)
and subtract 10% from the result (almost as easy and perfectly accurate)._{ }
To divide by 1.8 (from F to C), you divide by 2 and add 10% to the result.
(11% is better, since the exact percentage is 11.11111...%)
For example, to convert 20°C you double 20+40 which gives you 120,
subtract 12 (that's 10%) to obtain 108,
and subtract 40 to obtain 68°F.
The result is exact and not much more difficult to obtain
than the dubious approximations given by overly "simplified" formulas...
´
1
.
8
=
¸
1
.
8
=
With a simple calculator, exact conversions are performed either way
in only 5 keystrokes, since it's easy to account for the simple 40° translation
(the same in both scales) as you enter or read the data.
°C
40
35
30
25
20
15
10
5
0
5
10
15
20
25
30
35
37
40
100
°F
40
31
22
13
4
5
14
23
32
41
50
59
68
77
86
95
98.6
104
212
Alternately, you may memorize [part of] the above table and notice that
an increase of 5°C corresponds exactly to an increase of 9°F...
(S. P. of Piscataway, NJ. 20000715)
Can you determine the speed of a car by knowing the rpm of the motor,
the gear [ratio] and tire [diameter]? If so, what is the formula? [...]
If D is the diameter of the tires, the car moves forward a distance
pD with each revolution of the tires.
If there are F such revolutions per unit of time, it moves forward at a speed of
pD´F per unit of time.
If D is in inches and F is in revolutions per minute (rpm),
the above result will thus be in inches per minute!
Now, there are 63360 inches in a mile and 60 minutes in an hour,
so one inch per minute equals 60/63360 miles per hour or 1/1056 mph.
Therefore, we have:
V(mph) = D(in) ´ F(rpm) ´
(p/1056)
In the automotive industry, the coefficient 1056/p
is usually taken to be equal to 336
(it's actually 336.1352398...) and the formula corresponding to the above is memorized as:
(mph)(gear ratio)(336) = (tire diameter in inches)(rpm)
With the above numerical approximation, both formulas are identical if we consider that
(rpm)/(gear ratio) is the rate of rotation of the tires.
This means that the industry defines the
gear ratio as the number of rotations
of the driveshaft for each turn of the tire.
This ratio is normally more than 1, except in overdrive gear.
So, calling R this "gear ratio", your final answer is:
V(mph) = D(in) ´ F(rpm)
´ (p/1056) / R
or, approximately:
V = (D ´ F) / (336 ´ R)
where V is the velocity (in mph), D the tire diameter (in inches), F the engine rpm,
and R the gear ratio (the number of engine revolutions for each turn of the tires).
(C. W.. of Grandy, NC. 20001023)
If a car is doing 0  60 mph in 4.59 seconds, how far did the car travel?
Well, the most common way to answer this question is to assume that the car
has a constant acceleration, in which case the distance it travels starting at
zero speed is the same distance it would have traveled at a constant velocity
equal to half the final velocity.
Here that's 30 mph times 4.59 s or, with the proper conversion factors:
Note the foolproof way of converting units, by introducing unity factors
(like 5280 ft/mi) to cancel those units you don't want in the result.
Now, however, any engineer will tell you that acceleration is not constant,
so the above answer is wrong in practice.
It just gives you a rough idea...
In theory, a rocket vehicle could have a very large acceleration at the very beginning and
almost nothing after a small fraction of a second.
When that's the case the distance would be (almost) twice the above distance.
The opposite (silly) case is when you push by hand the above rocket vehicle for more than
4 seconds (moving it only a few feet, a few inches, or nothing at all),
then the rocket fires up and reaches 60 mph in a fraction of a second,
having traveled only a short distance in the process...
(Carlos of Jersey City, NJ. 20001115)
[...] If my car can do a 1/4 mile in 15.9 seconds at 94 mph,
how fast is it going at 0.1 or 1/10 of a mile?
It accelerates to 60 mph in 7.5 seconds.
What you are asking is to give a point on a curve (speed as a function of time)
about which we have preciously few properties.
(A total of 3 points including the trivial 0s,0mph
and the overall integral between 0 and 15.9 s, which is 1/4 mile.)
This is not enough to do a totally accurate job,
but we can make a fairly precise educated guess:
First, let's dispose of the issue of units.
If we measure times in seconds and speeds in mph,
the consistent unit of distance is then the mphsecond (mphs),
which is simply 1/3600 of a mile.
A quartermile is 900 of these units, a tenth of a mile is 360, and the distance from the
1/10 mile point to the 1/4 mile point is 540.
In particular, since a quartermile is 900 mphs,
the entire integral of our speed vs. time curve is 900.
For example, let's try a curve of the type V = C(1exp(t / T)),
where V is the velocity and t is the ellapsed time (whereas C and T are some constants).
With such a speed curve, the distance traveled in time t to reach speed V is
CtVT = 900 for the quarter mile.
The 1/4 mile in 15.9 s at 94 mph imposes [trust me!]
the following parameters on such a curve (forgive the ludicrous precision):
T=12.639196110577s and c=131.32606505623 mph.
With this formula, you do a distance of 360
(that's 1/10 of a mile with our "weird" units)
in 9.3475 seconds at 68.64 mph.
That's our first approximative answer to your question.
Note, however, that this approximation would mean a speed of only 58.775 mph
is reached after 7.5 s or, alternately,
that 060 mph is achieved in 7.715 seconds.
Our crude approximation, is thus not too far from what your car is actually doing,
but you may not like the 0.2s and/or 1.2 mph discrepancy.
I don't think I am far off in guessing that your car would do the 1/10 of a mile
in 9.3 s at 69 mph but I could definitely use some raw data
(acceleration curves of other cars, say) to finetune the above.
I hope somebody will be able to provide
this type of information.
A more elementary approach is to consider that,
from the time t=7.5 s to the time t=15.9 s,
the speed increases roughly linearly from 60 mph to 94 mph,
at a steady rate of
(9460) / (15.97.5) = 85/21, or about 4.04762 mph/s.
With this assumption, the speed V reached at a time t between 7.5 s and 15.9 s
is simply V=60+(85/21)(t7.5).
If t is the time when the 1/10 mile point is reached,
the distance of 1/41/10 mile (or rather 540 mphs with the "proper" unit
introduced above) is simply the area of a trapezoid, and we have:
(15.9t)(V+94)/2 = 540.
All told, this means that t is a solution of the equation
540 = (15.9t) [ 60+(85/21)(t7.5) + 94 ] / 2
The nonnegative solution of this quadratic equation is about 9.18431 s,
corresponding to a speed V of about 66.8174485 mph.
In other words:
Your car does the 1/10 mile in 9.2 s at 67 mph.
Overall, this more elementary approach may very well give a better, more reliable
estimate, given the data at hand.
Hsuan Liu (20020506; email)
With 200 hp applied to a 2350 lb vehicle, what's the 0 to 60 mph time?
Answer: at least2.57 seconds (see formula below), but this silly lower limit
would be the correct answer only if 200 hp was the actual
average mechanical power applied to the tires during acceleration
[assuming no skidding and deducting the energy lost to friction].
Such an acceleration is next to impossible with
current car technology:
The current record 060 mph time is 3.07 seconds; it was set by
a Ford RS200 Evolution (a 1986 "Group B" rally car)
at the Millbrook Proving Ground (UK), in May 1994.
Therefore, what you have in mind is most probably not this kind of theoretical
net power, but rather the rating most often
listed by manufacturers,
which is the maximum engine output.
The fastest production car is currently the McLaren F1 pictured at right,
whose 6.1 L engine is rated at 627 hp...
Here are a few actual examples which tell you what percentage [last column]
of the rated maximum power is the actual average net power propelling the car
over the 0 to 60 mph acceleration test whose duration is listed:
The last four lines of the above table give the answer to your question,
assuming the quoted power of 200 hp is one of the following:
(1) The rated maximum
power, if the average is 40% of that.
(2) The rated maximum
power, if we use the rule of thumbbelow.
(3) The rated maximum
power, if the average is 50% of that.
(4) The actual average power applied.
[Not a realistic assumption.]
Now, what's the formula and where does it come from?
Well, if M is the mass of the car, and V is the speed reached from a standing start
after applying an average power P for a time T, we have:
P T = ½ M V^{2}
The work done over the time of acceleration (PT) is equal to the vehicle's kinetic
energy (½MV^{2}, if we neglect the rotational energy of the wheels).
The above formula is valid as is if we use any consistent set of
units (for example, P in watts, T in seconds, M in kilograms,
and V in meters per second; 60 mph = 26.8224 m/s).
If, however, you insist on having P in horsepowers
(1 hp is exactly 745.699871582 W, or just about 745.7 W)
and M in pounds (1 lb is exactly 0.45359237 kg) at a speed of
60 mph (26.8224 m/s), the above formula becomes:
So, if (P/hp) is 200 and (M/lb) is 2350, the 060 time in seconds (T/s)
will be such that 200 (T/s) / 2350 = 0.2188....
In other words, T is 2.571016... s, as unrealistically advertised.
Divide this time by 0.4 to obtain T = 6.42754... s,
if the value of 200 hp reflects the maximum engine power and you estimate the
actual average net power to be only 0.4 = 40% of that.
Rule of Thumb:
If you live in the poundsandhorsepowers world and want to memorize a simple
ruleofthumb formula, just assume that the average power is about 43.76 % of the
rated power, so everything boils down to the following approximation:
(060 time in seconds) »
½ (test weight in lb) / (rated engine power in hp)
Recall that: test weight
=
curb weight + weight of the driver (160 lb)
Anonymous query via Google (20041123)
Power to Speed Ratio Horsepower to thrust conversion...
[at a given speed along the thrust]
The power P of a force (thrust) F applied to an object moving at
velocity V is obtained by dotting
F into V :
P = F . V
When thrust and velocity do not have the same direction, the two quantities on the
righthand side must be considered to be vectors (which is why they're shown in bold)
and the product is understood to be a scalar product, or "dot product".
We call speed V the magnitude of the velocity vector V, whereas
the "speed along a certain direction" (v) is the projection of the velocity vector
on that direction.
In what follows, "v" is the speed along the direction of the thrust,
and q is the angle between thrust and velocity:
P = F V = F v
= F V cos(q)
The above equation is directly applicable if a consistent system of units
is used:
The proper SI units for power, thrust and speed are, respectively, the watt (W),
the newton (N) and the meter per second (m/s).
With any other mix of units, an extra numerical factor appears, which can be obtained
by introducing standard ratios equal to unity
(like 5280 ft / mile).
For example, using the horsepower (hp),
the poundforce (lbf) and the "mile per hour" (mph) as units, we have:
(P/[hp]) [hp]/[W]
=
(F/[lbf]) [lbf]/[N] (v/[mph]) [mph]/[m/s]
(P/[hp]) [550 ft lbf]/[J]
=
(F/[lbf]) [lbf]/[N] (v/[mph]) [5280 ft / 3600 m]
This boils down to:
375 (P/[hp])
=
(F/[lbf]) (v/[mph])
The thrust (in lbf) is thus 375 times the
power to speed ratio (in hp per mph).
The thrust is zero for a car at cruising speed
(or at top speed, which is cruising speed at full throttle):
The engine's power is then entirely used to overcome friction and air resistance,
leaving no net power (P) to accelerate the vehicle.
Ahsin Sattar (20040629; email)
How does the volume [in cc] of an engine relate to its power [in hp]?
With a given technology,
an engine has to have a certain size to produce a given power.
Usually, only usable forms of power are considered:
mechanical or electrical power, but not power wasted as heat.
In an internal combustion engine, the "size" of each cylinder is universally
understood to be the maximum volume of its combustion chamber
(the total size of the engine is, of course,
the sum of the sizes of its cylinders).
This is most commonly given in cubic centimeters (cc),
liters (1 L = 1000 cc) or
cubic inches (1 cu in = 16.387064 cc).
In automotive parlance, the bore
is the crosssectional area of a piston and the stroke
is the distance it is allowed to travel. Multiplying
bore by stroke gives the swept volume which is the difference
between the maximum and the minimum volume of the combustion chamber;
the ratio of those two quantities is the compression ratio.
The power of an engine is best given in watts (W) but manufacturers most often
rate it in horsepowers, of which there are (unfortunately)
two flavors; the Imperial horsepower (hp)
is about 745.7 W whereas the "metric" horsepower (ch) is about 735.5 W.
The manufacturer's power rating corresponds to the maximum power output
at some optimal regime.
This regime is often stated using "rpm" as the frequency unit
("revolution per minute" = 1/60 Hz).
The powertovolume ratio depends very much of the engineering and technology involved.
For example, the French legal upper limit for a "125 cc" motorcycle is 15 ch,
which corresponds to a ratio of 0.12 ch/cc, or precisely 88.25985 W/cc.
Typically, 125 cc "street" motorcycles are rated around 70 W/cc,
but this can be as low as 50 W/cc.
On the other hand, highperformance 125 cc competition engines have been rated
as high as 34 ch (200 W/cc)...
Mark Nance (20050510; email)
What's the optimal gear ratio to maximize top speed on a flat road?
At uniform speed, no power is used to accelerate the car.
Instead, the engine's mechanical output (whatever is not directly wasted as heat)
is used to overcome mechanical resistance from the following sources:
Powertrain between the engine and the gearbox.
Powertrain between the gearbox and the wheels.
Periodic deformation of the tires.
Contact between the tires and the road.
Air resistance (assuming no wind with respect to the road).
Without lubrication [!] the
dry frictional forces for moving parts would
be roughly constant and the power lost to these would be proportional to how fast they
rub against each other, or against their bearings.
With good lubrication, such "raw" friction is greatly reduced
and what dominates (except at low speeds) are socalled viscous
forces of magnitude proportional to the speeds involved, entailing a smaller power loss
proportional to the square of that.
The relevant speeds in the powertrain are proportional either
to the engine's rpm (between engine and gearbox) or
to the car's speed (between gearbox and wheels).
For simplicity, we'll just assume that the power lost
inside the gearbox is approximately the sum of two components:
One of each of the two types presented above for the rest of the powertrain.
The deformation of the tires produces heat at a rate proportional to its
frequency: This power loss is thus proportional to the speed of the car
(similar to what would be unlubricated friction in the powertrain after the gearbox).
A rough idea of this effect's magnitude may be obtained by comparing
fuel efficiencies under controlled conditions, when the tires are properly
inflated and when they're not.
Air resistance is a complicated thing but it has two main components:
The viscous one is only significant at very low speeds.
The other component is essentially
a quadratic effect, involving a drag force
proportional to the square of the car's speed
(its power is thus proportional to the cube of the speed).
All told, let V be the top speed of the car for a gear ratio R, and P(x)
the power output [full throttle] of the engine at x rpm
(x = kRV for a constant k which is inversely proportional to the
diameter of the tires, as discussed above). We have
P( kRV ) =
a R^{2 }V^{2} + b RV + c V^{ 3} + d V^{ 2} + e V
In this equation, the constant coefficients k,a,b,c,d,e are characteristics
of the vehicle (a, b and c don't depend on tire size).
For a given R, the vehicle's top speed is obtained by solving this equation.
At the optimal gear ratio R for which the top speed V is greatest, a minute
increment dR wouldn't change V (i.e., dV = 0)
so the derivatives of both sides with respect to R are equal.
Dividing both sides of the resulting equation by V, we obtain:
k P'( kRV ) = 2 a R V + b
Solving the above two equations simultaneously gives the optimal gear
ratio R, for which the greatest top speed V can be obtained.
This requires detailed knowledge of your engine (the function P)
and of both parts of the powertrain.
Not to mention the vehicle's aerodynamics
(the drag coefficient we called "c").
You may acquire the relevant knowledge with 5 actual measurements of the
top speed of the vehicle, under engine conditions that you can reproduce
on the bench (for example, full throttle with 5 different
gears).
Reproducing the rpm and fuel intake on the bench allows you to measure
the actual engine output in these 5 cases.
This leaves you with 5 instances of our first equation,
which form a linear system
that yields easily the values of its 5 unknowns a,b,c,d,e
(the coefficient k being known from other considerations).
Once this is done, you still have to solve the two
nonlinear simultaneous equations mentioned above...
This requires only a plot of the P function near the expected solution,
from bench measurements or manufacturer's data.
You are awesome! Thank you for [posting] that information for
myself and everyone.
I have donated to your website to help keep it going.
Great information. Thank you so very much!
Joshua P. Gatcomb (20020715; email)
When I was in middle school, I found a way to derive the
Law of Cosines from
Heron's formula
[details attached, in 20 steps]...
Have you seen it done this way before?
Well, yes (in the backward sort of way described below).
Congratulations for rediscovering by yourself this connection between the two formulas, though.
Heron's formula states that
A^{2} = s (sa)(sb)(sc),
where A is the surface area of a triangle of sides a, b, c, and
of semiperimeter
s = ½ (a+b+c).
One quick way to derive Hero's formula is to start with the most commonly known expression
for the area A, namely: A = ½ ah,
where h = b sin g
is the height of the triangle (of base a)
and g is the angle between a and b.
We have:
A^{2}
=
(½ ab)^{2} sin^{2} g
= (½ ab)^{2} 
(½ ab cos g)^{2}
=
(½ ab)^{2}

(¼ [ a^{2} + b^{2} 
c^{2} ] )^{2}
[ Law of Cosines ]
[¼ c^{2} 
¼ (a  b)^{2 }]
[¼ (a + b)^{2} 
¼ c^{2 }]
=
[ ½(c  a + b)
½(c + a b) ]
[ ½(a + b  c)
½(a + b + c) ]
=
[ (s  a)
(s  b) ]
[ (s c)
s ]
The whole thing may be used backwards:
If you assume Heron's formula, the last equality holds and the above elementary
manipulations thus prove the equality on the second
line without invoking prior knowledge of the Law of Cosines.
So, Heron's formula does imply the Law of Cosines,
with a sign ambiguity
which must be lifted by other means (since we've only equated the squares of two
quantities here, they could be either equal or opposite).
However, this would really be a backward, because the
Law of Cosines is considered far
more fundamental than Heron's Formula (a.k.a. Hero's Formula).
Surface Area of a Spherical Triangle :
In spherical geometry,
the sides of a triangle are arcs
of great circles (the radius of a great circle is equal to the radius
of the sphere). Normally, the radius of the sphere is used as a
natural unit of length (so, the length of every side of a
spherical triangle is between
0 and 2p ).
There is a beautiful counterpart of Hero's formula which gives the
surface area A =
e R^{2}
of such a spherical triangle.
It is due to the Swiss mathematician
Simon Lhuilier (17501840):
e is (also) equal to the "spherical excess"
of that spherical triangle, namely the angle by which the sum of its three angles
exceeds p.
The locution "spherical excess"
(French: excès sphérique)
was coined around 1626 by the Frenchborn Dutch mathematician
Albert
Girard (15951632) who, incidentally, contributed to
the popularity of Fibonacci numbers
(0,1,1,2,3,5,8,13,21,34,55,89...) by defining them inductively:
Girard is also remembered for discovering
the GirardWaring expansions
more than 130 years before Waring.
(20020716) Brahmagupta's Formula
How does Heron's formula generalize to quadrilaterals?
Brahmagupta's formula states that
A^{2} = (sa)(sb)(sc)(sd),
where A is the surface area of a cyclic quadrilateral
(that's to say a quadrilateral inscribed in a circle)
of sides a, b, c, d, and of
semiperimeter
s = ½ (a+b+c+d).
That formula was devised around AD 620 by the Indian mathematician
Brahmagupta (598670),
as a generalization of Heron's formula (d = 0).
Brahmagupta didn't give a proof of his formula and
never mentioned that it was only valid in the special case of quadrilaterals
that can be inscribed in a circle.
According to Pierre Dedron and Jean Itard
(Mathématiques et mathématiciens, p. 126, Magnard 1972)
alBiruni
(9731048) was the first mathematician who properly insisted
that the formula was only good for cyclic quadrilaterals.
The formula can be generalized to any quadrilateral by introducing
an angle q equal to the halfsum of either pair
of opposite angles:
(20031208; email, name withheld ) _{ } Bretschneider's Formula
Does the above relation give the area (A) of a quadrilateral in terms of
its sides (a,b,c,d)
and diagonals (p,q) without any restrictions?
The short answer is: Yes.
This formula was established in 1842,
by the German mathematician Carl Anton Bretschneider (18081878),
an author who earned a living as a highschool teacher in
the town of Gotha (Thüringen).
Bretschneider's Formula does hold for convex
quadrilaterals and concave chevrons, and it's also good for [the
signed areas of] crossed quadrilaterals.
A "crossed quadrilateral" is a butterflyshaped figure
which appears as two triangles sharing a single point
(this point where two edges cross is not counted as a vertex).
Travelling along the 4 edges of such a quadrilateral,
one encircles a triangle clockwise and the other triangle counterclockwise.
The area of a crossed quadrilateral is
thus usually defined as the difference
between the areas of those two triangles.
Bretschneider's Formula does give [the absolute value of] this quantity.
Proof:
Let's call a and b
the angles from a diagonal (p) to each of the two adjacent sides a and b
(for convex quadrilaterals, these two angles have opposite signs,
but this need not be so for other quadrilaterals).
The area of the quadrilateral is either
the sum or the difference of the two triangles
of base p. More precisely:
4 p^{2}q^{2}

(a^{ 2}  b^{ 2} +
c^{ 2}  d^{ 2 })^{ 2} =
4 p^{2}q^{2}

( 2 b p cos b

2 a p cos a )^{ 2} =
4 p^{2} [
a^{ 2}
+ b^{ 2}

2 a b cos (ab)

( b cos b

a cos a )^{ 2 } ]
=
4 p^{2} [
a^{ 2} sin^{2 }a
+
b^{ 2} sin^{2 }b

2 a b sin a sin b ]
=
( 2 p [ a sin a

b sin b ] )^{ 2 }
=
(4A)^{ 2}
Bretschneider's Formula
is just the simplest way of expressing the area of a quadrilateral in terms of
[ the squares of ]
its sides and diagonals. There are infinitely many other ways to do so,
since those 6 quantities are linked by the following polynomial relation.
( p is the diagonal joining the corner of a
and b to the corner of c and d,
while q is the other diagonal.)
(That's a quadratic equation in
the square of every length involved.)
(20070920) Vectorial area of a triangle...
or a quadrilateral :
The area of a quadrilateral is
half the crossproduct of its diagonals.
The area of a triangle is half the length of the crossproduct
of two of its sides.
A straight consequence of this wellknown fact is that the
area of a (planar) quadrilateral is half the length of the
crossproduct of its diagonals
(itself equal to twice the
vectorial area).
Proof :
A planar quadrilateral ABCD consists of two triangles ABC
and ACD whose
area vectors are collinear (so that the sum of
the lengths is the length of the sum).
Twice the area vector of ABCD is thus:
2 S =
AB´AC + AC´AD
=
AB´AC 
AD´AC
=
DB ´ AC
Kristiana KandereGrzybowska
(20040517)
_{ } Parabolic Segments
[In my postdoctoral work in cell biology, I need to estimate]
the area of a modified triangle with one convex side and two concave sides.
[...]
The simplest approach is to consider that such sides are
curved roughly in the shape of parabolas (possibly asymmetrically with
respect to the straight sides).
The surface between a parabolic arc and its chord is
called a parabolic segment and has an area
equal to 2/3 of the length D of the chord multiplied by the "height" H
(the largest distance from the chord to a point on the arc).
To obtain the area of the "curved" triangle,
the above result for each of the three sides is simply added to,
or subtracted from, the area of the "straight" triangle.
Two millennia before "Calculus"...
The formula for the area of a parabolic segment predates "Calculus"
by almost two millennia.
The squaring of the parabola is attributed to
Archimedes of Syracuse, who first showed in an
"elementary" way (using an intuitive notion of limit) that a parabolic segment has 2/3
the area of its circumscribed
Archimedes
Triangle (whose sides are the chord and the two extreme tangents).
The surface area of a parabolic segment depends only
on the chord length (D) and the height (H) because
all parabolic segments of given D and H are obtained from each other by a
shear linear transform (about a direction parallel to the chord)
which preserves areas (by the first Cavalieri principle
concerning planar areas).
schmeelke (20020127)
If I know the dimensions of a cylindrical tank on its side
[the axis of revolution is horizontal] and can measure the depth of the liquid inside,
how do I calculate the volume of liquid present?
Let R be the radius of the tank, L its length and H the height of the liquid in it.
Consider a (circular) vertical crosssection of the tank and call
q the angle (from the center O of the circle)
between the vertical and the line OF,
where F is one of the two points where the horizontal line representing the surface of
the liquid meets the circle representing the wall of the tank [see figure at right].
We clearly have RH = R cos(q),
which means q is equal to
arccos(1H/R)
and is thus readily obtained using the proper inverse trigonometric function
on a scientific caculator.
[Angle q must be expressed in radians
and is thus between 0 and p;
don't forget to multiply a result in degrees by p/180,
if applicable.]
The surface area corresponding to the liquid in that cross section is obtained
as the difference of areas between a circular sector (a pie portion) and a triangle, namely
q R^{2} 
(RH)Ö(H(2RH)).
Just multiply this area by the tank's length L to obtain the formula for the volume V
of the liquid in the tank, namely:
V = L [ R^{2} arccos(1H/R) 
(RH)Ö(H(2RH)) ]
This formula is perfectly valid for the whole range of H (from 0 to 2R),
although the above "visual" explanation (involving the "difference of two areas")
assumed that H was less than R (tank at most halffull).
The formula could have been obtained symbolically without splitting cases.
I'll leave it up to you to "visualize" the other case
(where the area of a circular sector is to be added to that of a triangle of height HR),
should you feel the urge to do so.
Follow up : What if the tank is spherical?
[Volume of a Spherical Cap]
zchas40 (20020514)
[...] How much water is 2 cm in a hemispherical bowl 3 cm in radius?
Rick94602 (20020516)
[...] What's the volume of the bottom part of a sphere?
The vertical cross section [pictured at right] is the same
as above.
The portion of the sphere involved is called a spherical cap and
its volume V is obtained as:
ó õ
^{ q} _{0}
p R^{3}
sin^{3}a da
=
ó õ
^{ R} _{RH}
p
( R^{2}  u^{2 }) du
In other words,
V = p H^{2} ( R  H/3 ).
[ For the record, the surface area is
S = 2pRH. ]
If R is 3 cm and H is 2 cm, the volume of water in the "hemispherical bowl"
expressed in cm^{3}
(cc or mL, same thing) is 28p/3 or about 29.32 cc.
What about an ellipsoid of revolution on its side? [horizontal axis]
This will be useful to make short order of the next question...
We may observe that squeezing or stretching a spherical tank along any
horizontal direction turns the sphere into an ellipsoid with principal
radii R, R, and W (for some W).
Furthermore, for any height H of the liquid, the volume in the elliptical tank is
simply W/R times what it would be in the spherical tank.
In other words,
V = p H^{2} W ( 1  H/3R ).
Note that stretching along two horizontal directions would be just as easy,
if you ever need to derive a more general formula for any elliptic tank, as long as
one of the principal axes of the ellipsoid is vertical.
For the record, if R is the vertical radius and W and L are the horizontal ones,
the general formula is simply:
V = p WL ( 1  H/3R ) H^{2}/ R
Daniel Trottier (20020426; email)
[...] A horizontal vessel consists of a cylinder equipped with bell
ends. These caps sealing the ends of the cylinder are elliptical in shape.
(Removing the cylinder middle section, you would have an ellipsoid of revolution.)
What is the volume of liquid in the vessel, as a function of the height of the liquid?
If L is the length of the cylindrical section of radius R, and W is the "width"
of the end caps [W would be equal to R if those caps were spherical],
then the volume V corresponding to
a height H of liquid is simply obtained by adding the two volumes described in the
two previous articles, namely:
V = L [ R^{2} arccos(1H/R) 
(RH)Ö(H(2RH)) ] +
p H^{2} W ( 1  H/3R )
As usual, we assume that the result of an inverse trigonometric function (like arccos)
is given in radians.
Do not forget to multiply a result in degrees by p/180,
whenever applicable.
(20031208)
What's the volume of a truncated cylinder?
What's the volume of a [triangular] prism with tilted [nonparallel] bases?
The [American Heritage] dictionary defines a cylinder as
"the surface generated by a straight line
intersecting and moving along a closed planar curve".
When the planar curve is a polygon, this is also called a prism.
However, both names are more commonly used for
the solid bounded by two parallel planes and such a surface.
The two planar faces are called bases and the cylindrical surface
between the planes is called the lateral surface.
It turns out that a general expression can be given for the volume of such a solid
even when the
bases are not parallel to each other.
We only need to know the surface area (S) of a crosssection
[namely, the intersection with a plane perpendicular to the lateral direction]
and the average height (h) of such a solid.
The volume (V) of the solid is then simply given by:
V = h S
The only delicate part is to specify how the average height (h) is obtained.
Well, we only need to know the position of the centroid of the crosssection
(that would be the barycenter or center of mass of the surface if it was cut
from an homogeneous planar sheet).
This is particularly easy to find if the crosssection is triangular,
or if it is symmetrical enough
(circle, ellipse, regular polygon, parallelogram, etc.).
Draw a line parallel to the lateral direction through this centroid.
This line intersects the two bases at two points and the distance
between these two point is
the quantity h used in the above formula.
Of course, if either base is tilted by an angle q,
its surface area A is S/cos(q).
If the two bases are parallel to each other,
the distance d between their respective planes
happens to be h cos(q).
Thus we have
V = d A
This is a more commonly quoted formula, but it's only good for parallel bases...
(20040324) _{ } Conical Frustum
What's the volume (V) of the part of a cone between parallel bases of areas B and b ?
V = (H/3) [ B + b +
Ö
Bb
]
In this, B and b are the surface areas
of the two homothetic planar
bases (which need not have any specific shape),
whereas H is the solid's height
(the distance between the two planes of the bases).
The above French cartoon by Gotlib
[text by René Goscinny] is from a Dingodossier
(crazy file) about the flu, featuring
an average student without a first name:
l'élève Chaprot.
Chaprot's ability to repeat this formula flawlessly is
presented as a sure symptom of the flu...
Comical Frustum:
The formula was no longer taught systematically to French schoolchildren during the
golden age of French comic strips (spearheaded by Pilote, the weekly
magazine which published the above in 1967).
However, it had been a staple of French elementary mathematics at a time when
authors Goscinny (19261977)
and/or Gotlib
(Marcel Gotlieb, b.19340714)
were schoolboys...
Incidentally, René Goscinny
was born in Paris but was educated in
the French School of Buenos Aires, where his father (d.1942) taught mathematics.
The French Lycée of Warsaw (Poland) has been named after René Goscinny,
in part because of his Polish ancestry:
He was the grandson of Rabbi Abraham Goscinny of Warsaw.
René Goscinny was also the grandson of Abraham Beresniak,
who authored a 1939 Yiddish/Hebrew dictionary published in Paris.
In spite of this, Yiddish is not
one of the 107 languages currently spoken by Astérix,
the most famous of the 2120 characters created by René Goscinny,
in 18 different series
(a total of over 500 000 000 comic books sold).
The Polish roots of René Goscinny were pointed out to us by the current (2005)
Consul General of France in Los Angeles,
Philippe Larrieu,
formerly stationed in Warsaw (19941998). We thank his Deputy,
Olivier Plançon,
for many fun discussions about French comic books of the golden age.
(20040723) _{ } Volume of a Sphere,
using Cavalieri's Principle (1635)
How did they compute volumes before the advent of Calculus?
Here is our version of what's called the [second] principle of Cavalieri.
(The first principle of Cavalieri
deals with planar areas instead of volumes.)
If an horizontal plane always intersects two given solids
in sections of equal areas, then the solids have equal volumes.
Two such solids are said to be Cavalieri congruent.
This oldfashioned "principle" can be made to sound trivial nowadays
(the integrals of equal functions are equal) but it helped
define the very concept of integration.
Some direct applications are still interesting.
For example, we can deduce the volume of a sphere from the formulas giving the volume
of a cylinder and the volume of a cone:
The section of a sphere of radius R on a plane at a distance z (<R)
above its equatorial plane is a disk of area:
p ( R^{ 2}  z^{ 2 })
This is also clearly the area of a ring with outside radius R and
inside radius z,
which corresponds to a solid that consists of a cylinder of radius R and height R from
which a cone of the same base and height is removed
(the apex of the cone being on the
equatorial plane). The volume of this solid (2/3 the volume of the cylinder)
is thus the volume of the hemisphere.
For the whole sphere, this does gives the familiar formula of
which Archimedes of Syracuse was most proud of:
V_{sphere} = (4p/3) R^{3}
This was obtained two millenia before Cavalieri's principle got its name.
(20070921) _{ } Volume of a Tetrahedron
A tetrahedron's volume is ^{1}/_{6}
the determinant of 3 edges.
Proof : A tetrahedron ABCD can be
construed as a cone of apex D and base ABC.
As such, its volume is a third of the product of the base area by
the height, which is also the dot product of AD by the
vectorial area of ABC
(itself equal to half
the cross product of AB by AC). All told, the volume is:
V = ^{1}/_{6}
 ( AB ´ AC ) . AD 
= ^{1}/_{6}  det ( AB, AC, AD ) 
John Hanson (20090211; email)
Right Section of a Cone
What's the volume of a conical heap limited by a vertical plane?
A vertical cone of height H, radius R and apex (0,0,0)
has cartesian equation:
( x^{2} + y^{2 }) / R^{2} = z^{2} / H^{2}
On the vertical wall where x is a positive constant not exceeding R,
the above can be interpreted as the cartesian equation (in y and z) of the
hyperbola where the plane and the cone intersect.
Between the altitude
H [of the cone's base] and the altitude xH/R
[of the hyperbola's apex] the surface area of the hyperbolic segment is:
ó õ
^{ H/R} _{xH/R}
2
Ö
(zR/H)^{2}  x^{2}
dz
Srikanth (20040717; email)
_{ } Wedge of a Cone
What's the volume of the portion of a circular cone
included between two halfplanes?
In particular, when the intersection of the planes goes through
the cone's axis. [ As is the case for the red section
of the grey cone pictured at right,
seen from the direction (Ox) shared by the two halfplanes. ]
This solid may be seen as the difference of two cones with the same apex:
The larger cone is of height H. Its base is a halfcircle of area
pR^{2 }/ 2.
The height of the smaller cone is h = H cos q.
Its base is a segment of a conic section of area A
(a parabola when the plane's slope is H/R,
an ellipse for a less inclined plane, an hyperbola for a steeper one).
The volume (V) of the wedge thus reduces to the computation of the area A:
V = (H/3) [ pR^{2 }/ 2
 A cos q ]
With suitable choices of coordinate systems, a point of
coordinates (x,y) in the inclined plane has spatial coordinates
(x, y cos q, y sin q)...
Plug these into the cone's spatial equation to obtain the equation of
its intersection with the plane:
H^{2 }( x^{ 2} + y^{ 2 }cos^{ 2 }q )
= R^{2 }( H  y sin q )^{ 2}
The curve's apex (x=0, y>0) is at
y = r = [ (sin q) / H +
(cos q) / R ]^{ 1}
The segment's area A is thus given by either of the following expressions, where
x and y are positive quantities related to each other by the above equation:
A = 2
ó õ
^{ r} _{0}
x dy =
2
ó õ
^{ R} _{0}
y dx
For example, in the special case of a parabolic segment
(tan q = H/R) we have:
y = (R^{2}x^{2}) / 2R cos q
and
A = 2R^{2} / 3 cos q
Therefore, in this parabolic case, the volume boils down to:
V =
HR^{2} ( 3p
 4 ) / 18
[ Thanks to Steve Battison for catching an
embarrassing typo in the last step. ]
In the elliptic case (tan q < H/R)...
In the hyperbolic case (tan q > H/R)...
Generally, any conical wedge bounded by planar sections
of respective areas A_{0} and A_{1} [not necessarily joining
within the cone] and whose respective planes are at distances
h_{0} and h_{1} from the cone's apex,
has the following volume V:
V =
^{1}/_{3}  h_{0 }A_{0}
 h_{1 }A_{1 }
In this, the cone is understood to be singlesided [its surface is generated
by halflines originating at the apex] but it need not be a circular one...
suzjor (Riverside, CA. 20001124)
It takes you 6 hours to do a job. It takes a friend 3 hours to do the same job.
How long would it take both of you working together to do the job?
2 hours. In that amount of time,
you'll have completed 1/3 of the job and your friend 2/3 of it.
Of course, the assumption is that the "job" is of such a nature
that it can be "distributed" efficiently.
This result is obtained by adding the rates [the number of jobs per hour]
of both workers to obtain the rate for the team.
In this case, your rate is 1/6 and that of your friend is 1/3,
so the combined rate of the team is 1/6+1/3, or 1/2; the team does half a job in an hour,
so it takes 2 hours for the whole job.
seriouslyman (20001125)
A clerk is assigned a job that she can complete in 8 hours.
After she has been working for 2 hours, another clerk,
who is able to do the job in 10 hours, is assigned to help her.
How long will it take to complete the job?
We assume that the "job" is of such a nature that it can be distributed between the
two clerks. (A modeling assignement would not be such a "job"...)
The first clerk completes 1/8 of a job per hour, the second 1/10 of a job per hour
(if you want to be silly, you could say their respective "speeds" are 1/8 jph and 1/10 jph).
After t hours (with t>2), the first clerk will have completed a fraction t/8 of the job,
and the second clerk a fraction (t2)/10. The whole job will be completed when t is such that:
1 = t/8 + (t2)/10 Solve for t and you get t=16/3 hours, that's 5 hours and 20 minutes.
Check your answer:
The first clerk has worked 16/3 hours and has completed 16/24=2/3 of the job.
The second clerk has worked 16/32=10/3 hours and had completed 1/3 of the job.
(S. H. of Mays Landing, NJ. 20001123twice)
Let us say that you want to prepare an a % solution [of peroxide],
using 3% and 30% solutions. In what proportion should you prepare the mix?
(Here, a is a given number, between 3 and 30.)
To get an a% solution from a strong solution at p% and a weak one at q%,
you should mix (pa) parts of the weak one and
(aq) parts of the strong one.
This will indeed give you (pq) parts of a solution rated at:
a =
(pa)q + (aq)p
(pa) + (aq)
For example, if p=30% and q=3%, you obtain a solution at a = 12%
by mixing 3012 = 18 parts of the stronger solution with
123 = 9 parts of the weaker one
(that's 2 parts of the strong solution and 1 part of the weak one,
if you must have proportions expressed in lowest terms).
(20001111) Alcohol Content
Relating ethanol content by volume (x) to content by mass (y).
Metallic alloys and almost all other chemical mixtures are
normally rated by mass.
The percentage given is the ratio of the mass of the ingredient
of interest to the mass of the whole mixture.
Alcohol solutions, on the other hand, are usually rated by volume.
The volumes involved are understood to be volumes of pure alcohol
or pure water before mixing. Therefore, the
basic relevant computations
don't depend on the wellknown fact that
mixing one liter of pure alcohol with one liter of water yields only
1.92 liters of liquid!
Let d be the relative density of alcohol with respect to water.
The value of d is practically equal to 0.79 (nearly 4/5).
More precisely, the relative density d is the ratio of the
absolute densities of
alcohol
and water:
At 20°C, d = (789.45 g/L) / (998.2071 g/L) = 0.79087
At 25°C, d = (785.22 g/L) / (997.0479 g/L) = 0.78754
The relative density is exactly 0.79
at a temperature of about 21.3°C.
A rating of x by volume (ABV = alcohol by volume) would be obtained
by mixing x volumes of pure ethanol with (1x) volumes of water.
This yields a liquid whose rating by weight is y
(ABW = alcohol by weight) where:
y = xd / (xd + 1  x)
or
1/y = 1 + 1/xd  1/d
Rating alcohol (d = 0.79) by volume (x) or by weight (y)
(1/y  1) d = (1/x  1)
For a dilute alcohol (x and/or y are small) each bracket is well approximated
by its first term and, therefore, y is approximately equal to xd.
For example, 5% by volume is approximately 4% by weight (using d = 4/5).
Of course, 100% by volume is 100% by weight.
The above formula does hold for x = y = 1
(the approximation y = xd isn't applicable;
it would be off by 20% or 21%).
US proof is twice the percent of alcohol by volume :
Vodka is normally 40% by volume (sold as "80 proof", in the US)
which corresponds
to 34.5% by weight (using d = 0.79).
Everclear is 95% by volume
(190 proof) which is 93.75% by weight. For an Nproof spirit:
x = N / 200
y = 0.79 / ( 200 / N  0.21 )
(V. V. of Emeryville, CA.
20001111)
How much pure ethanol must a nurse add to 10 cubic centimeters of a 60% ethanol solution
to strengthen it to a 90% alcohol solution?
The computation would be easy if the nurse
had 6 cc of pure alcohol and
4 cc of water instead of the 10 cc mix.
She would simply need to add z cc of pure alcohol to
make the volume rating (defined above) equal to 0.9:
0.9 = (z + 6) / (z + 6 + 4)
Solving for z, we obtain z = 30 cc.
That's what they must teach in nursing school
(possibly with the nice shortcut that's valid either for
ratings by mass or for mixing liquids whose volumes don't change significantly upon mixing).
What a stellar nurse would do :
Because water and ethanol shrink significantly upon mixing,
there's more stuff (both water and ethanol) in her 10 cc mixture than would
have been obtained from mixing 6 cc of ethanol with 4 cc of water...
Let's assume that the original solution had been
prepared at 20°C and that it's being adjusted at the same temperature.
We use a single
data source, according to which, at 20°C :
The specific gravity of pure water is 0.99823
The specific gravity of pure ethanol is 0.78934
The relative density is thus d = 0.78934 / 0.99823 = 0.79074
By weight, a 60% ABV solution is y = 1 / (1+2/3d) = 54.2566...%
The specific gravity of a 54% solution by weight is 0.90485
The specific gravity of a 55% solution by weight is 0.90258
The interpolation 0.90258(100y54)+0.90485(55100y) is 0.90427
The unmixed average 0.99823(0.4)+0.78934(0.6) is 0.872896
The mix is 0.90427/0.872896 = 1.035942 denser than its components.
Increasing her previous estimate by that last factor, the nurse will add:
(30 cc) 1.035942 = 31.08 cc
We leave it to the reader to work out what the volume of the resulting 90% ABV solution will
be at 20°C (no, it's not 41.08 cc).
dva1270 (20020329)
_{ } Averages Galore...
A vehicle goes from A to B at 60 mph, and returns at 40 mph.
Why is the average rate of speed
equal to 48 mph? [not 50 mph]
That's because average speeds are harmonic, not arithmetic... Read on.
The average of several things is whatever single value you could replace all of these
things with, and still obtain the same effect.
Obviously, the way you would actually compute such an average may very well
depend on exactly which type of effect you are interested in.
In various cases, there may be hidden assumptions,
which must be made explicit before a reliable computation can be made:
Arithmetic Mean:
For example, if you are receiving two checks, the only thing you normally care about
is the total amount of money you are receiving.
In this case, it does make sense to consider that an average check is equal to
half the sum of both, because two such average checks would have exactly the
same effect on your bank account as the two checks you actually got.
This is called an arithmetic average and it's, by far, the most common form
of averaging.
It's clearly not the only one, though...
Harmonic Mean:
In the case of speeds, the important thing is how long [how much time]
it takes to accomplish a journey whose legs are traveled at different speeds.
In other words, the average speed has to be whatever uniform speed
would allow the journey to be completed in the same amount of time as the actual one
[where the speeds on different legs may have been widely different].
This is clearly equal to the total distance divided by the total duration of the journey.
That number is obtained as the socalled harmonic average of the speeds,
namely the speed whose reciprocal is the arithmetic average
of the reciprocals of the speeds involved.
Using the numerical example in the question, 48 is the harmonic average of
60 and 40, because:
1/48 = ( 1/60 + 1/40 ) / 2
Geometric Mean:
The geometric mean of two nonnegative numbers is the square root
of their product.
More generally, geometric averaging is appropriate
for successive relative increases:
Consider, for example, a firstyear increase of
10% (a factor of 1.1),
a secondyear increase of 21% (a factor of 1.21),
and a thirdyear increase of 700% (a factor of 8).
Over three years, the increase corresponds to a factor of
1.1 ´ 1.21 ´ 8 = 10.648.
The yearly average is the cube root of that, namely
a factor of 2.2, corresponding to an average yearly increase of 120%.
The logarithm of a geometric average is the arithmetic average of the logarithms:
Log(2.2) = [ Log(1.1) + Log(1.21) + Log(8) ] / 3
The fact that people often express a relative increase in terms of a
percentage difference is best discarded,
except for the input and output of data.
Quadratic Mean:
This is the [positive] quantity whose square is the arithmetic mean of
the squares of the quantities under consideration.
It's commonly called the root mean square (RMS),
as it's the square root of the mean of the squares.
This averaging is common in physics.
For example, the RMS speed of a [large] set of molecules may be used to
define
their temperature (which is classically proportional
to their average kinetic energy).
Also, the RMS value [over time] of the current through a pure resistor
equals the continuous current that would dissipate the same heat as the
varying current observed.
Hölder Mean:
The German mathematician
Otto
L. Hölder (18591937) investigated, for any given exponent p,
the quantity H whose pth power is the arithmetic mean of the pth powers
of the quantities under consideration.
(For arbitrary exponents,
this must be restricted to nonnegative quantities.)
Exponents 1, 1, 0 and 2 correspond, respectively, to the aforementioned
arithmetic, harmonic, geometric and quadratic means._{ }
The case p = 0 for the geometric mean
is defined by continuity as the limit
of the pexponent Hölder mean H_{p},
when p tends to zero.
If p < 0 and anya_{i} is zero,
then H_{p} = 0 (defined by continuity).
H_{p}(a,b)
H_{p}(a,b) = ab
For an infinite exponent,
the Hölder mean of a finite set of nonnegative quantities is defined as either
their maximum (p = +¥)
or their minimum (p = ¥).
Generalized Mean:
A bijective function f can be introduced so that the
following relation defines the "mean" m of n quantities
a_{1}, a_{2} ... a_{n}:
f (m) =
(1/n) [ f (a_{1 }) +
f (a_{2 })
+ ... + f (a_{n }) ]
f (x) is x for the ordinary arithmetic mean,
x^{2} for the quadratic mean,
1/x for the harmonic mean,
Log(x) for the geometric mean [only when x>0], etc.
Averaging over a spheroid: Nautical mile as an "average" minute of latitude.
Batting Averages:
...Farey Series, etc.
(20030108) 30.436875 or 30.458729474253406983...
What's the average number of days in a Gregorian month?
This is a good way to expand the discussion in the
previous article.
The average or mean value of anything is often ultimately defined
as the "expected value" of that thing according to some probability distribution.
For most pratical purposes, the distribution to use is clearly understood and
may be left unspecified.
For example, when asked out of context about the average of
28, 29, 30 and 31, we would normally assume that these 4 numbers
are equally likely and would readily equate their mean (their expected value)
with their arithmetical average (namely 29½).
On the other hand, this approach is clearly unacceptable if we have to work out
the mean number of days in a month, since durations of 28 or 29 days
should clearly weigh less than months of 30 or 31 days,
on account of their lesser likelihood.
Our modern secular calendar (the Gregorian calendar) repeats with a period
of 400 years which includes 97 leaps years of 366 days and 303 regular years
of 365 days.
Therefore, this period consists of 146097 days, or 4800 months.
The average number of days in a month would thus seem to be
146097/4800, namely 30.436875.
But is is really so?
Well, yes and no. The question is not precise enough to have a definite answer.
It depends on exactly what type of event you expect the "mean" to predict.
If each month is equally likely, the above is indeed the correct answer.
However, it would be at least as reasonable to assume that what you really want is the
average duration of the current month for a random Gregorian date.
In this case, each day is equally likely in the Gregorian cycle of 146097 days and the
result must be somewhat higher than the above, because longer months are more likely.
More precisely, in a Gregorian cycle there are 97 months of 29 days (1 per
leap year), 303 months of 28 days (1 per regular year), 1600 months of 30
days (4 in any year) and 2800 months of 31 days (7 in any year).
Therefore, there are 2813 days (97 times 29) which fall in a month of 29 days,
8484 days (303 times 28) falling in a month of 28 days, 48000 days falling
in a month of 30 days, and 86800 (2800 times 31) falling in a month of 31 days.
The "mean" is obtained by summing up all possible values "weighted" with their
respective probabilities.
For example, the value 29 is weighted by the
probability of 29, which is 2813/146097 (roughly 1.925433%).
Work it out and what you obtain for the mean length of a month, in that sense,
is 4449929/146097, or about 30.4587294742534...
This is slightly more than our previous result, as predicted.
Neither answer is better that the other; they are just answers to different questions.
Which answer you pick depends on which question you feel is
"intended" when pople ask about the mean length of a month.
Given enough time to think it over,
most mathematicians would probably find the latter interpretation of the question
more "natural", but they would still acknowledge that other interpretations
are better adapted to specific contexts.
For example, our earlier value of 30.436875 days is the only
acceptable conversion factor between durations expressed in days and
durations expressed in average calendar months;
any other number leads to a systematic bias for long durations...
(20090208) ArithmeticGeometric Mean (AGM)
The AGM of a and b is the AGM of
(ab)^{½} and (a+b) / 2
The AGM of two positive numbers can be obtained by replacing them iteratively by their
geometric mean and their arithmetic mean until those
two bounds coincide at a given precision.
For example, starting with 4 and 9:
WHILE A+B <> A+A
X=SQRT(A*B) : B=(A+B)/2 : A=X
WEND
The above is an efficient algorithm with a quadratic convergence
(i.e., the number of correct digits roughly doubles with each iteration).
The arithmeticgeometric mean is tightly related to the complete elliptic integral
of the first kind (K) as follows.
(See also: Gauss's constant.)
K(k) ^{ } =
ò
^{ p}/_{2}
_{0}
dj
^{ } =
^{p}/_{2}
Ö
1  k^{2} sin^{2 }j
agm ( 1k , 1+k )
The above expression of K is called its normal trigonometric form.
The substitution x = sin j.
shows it to be equivalent to the Legendre normal form
that appears elsewhere on this site.
Conversely,
agm ( a,b ) =
p ( a + b )
/
4 K( a  b /
a + b )
wankman (20020410)
How far away is the ocean horizon line?
Is there a formula for figuring that out?
Suppose the Earth is a perfect sphere of radius R (that's not quite true, but close enough).
If your eyes are at a height H above the surface of the ocean,
then the horizon is at a distance D which is such that your line of sight
is tangent to the sphere at that distance. See figure at right.
So, you have a right triangle whose vertices consist of your own eye,
a point on the horizon and the center of the Earth.
The hypotenuse is R+H and the sides are R and D.
Therefore, (R+H)^{2} = R^{2} + D^{2},
which boils down to:
D^{2} = (2R+H) H
In practice, H is much smaller than the diameter of the Earth (2R),
so that (2R+H) and 2R are virtually the same and we have
2RH » D^{2}.
The distance D to the horizon seen from an altitude H is thus^{ }:
D »
Ö
2RH
In this, R is the (conventional) radius of the Earth:
R = 6371000 m or R = 20902231 ft
(choose whichever unit of length you use to express H
in order to have the result D expressed in the same unit).
For example, at an altitude of 2 m [standing on a small boat],
the horizon is at 5048 m
(about 5 km, 3.2 statute miles, or 2.7 nautical miles).
If you climb on top of a tall mast at 20 m [ten times higher],
the horizon is at 15964 m [3.16 times farther].
From a hilltop on the seashore at 200 m,
the horizon is 10 times farther than if you are standing on your small boat.
From a high mountain at 2000 m,
you could see the ocean 31.6 times farther than when standing on a small rock on the beach:
The horizon line is then almost 100 miles away! [86 nautical miles]
For a quick estimate, compute the square root of your altitude in feet.
Add 20% or 25% to that result [more precisely 22.455%] to
obtain the distance to the horizon in statute miles
(for a result in nautical miles, add 5% instead, or
6.41% to be more precise).
For example, if your eyes are 9 feet above the ocean
[the square root is 3], the horizon is about 3.7 miles away
(roughly 3.2 nautical miles).
The distance to the horizon in kilometers is roughly 3.5 times
the square root of the altitude of your eyes in meters
(a more precise coefficient is 3.5696).
(R. S. of Austin, TX. 20001025)
Distances along a great circle How can I find the distance (in feet) between two exact locations by
entering decimal lat/long coordinates to 6 place values?
lancelizabeth^{ }
(20010711)
Using latitude and longitude, what is the formula that calculates
the distance [in meters] between two points on Earth?
Assuming the Earth to be a perfect
sphere whose conventional radius is
R = 6371000 m
» 20902231 ft,
a point of latitude j
and longitude q has the following
cartesian coordinates
(in a wellchosen coordinate system) :
[x, y, z] =
[ R cos j cos q ,
R cos j sin q ,
R sin j ]
In this, a northern latitude is positive, a southern latitude is negative.
An eastern longitude is positive and a western one is negative.
(You could use any other convention, provided you do so consistently.)
The conventional radius of the Earth (R = 6371 km)
is very close to the
average distance from the Earth's center
to a point "at sea level".
If you have two such points, the scalar product ("dot product") of two such vectors
gives you R^{2} times the cosine of the angle d between
them (along a great circle, which is the shortest route on a sphere's surface).
In other words, two points whose spherical coordinates
(latitude,longitude) are respectively
(j,q) and
(j',q'^{ })
are separated by an angle d given by:
cos d =
cos j cos q
cos j' cos q'
+
cos j sin q
cos j' sin q'
+
sin j
sin j'
Angular separation (d)
of distant points on a sphere :
cos d =
cos j cos j'
cos (q'q)
+
sin j sin j'
From the above cosine of d,
you may obtain the angle d by using
the inverse trigonometric function on a scientific calculator.
[However, you may prefer to use the better formula qelow,
which remains accurate for nearby points!]
If d is expressed in radians,
the distance between the two points along a great circle
is simply equal to dR.
If d is in degrees,
the distance is pdR/180.
A Better Formula (especially for small distances):
The above formula is not wellsuited numerically
to the computations of distances between
nearby points, because d
is small and cos d is so close to unity
that the use of the reverse trigonometric function arccos
would cause an unacceptable loss of accuracy in the final result.
To skirt the difficulty, we may use the following equivalent formula,
which turns an accurate knowledge of
(j'j)
and
(q'q)
into a similarly accurate floatingpoint value for d:
Angular separation (d)
of distant or nearby
points on a sphere :
sin^{2}_{ }
d
= sin^{2 }_{ }
j'j
cos^{2 }_{ }
q'q
+ cos^{2 }_{ }
j'+j
sin^{2 }_{ }
q'q
2
2
2
2
2
To obtain this relation from the above, the identity
cos x = 1  2 sin^{2 }(x/2) may be used,
which allows the following transformation of our original formula:
1  2 sin^{2 }(d/2)
= cos j cos j'
{ 1  2 sin^{2 }([q'q]/2) }
+ sin j sin j'
= cos (j'j)
 2 sin^{2 }([q'q]/2)
cos j cos j'
= 1  2 sin^{2 }([j'j]/2)  2 sin^{2
}([q'q]/2)
cos j cos j'
Therefore:
sin^{2 }(d/2)
= sin^{2 }([j'j]/2)
+ sin^{2 }([q'q]/2)
cos j cos j'
= sin^{2 }([j'j]/2)
+ sin^{2 }([q'q]/2)
{ ½ cos(j'+j)
+ ½ cos(j'j) }
(Ukulele8421. 20100209)
Distance through the Earth
What's the length of a straight line
from Grand Rapids to Melbourne?
In the above spherical approximation, the two cities
are adequately represented by two points (at zero elevation) on a perfect
sphere of radius R = 6371 km.
Each point is located in space by the cartesian coordinates of a vector:
U =
(R+h) [ cos j cos q ,
cos j sin q ,
sin j ]
The distance d
between the two points through the Earth is simply
the length of the difference between two such vectors.
The square of that distance is:
 U  U' ^{ 2}
= 2 R^{2}  2 U . U'
(assuming h=0)
That expression would be easy enough to spell out but it suffers from the
same flaw as the first formula of the
previous section
(i.e., numerical loss of precision at small distances).
Instead, it's much better to compute the angular separation d
between the two cities (using the
second formula given above) with the following
relation, obtained from elementary geometry (in the plane of the great circle):
d = 2 R sin (d/2)
Happily, the formula
we are recommending gives sin (d/2) directly,
via its square!
In the case of Grand Rapids and Melbourne, we obtain numerically:
sin (d/2) = 0.94462778
Multiply this into the diameter of the Earth (2R) to obtain:
d = 12036.448 km.
At this level of precision (1 m) the difference of altitudes come
back into play:
If Grand Rapids was at the same altitude (31 m) as
Melbourne, we would obtain an exact result (for a spherical Earth)
by using R = 6317.031 km. This yields
d = 12036.506 km. The difference in elevations
of 164 m yields (in the main) an additional
correction equal to
that difference multiplied into the negative cosine of the angular separation
(d = 141.6877°).
All told, the distance between the two cities would be 12036.635 km
on a spherical Earth.
The correction in distance due to the oblateness of the Earth varies between
a negative minimum of 28.5 km and a positive maximum
of +14.3 km obtained respectively when the two cities are
located either on the two poles or on opposite points of the equator.
The true oblateness correction (7.372 km)
can be worked out as follows:
Full precision, using the Reference Ellipsoid (IUGG, 1980) :
In geography, latitude is
always defined as the angle
j
from the plane of the equator to the local vertical
(northward angles being positive).
If the surface of the Earth is approximated by a
spheroid of equatorial radius
a = 6378137 m and polar radius
b = 6356752.3141 m, then the local vertical
is perpendicular to the surface of that ellipsoid. It's a
simple exercise in calculus
to express the cartesian coordinates
of the [elliptic] meridian as functions of the parameter
j:
Geodetic
Latitude (j) and Elevation (h).
Cartesian coordinates of a point of elevation h
at [geodetic ] latitude j
and longitude q :
x =
( a^{ 2}/ R_{j} + h )
cos j
cos q
y =
( a^{ 2}/ R_{j } + h )
cos j
sin q
z =
( b^{ 2}/ R_{j } + h )
sin j
where R_{j}^{2} =
a^{2} cos^{2} j
+
b^{2} sin^{2} j
Those formulas can be used to compute the components of the vector which
goes from Grand Rapids to Melbourne.
The straightline distance is its length, namely :
Ö
= 12029.263 km
(x'x)^{2} +
(y'y)^{2} +
(z'z)^{2}
Incidentally, the above formulas also give the relation
(at zero elevation) between the geodetic
latitude j
and the (useless) geocentric latitude
j_{0 }:
tg j_{0}
=
( b^{ 2}/ a^{ 2 })
tg j
To see that, notice that we have
tg j_{0} = z/x when
q = 0 and h = 0.
(Heather of Canada. 20001106)
About geodetic latitude...
How do I find the radius of the Earth at 51° North?
The term is ambiguous. I'll compute two related quantities:_{ }
The distance from the center of the Earth to a point at a
latitude L of 51°.
The radius of that latitude's parallel.
Before we do a precise computation, let's do the usual rough one by considering
that the Earth is a perfect sphere whose radius is the "conventional" radius
which is defined to be exactly R=6371000m. With this approximation,
we have:_{ }
A (constant) distance to the center equal to R, namely 6371 km.
A parallel of radius R cos(L), or about 4009.4 km when L is 51°.
A much better approximation is to consider that the Earth matches exactly
the regular shape against which its irregularities are charted.
That shape is called the "reference ellipsoid" and its dimension have been
precisely defined once and for all by the
IUGG in 1980:
The meridian is a perfect ellipse whose equatorial radius is exactlya = 6378137m and whose polar radius is "approximately"
b = 6356752.3141m. So far so good.
Now, what does a latitude L correspond to on that ellipse?
It is not the socalled "geocentric" latitude
(the angle between the line to the center of the Earth and the plane of the equator).
Instead the latitude L is the geodetic latitude:
the angle between your local vertical (which is the line perpendicular to the
ellipsoid's surface) and the equatorial plane.
L is the latitude you would get from local astronomical observations.
With the spheric approximation, we did not have to worry about this fine point,
because both definitions amount to the same thing.
In the case of the ellipse they are slightly different.
How different? Well, if we call C the geocentric latitude,
the difference LC is given by:
L  C = 692.74" sin(2L)  1.16" sin(4L).
For L = 51°, LC is thus about 678.07378" or about 0.1888358°.
So the geocentric latitude C is about 50.811646°.
Let x be the radius of your parallel and y its distance to the equatorial plane.
We have y = x tan(C) = (1.22663011194) x,
while the equation of the ellipse gives you (with the above values
for a and b)
x^{2}/a^{2 }+ y^{2}/b^{2} = 1.
Replace y by x tan(C) and solve for x to obtain the value of x,
whereas the distance to the center of the Earth is simply x/cos(C)
(its square is x^{2}+y^{2}):_{ }
The distance to the center of the Earth is 6365264.58 m (that's about 5.7 km
less than the conventional radius of the Earth R).
The radius of the parallel is x = 4022.031 km
(that's about 12.6 km more than what we got from the spherical approximation).
yonda1234 (20010428)
Applying Girard's Formula
How can I find the area of a spherical polygon formed by a series of points
on the Globe, given by their latitudes and longitudes?
Here, we assume the Globe is a perfect sphere. For better precision,
we should consider it to be an oblate ellipsoid instead.
I've worked out elsewhere
the formula to use for the surface area of a patch in this case.
The result is surprisingly simple for a polygonal patch where vertices are connected
by some special curves on the ellipsoid.
A basic result of spherical geometry (called Girard's Theorem)
states that the surface area of a spherical triangle on a sphere of radius R
is equal to R^{2}(A+B+Cp) ,
where A, B and C are the inside angles of the triangle.
You could dissect your polygon into triangles and add up the results obtained from
Girard's formula, but there's a more practical way to proceed
(whose validity may be established using such a virtual dissection):
Consider the sum S
of all the inside angles at each of the n vertices of your polygon
(see below for a recipe to compute these angles from your list of spherical coordinates).
The surface area of your polygon is simply:
R^{2} [ S  (n2)p ]
The polygon need not be convex,
it's only assumed that the polygonal line does not intersect itself
(or else, the formula would only hold for a critical matching definition of
inside angles and inside area; see footnote below).
Therefore, the only real problem is to compute the relevant angles from your
latitude/longitude list.
Well, on a sphere of unit radius, a point of latitude u and longitude v
has cartesian coordinates (cos(u)cos(v), cos(u)sin(v), sin(u)).
To compute the angle at point B in your polygon,
you need to consider this vector for B as well as for the previous point (A, say)
and the next one (C, say). The cross products BxA and BxC give you the directions
perpendicular to the planes of the dihedral angle you're after.
Divide these two vectors by their length to obtain two unit vectors U and V.
The scalar product of U and V is the cosine of the angle x you want
(if your polygonal line is an approximation to a smooth curve,
there won't be too much bending at each corner and x will be close to
p).
Note, however, that the positive angles x and
(2px)
have the same cosine; you must choose whichever corresponds to the
inside of your polygonal line. That's it. I hope this helps...
Footnote:
The only consistent way to define which of the two possible angles is the inside
one is to always put it on the same side with respect to
your progression along the polygon line (the usual convention is that the inside
is to your left).
That's critical when the line selfintersects,
in which case the inside angle occasionally faces "outside".
(20021201) Kepler's Third Law
The relation between distances and periods of orbiting bodies.
Johannes Kepler (15711630)
was the first to observe that the square of the period of
a planet around the Sun is proportional to the cube of the size of its orbit.
More precisely, the size of an elliptical orbit should be defined as its
major radius (a),
also called its semimajor axis (which happens to be equal
to the average distance of the planet to the Sun).
This is now known as Kepler's Third Law of planetary motion.
It occurred to Johannes Kepler on March 8, 1618; more than twelve years
after he had formulated his first two laws.
[He did not accept his own idea until better computations finally allowed him to formulate
the third law, on May 15, 1618.]
Kevin Brown
remarks that the introduction of logarithms
by John Napier (in 1614) is likely to have been instrumental in Kepler's
discovery of the third law...
Kepler discovered his first two laws in 1605 and published both in 1609:
The first law is that orbits are ellipses of which the Sun is a focal point,
whereas the second law states that a straight line from the Sun to the planet
"sweeps out equal areas in equal time intervals".
Oddly enough, the second law was actually discovered first...
Kepler's third law was later explained by Newtonian mechanics
and may be more precisely
expressed in terms of the following relation, which involves Newton's constant of gravity (G),
the masses of the two bodies that orbit each other (M and m),
their average distance
(R = the major radius of the orbit of one body,
if the other is considered fixed)
and the period of revolution (T):
4 p^{2} R^{ 3}
= G (M + m) T^{ 2}
This relation holds for SI units,
or within any other consistent system of physical units.
However, it's worth noting that, for the motion of a small mass around the Sun,
we may choose to express distances in astronomical units and times in
sidereal years, which makes the relation boil down to: