public function permutationMatrix() : Matrix
{
$size = count($this->permutation);
$permutationMatrix = [];
for ($i = 0; $i < $size; $i++) {
$permutationMatrix[] = [];
for ($j = 0; $j < $size; $j++) {
$permutationMatrix[$i][] = $this->permutation[$i] === $j ? 1 : 0;
}
}
return new Matrix($permutationMatrix);
}