Arithmetic-geometric mean: Difference between revisions

m
Line 1,002:
-- sufficiently close, as determined by "eq".
agm :: (Floating a) => a -> a -> ((a, a) -> Bool) -> a
agm a g eq = snd .$ head . dropWhile (not .until eq) $ iterate step (a, g)
where step (a, g) = ((a + g) / 2, sqrt (a * g))
 
Anonymous user