function U=dr011(x1,x0) %This function calculates %the dressed cnot (-1)^{n_{x1}(1) n_{x0}(0)}, %that acts on 3 bits called 2,1,0. %See related function dr11(). U=kron(eye(2),dr11(x1,x0));