%
% f=maxentobj(x)
%
function f=maxentobj(x)
global A;
global b;
global w;
global alpha;
global K;
%
%
%
xabs=abs(x);
f=norm(A*x-b)^2+alpha^2*xabs'*log(w.*xabs)+K*norm((xabs-x)/2)^2;

