sigma_x = [0 1;1 0];
sigma_y = [0 -1i;1i 0];
sigma_z = [1 0;0 -1];
sigma_i = [1 0;0 1];
id      = sigma_i;