expected = { [1] };
actual = permutation_matrices(1);

unit_test_equals("permutation_matrices(1) works", ...
		 true, ...
		 isequal(actual,expected));
