11 lines
201 B
Text
11 lines
201 B
Text
|
$define EPS = 1e-14
|
||
|
|
|
||
|
|
local function agm(a, g)
|
||
|
|
while math.abs(a - g) > math.abs(a) * EPS do
|
||
|
|
a, g = (a + g) / 2, math.sqrt(a * g)
|
||
|
|
end
|
||
|
|
return a
|
||
|
|
end
|
||
|
|
|
||
|
|
print(agm(1.0, 1.0 /math.sqrt(2)))
|