A=[1 0; 0 1;0 -1] b=[-0.2;12;-8] x0=[-1;10] i=0; for w=logspace(-2,2,200) x=fmincon(sprintf('-abs((x(1)^2+x(2)^2)/(-(%f)^2-2*%f*sqrt(-1)*x(1)+x(1)^2+x(2)^2))',w,w),x0,A,b); m=abs((x(1)^2+x(2)^2)/(-w^2-2*x(1)*(w*sqrt(-1))+x(1)^2+x(2)^2)-1); i=i+1; W2(i,1)=w; W2(i,2)=m; end