% Return standard single-qubit Pauli matrices I,X,Y,Z.

function [I,X,Y,Z] = pauli()
I = [1 0;0 1];
X = [0 1;1 0];
Y = [0 -1i;1i 0];
Z = [1 0;0 -1];
end