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)));

 