function U=dr101(x2,x0)



%This function calculates 

%the dressed cnot (-1)^{n_{x2}(2) n_{x0}(0)},

%that acts on 3 bits called 2,1,0. 

%See related function dr11()



sig_0 = sig(x0);

sig_2 = sig(x2);



U=1/2*(eye(8) + kron(sig_2,eye(4)) + kron(eye(4),sig_0) - kron(sig_2,kron(eye(2),sig_0)));