Anonymous user
Arithmetic-geometric mean: Difference between revisions
Adding SequenceL Version
(Adding SequenceL Version) |
|||
Line 2,126:
<pre>
1.456791
0.847213
</pre>
=={{header|SequenceL}}==
<lang sequencel>import <Utilities/Math.sl>;
agm(a, g) :=
let
iota := 1.0e-15;
arithmeticMean := 0.5 * (a + g);
geometricMean := sqrt(a * g);
in
a when abs(a-g) < iota
else
agm(arithmeticMean, geometricMean);
main := agm(1.0, 1.0 / sqrt(2));</lang>
{{out}}
<pre>
0.847213
</pre>
| |||