function g=maxentgrd(x)
global A;
global b;
global w;
global alpha;
global K;

xabs=abs(x);
g=2*A'*(A*x)-2*A'*b+alpha^2*(log(w.*xabs)+ones(size(w)))-K*2*((xabs-x)/2);
