% U_CNOT   4x4 unitary matrix realizing a CNOT gate
function U=U_CNOT

U=[1 0 0 0;0 1 0 0;0 0 0 1;0 0 1 0];