13 lines
253 B
Text
13 lines
253 B
Text
print agm(1, 1/sqr(2))
|
|
print using("#.#################",agm(1, 1/sqr(2)))
|
|
|
|
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
|