Assistance in coding needed

Hello, I am trying to include the following simple code as one of my constraints in my optimization model:

if Xij = 1, then Xik=0 for all k not equal j.

Please assist. Cheers

Answer 1

[dpb]

One way; may be more better...presuming X is square...

is1=(X==1);
X(is1 & ~eye(size(X)))=0;

Oh!...easy enough if isn't square, too...

is1=(X==1);
X(is1 & ~eye(min(size(X))))=0;