function y = myobj(x) A = reshape(x(1:4),[2 2]); C = x(5:6); y = 0; for i=2 for j=1:2 y = y + C(i) * (A(j,1)*F([1 j]) - A(j,i)*F([i j])); end end return function r = F(x) r = length(find(x~=1)); return