function y=g2max(g1,Pmax,A,N);[l,n]=size(g1);g2=zeros(1,n);for i=1:n   g21=Pmax/(N(2)+g1(i)*A(1,2)*(N(1)+A(2,1)*Pmax));   g22=(Pmax-g1(i)*N(1))/(g1(i)*A(2,1)*(N(2)+A(1,2)*Pmax));   g2(i)=min(g21,max(0,g22));end;y=g2;