function CON = intersect2(CONfeas,CONadd, varargin)        h1 = CONfeas;    h2 = CONadd;        H.A = [h1.CI; h2.CI];    H.B = [h1.DI; h2.DI];        linstart = size(H.B,1) + 1;        H.A = [H.A; h1.CE; h2.CE];    H.B = [H.B; h1.DE; h2.DE];    linend = size(H.B, 1);        if (linend >= linstart)        H.lin = linstart: linend;    end        h = cddmex('reduce_h', H);        CE = h.A(h.lin,:);    DE = h.B(h.lin,:);    CI = h.A;    DI = h.B;    CON = linearicon(CE,DE,CI,DI);