Inspired by Matt Parker’s video “Generating π from 1,000 random numbers” (https://youtu.be/RZBhSi_PwHU)

generate_pi.raku

unit sub MAIN(UInt :s(:$sides) = 120, UInt :r(:$rolls) = 1_000);

sub roll-die { (1..$sides).roll }

my $coprimes = 0;

for ^$rolls {
    my $a = roll-die();
    my $b = roll-die();

    $coprimes++ if $a gcd $b == 1;
}

say "π ≈ &sqrt(6 * $rolls / $coprimes)";