%0-1 inter

clc,clear;

close all;

N = 10;

M=6;

A = rand(N);

A =(A*A’);

A = (A+A’)/2;

% x0 = [ones(M,1);zeros(N-M,1)];

rho=0.1;

cvx_begin

variable x(N,1)

variable x1(N,1)

minimize (x’*A*x+rho*(square_pos(norm(x.*(1-x1),2))+square_pos(norm(x-x1,2))))

subject to

for ii=1:N

x(ii)>=0;

x(ii)<=1;

x1(ii)>=0;

x1(ii)<=1;

end

sum(x)==M

cvx_end