function y=GG(P_tot)
global DD1 DD2 hh P N N_R;
I=eye(N_R);
y=0;
for jj=1:N,
    D1(1:N_R,1:N_R)=DD1(jj,:,:);
    D2(1:N_R,1:N_R)=DD2(jj,:,:);
    h(1:N_R,1)=hh(jj,:);
    y=y-log2(1+2*P(1,jj)*((P_tot(jj)/2)-P(1,jj))*((h'/(2*P(1,jj)*D1+I+((P_tot(jj)/2)-P(1,jj))*D2))*h));
end
	y=real(y);
return