15 lines
257 B
Text
15 lines
257 B
Text
|
print agm(1, 1/sqr(2))
|
||
|
|
print using("#.#################",agm(1, 1/sqr(2)))
|
||
|
|
end
|
||
|
|
|
||
|
|
function agm(a,g)
|
||
|
|
do
|
||
|
|
absdiff = abs(a-g)
|
||
|
|
an=(a+g)/2
|
||
|
|
gn=sqr(a*g)
|
||
|
|
a=an
|
||
|
|
g=gn
|
||
|
|
loop while abs(an-gn)< absdiff
|
||
|
|
agm = a
|
||
|
|
end function
|