|
In mathematics, the arithmetic-geometric mean (AGM) of two positive real numbers x and y is defined as follows: First compute the arithmetic mean of x and y and call it a1. Next compute the geometric mean of x and y and call it g1; this is the square root of the product xy: Then iterate this operation with a1 taking the place of x and g1 taking the place of y. In this way, two sequences (an) and (gn) are defined: These two sequences converge to the same number, which is the arithmetic-geometric mean of x and y; it is denoted by M(x, y), or sometimes by agm(x, y).
[edit] ExampleTo find the arithmetic-geometric mean of a0 = 24 and g0 = 6, first calculate their arithmetic mean and geometric mean, thus: and then iterate as follows:
The first four iterations give the following values:
The arithmetic-geometric mean of 24 and 6 is the common limit of these two sequences, which is approximately 13.45817148173. [edit] PropertiesM(x, y) is a number between the geometric and arithmetic mean of x and y; in particular it is between x and y. If r > 0, then M(rx, ry) = r M(x, y). There is a closed form expression for M(x,y): where K(x) is the complete elliptic integral of the first kind. The reciprocal of the arithmetic-geometric mean of 1 and the square root of 2 is called Gauss's constant. named after Carl Friedrich Gauss. The geometric-harmonic mean can be calculated by an analogous method, using sequences of geometric and harmonic means. The arithmetic-harmonic mean can be similarly defined, but takes the same value as the geometric mean. [edit] Implementation in PythonThe following example code in Python computes the arithmetic-geometric mean of two positive real numbers:
from math import sqrt
def avg(a, b, delta=None):
if None==delta:
delta=(a+b)/2*1E-10
if(abs(b-a)>delta):
return avg((a+b)/2.0, sqrt(a*b), delta)
else:
return (a+b)/2.0
[edit] See also[edit] References
Directorio de Enlaces Directorio dmoz Directorio espejo dmoz Pedro Bernardo |