RosettaCodeData/Task/Arithmetic-geometric-mean/Oforth/arithmetic-geometric-mean-1.fth

: agm   \ a b -- m
   while( 2dup <> ) [ 2dup + 2 / -rot * sqrt ] drop ;