8 lines
199 B
Text
8 lines
199 B
Text
F agm(a0, g0, tolerance = 1e-10)
|
|
V an = (a0 + g0) / 2.0
|
|
V gn = sqrt(a0 * g0)
|
|
L abs(an - gn) > tolerance
|
|
(an, gn) = ((an + gn) / 2.0, sqrt(an * gn))
|
|
R an
|
|
|
|
print(agm(1, 1 / sqrt(2)))
|