Pij = sdpvar (6, 1);Qij = sdpvar (6, 1);I = sdpvar (6, 1);V = sdpvar(6,1);for iofbranch = 1:6 Initialbus =Branch(iofbranch,2); Finalbus =Branch(iofbranch,3);
constraint = [ constraint, rcone( [Pij(iofbranch,1) ; Qij(iofbranch,1)], V(Initialbus,1), 1/2*I(iofbranch,1) ) ];
end