Arithmetic-geometric mean: Difference between revisions

LONG REAL a := x, g := y;

LONG REAL epsilon := a + g;

LONG REAL next a := (a + g) / LONG 2.0, next g := long sqrt (a * g);

LONG REAL next epsilon := ABS (a - g);

WHILE next epsilon < epsilon

DO

print ((epsilon, " ", next epsilon, newline));

epsilon := next epsilon;

a := next a; g := next g;

next a := (a + g) / LONG 2.0; next g := long sqrt (a * g);

next epsilon := ABS (a - g)

OD;