RosettaCodeData/Task/Arithmetic-geometric-mean/BASIC/arithmetic-geometric-mean-2...

14 lines
210 B
Plaintext

*FLOAT 64
@% = &1010
PRINT FNagm(1, 1/SQR(2))
END
DEF FNagm(a,g)
LOCAL ta
REPEAT
ta = a
a = (a+g)/2
g = SQR(ta*g)
UNTIL a = ta
= a