Dmat=[0,0,3;4,5,6;7,8,9;10,10,10]

[C,M] = size(Dmat);
for i=1:M
    temp=0;
     for k=1:C
         if(Dmat(k,i)==0) temp=temp+1; end
     end
             Umat(k,i)=Dmat(k,i);
         else
             for r=1:C
                 if(r~=k)
                     Umat(r,i)=0.0;
                 else
                     Umat(r,i)=1;
        
                 end
             end
         end
     end
 end
  Umat