Arithmetic-geometric mean/Calculate Pi: Difference between revisions

Content added Content deleted

Inline

@@ Line 1,405: / Line 1,405: @@
 multi sqrt(Int $n) {
-    my $guess = 10**($n.chars div 2);
+  (10**($n.chars div 2), { ($_ + $n div $_) div 2 } ... * == *).tail
-    my $iterator = { ( $^x   +   $n div ($^x) ) div 2 };
-    my $endpoint = { $^x == $^y|$^z };
-    return min (+$guess, $iterator … $endpoint)[*-1, *-2];
 }
 multi sqrt(FatRat $r --> FatRat) {
-    return FatRat.new:
+  return FatRat.new:
     sqrt($r.nude[0] * 10**(number-of-decimals*2) div $r.nude[1]),
 **number-of-decimals;
@@ Line 1,420: / Line 1,417: @@
 my FatRat $g = sqrt(1/2.FatRat);
 my $z = .25;
 for ^10 {
-    given [ ($a + $g)/2, sqrt($a * $g) ] {
+  given [ ($a + $g)/2, sqrt($a * $g) ] {
-	$z -= (.[0] - $a)**2 * $n;
+    $z -= (.[0] - $a)**2 * $n;
-	$n += $n;
+    $n += $n;
-	($a, $g) = @$_;
+    ($a, $g) = @$_;
-	say ($a ** 2 / $z).substr: 0, 2 + number-of-decimals;
+    say ($a ** 2 / $z).substr: 0, 2 + number-of-decimals;
-    }
+  }
 }</syntaxhighlight>
 {{out}}