foreach my $x (3 .. $rep ) {
    foreach my $y (2 .. $x-1 ) {
        foreach my $z (1 .. $y-1 ) {
            # expanded form of x*(x - a)*(x - b)*(x - c)
            # coeffs are in an array
            $haystack{"$x, $y, $z"} =  [
                poly_roots(
                    poly_derivative(
                        1, -$x - $y - $z, $x*$y + $x*$z + $y*$z, -$x*$y*$z, 0
                    )
                )
            ];
        }
    }
}