function grad = gvalue(n,G_gamma, Ahandle, AThandle,mu2,gamma)
% Compute the gradient from the information 
% in the subroutine fvalue
U = G_gamma(1:n);
V = G_gamma(n+1:2*n);
grad = [U - gamma*AThandle(Ahandle(U)) + gamma*mu2*V; gamma*mu2*U + V];
end

