8 lines
152 B
Crystal
8 lines
152 B
Crystal
def agm (a : Float64, g : Float64)
|
|
while a - g >= Float64::EPSILON
|
|
a, g = (a + g) / 2, Math.sqrt(a * g)
|
|
end
|
|
g
|
|
end
|
|
|
|
p agm(1, 1 / Math.sqrt(2))
|