function [G,H]=grextrat(x)
Ps=0.4;   % R$/kg-soluto
Px=0.01;  % R$/kg-solvente
F=1e4;    % kg-agua/h
xo=0.02;  % kg-soluto/kg-agua
s=7e-4;   % kg-solvente/kg-agua
%s=0;      % kg-solvente/kg-agua
m=4;      % kg-agua/kg-solvente

b=Px*F*xo/m;
c=Ps*F;
d=Px*F/m;

G(1,1)=-b./x(1).^2+d./x(2);
G(1,2)=c-d*x(1)./x(2).^2;
G(2,1)=-1;
G(2,2)=1;
H=[];
