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