function [gf_nlp,gg_nlp,C,B] = grad16(x,y)

if size(x,2) > 93
  y = x(94:128);         % NLP com integralidade relaxada
end

% * * * exemplo 16 - Daichendt e Grossmann(1994) - Trocadores de Calor * * *  

Fh=[15 3 8 9 5]';
Fc=[8 15 8 4 6]';
TINh=[500 460 440 350 350]';
TINc=[300 340 340 380 460]';
TOUTh=[340 400 400 310 320]';
TOUTc=[340 360 400 460 560]';
TINhu=580;
TINcu=300;
TOUThu=580;
TOUTcu=320;
CHU=125;
CCU=10;
EMAT=1;
C=300;
h=1.6;
U=(2/h)^(-1);
CF=900;
alfa=10^3;
n=5;


%  * * *  Variaveis binarias * * *

Y1=zeros(5);
Y2=zeros(5);
Y3=zeros(5);
Y4=zeros(5);
Y5=zeros(5);
Y6=zeros(5);
Y7=zeros(5);
Y1(1,5)=y(1);
Y2(1,4)=y(2);
Y3(1,4)=y(3);
Y3(2,4)=y(4);
Y4(1,3)=y(5);
Y4(1,4)=y(6);
Y4(2,3)=y(7);
Y4(2,4)=y(8);
Y4(3,3)=y(9);
Y4(3,4)=y(10);
Y5(1,2)=y(11);
Y5(1,3)=y(12);
Y5(2,2)=y(13);
Y5(3,2)=y(14);
Y6(1,2)=y(15);
Y6(1,3)=y(16);
Y6(4,2)=y(17);
Y6(4,3)=y(18);
Y6(5,2)=y(19);
Y6(5,3)=y(20);
Y7(1,1)=y(21);
Y7(2,1)=y(22);
Y7(3,1)=y(23);
Y7(4,1)=y(24);
Y7(5,1)=y(25);
Yh=y(26:30);
Yc=y(31:35);

%     * * * Variaveis continuas  * * *

DT1=zeros(5);
DT2=zeros(5);
DT3=zeros(5);
DT4=zeros(5);
DT5=zeros(5);
DT6=zeros(5);
DT7=zeros(5);
DT8=zeros(5);
DT1(1,5)=x(1);
DT2(1,5)=x(2);
DT2(1,4)=x(3);
DT3(1,4)=x(4);
DT3(2,4)=x(5);
DT4(1,4)=x(6);
DT4(2,4)=x(7);
DT4(1,3)=x(8);
DT4(2,3)=x(9);
DT4(3,3)=x(10);
DT4(3,4)=x(11);
DT5(1,3)=x(12);
DT5(1,4)=x(13);
DT5(2,3)=x(14);
DT5(2,4)=x(15);
DT5(3,3)=x(16);
DT5(3,4)=x(17);
DT5(1,2)=x(18);
DT5(2,2)=x(19);
DT5(3,2)=x(20);
DT6(1,2)=x(21);
DT6(1,3)=x(22);
DT6(2,2)=x(23);
DT6(3,2)=x(24);
DT6(4,2)=x(25);
DT6(4,3)=x(26);
DT6(5,2)=x(27);
DT6(5,3)=x(28);
DT7(1,2)=x(29);
DT7(1,3)=x(30);
DT7(4,2)=x(31);
DT7(4,3)=x(32);
DT7(5,2)=x(33);
DT7(5,3)=x(34);
DT7(1,1)=x(35);
DT7(2,1)=x(36);
DT7(3,1)=x(37);
DT7(4,1)=x(38);
DT7(5,1)=x(39);
DT8(1,1)=x(40);
DT8(2,1)=x(41);
DT8(3,1)=x(42);
DT8(4,1)=x(43);
DT8(5,1)=x(44);
DTh=x(45:49);
DTc=x(50:54);

Q1=zeros(5,4);
Q2=zeros(5,4);
Q3=zeros(5,4);
Q4=zeros(5,4);
Q5=zeros(5,4);
Q6=zeros(5,4);
Q7=zeros(5,4);
Q4(1,4)=x(55);
Q4(2,4)=x(56);
Q6(1,3)=x(57);
Q6(4,3)=x(58);

Th=zeros(5,8);
Th(:,1)=TINh;
Th(1,2)=x(59);
Th(1,3)=x(60);
Th(1,4)=x(61);
Th(1,5)=x(62);
Th(1,6)=x(63);
Th(1,7)=x(64);
Th(1,8)=x(65);
Th(2,3)=x(66);
Th(2,4)=x(67);
Th(2,5)=x(68);
Th(2,6)=x(69);
Th(2,7)=x(70);
Th(2,8)=x(71);
Th(3,4)=x(72);
Th(3,5)=x(73);
Th(3,6)=x(74);
Th(3,7)=x(75);
Th(3,8)=x(76);
Th(4,6)=x(77);
Th(4,7)=x(78);
Th(4,8)=x(79);
Th(5,6)=x(80);
Th(5,7)=x(81);
Th(5,8)=x(82);

Tc=zeros(5,8);
Tc(:,8)=TINc;
Tc(1,1)=x(83);
Tc(2,1)=x(84);
Tc(3,1)=x(85);
Tc(4,1)=x(86);
Tc(5,2)=x(87);
Tc(3,4)=x(88);
Tc(4,5)=x(89);
Tc(2,5)=x(90);
Tc(2,6)=x(91);
Tc(2,7)=x(92);
Tc(3,7)=x(93);

Q1(1,5)=Fh(1)*(Th(1,1)-Th(1,2));
Tc(5,1)=Q1(1,5)/Fc(5)+Tc(5,2);

Q6(1,2)=Fh(1)*(Th(1,6)-Th(1,7))-Q6(1,3);
Q6(4,2)=Fh(4)*(Th(4,6)-Th(4,7))-Q6(4,3);
Q6(5,2)=Fc(2)*(Tc(2,6)-Tc(2,7))-Q6(1,2)-Q6(4,2);
Q6(5,3)=Fh(5)*(Th(5,6)-Th(5,7))-Q6(5,2);
Tc(3,6)=(Q6(1,3)+Q6(4,3)+Q6(5,3))/Fc(3)+Tc(3,7);

Q5(2,2)=Fh(2)*(Th(2,5)-Th(2,6));
Q5(3,2)=Fh(3)*(Th(3,5)-Th(3,6));
Q5(1,2)=Fc(2)*(Tc(2,5)-Tc(2,6))-Q5(2,2)-Q5(3,2);
Q5(1,3)=Fh(1)*(Th(1,5)-Th(1,6))-Q5(1,2);
Tc(3,5)=Q5(1,3)/Fc(3)+Tc(3,6);

Q4(1,3)=Fh(1)*(Th(1,4)-Th(1,5))-Q4(1,4);
Q4(2,3)=Fh(2)*(Th(2,4)-Th(2,5))-Q4(2,4);
Q4(3,3)=Fc(3)*(Tc(3,4)-Tc(3,5))-Q4(1,3)-Q4(2,3);
Q4(3,4)=Fh(3)*(Th(3,4)-Th(3,5))-Q4(3,3);
Tc(4,4)=(Q4(1,4)+Q4(2,4)+Q4(3,4))/Fc(4)+Tc(4,5);

Q3(1,4)=Fh(1)*(Th(1,3)-Th(1,4));
Q3(2,4)=Fh(2)*(Th(2,3)-Th(2,4));
Tc(4,3)=(Q3(1,4)+Q3(2,4))/Fc(4)+Tc(4,4);

Q2(1,4)=Fh(1)*(Th(1,2)-Th(1,3));
Tc(4,2)=Q2(1,4)/Fc(4)+Tc(4,3);

Q7(1,1)=Fh(1)*(Th(1,7)-Th(1,8));
Q7(2,1)=Fh(2)*(Th(2,7)-Th(2,8));
Q7(3,1)=Fh(3)*(Th(3,7)-Th(3,8));
Q7(4,1)=Fh(4)*(Th(4,7)-Th(4,8));
Q7(5,1)=Fh(5)*(Th(5,7)-Th(5,8));
Tc(1,7)=(Q7(1,1)+Q7(2,1)+Q7(3,1)+Q7(4,1)+Q7(5,1))/Fc(1)+Tc(1,8);

Qh=Fc.*(TOUTc-Tc(:,1));
Qc=Fh.*(Th(:,8)-TOUTh);

omega=max(Fh.*(Th(:,1)-TOUTh),Fc.*(TOUTc-Tc(:,8)));

gf_nlp(1)=-(C/(3*U))*Q1(1,5)*((0.5*x(2)*(x(1)^2)+0.5*x(1)*(x(2)^2))^(-4/3))*...
(x(1)*x(2)+0.5*(x(2)^2));

gf_nlp(2)=-(C/(3*U))*Q1(1,5)*((0.5*x(2)*(x(1)^2)+0.5*x(1)*(x(2)^2))^(-4/3))*...
(x(1)*x(2)+0.5*(x(1)^2));

gf_nlp(3)=-(C/(3*U))*Q2(1,4)*((0.5*x(3)*(x(4)^2)+0.5*x(4)*(x(3)^2))^(-4/3))*...
(x(3)*x(4)+0.5*(x(4)^2));

gf_nlp(4)=-(C/(3*U))*Q2(1,4)*((0.5*x(3)*(x(4)^2)+0.5*x(4)*(x(3)^2))^(-4/3))*...
(x(3)*x(4)+0.5*(x(3)^2))-(C/(3*U))*Q3(1,4)*((0.5*x(4)*(x(6)^2)+0.5*x(6)*...
(x(4)^2))^(-4/3))*(x(4)*x(6)+0.5*(x(6)^2)); 
 
gf_nlp(5)=-(C/(3*U))*Q3(2,4)*((0.5*x(5)*(x(7)^2)+0.5*x(7)*(x(5)^2))^(-4/3))*...
(x(5)*x(7)+0.5*(x(7)^2)); 

gf_nlp(6)=-(C/(3*U))*Q3(1,4)*((0.5*x(4)*(x(6)^2)+0.5*x(6)*(x(4)^2))^(-4/3))*...
(x(4)*x(6)+0.5*(x(4)^2))-(C/(3*U))*x(55)*((0.5*x(6)*(x(13)^2)+0.5*x(13)*...
(x(6)^2))^(-4/3))*(x(6)*x(13)+0.5*(x(13)^2)); 

gf_nlp(7)=-(C/(3*U))*Q3(2,4)*((0.5*x(5)*(x(7)^2)+0.5*x(7)*(x(5)^2))^(-4/3))*...
(x(5)*x(7)+0.5*(x(5)^2))-(C/(3*U))*x(56)*((0.5*x(7)*(x(15)^2)+0.5*x(15)*...
(x(7)^2))^(-4/3))*(x(7)*x(15)+0.5*(x(15)^2));

gf_nlp(8)=-(C/(3*U))*Q4(1,3)*((0.5*x(8)*(x(12)^2)+0.5*x(12)*(x(8)^2))^(-4/3))*...
(x(8)*x(12)+0.5*(x(12)^2));

gf_nlp(9)=-(C/(3*U))*Q4(2,3)*((0.5*x(9)*(x(14)^2)+0.5*x(14)*(x(9)^2))^(-4/3))*...
(x(9)*x(14)+0.5*(x(14)^2));

gf_nlp(10)=-(C/(3*U))*Q4(3,3)*((0.5*x(10)*(x(16)^2)+0.5*x(16)*(x(10)^2))^(-4/3))...
*(x(10)*x(16)+0.5*(x(16)^2));

gf_nlp(11)=-(C/(3*U))*Q4(3,4)*((0.5*x(11)*(x(17)^2)+0.5*x(17)*(x(11)^2))^(-4/3))...
*(x(11)*x(17)+0.5*(x(17)^2));

gf_nlp(12)=-(C/(3*U))*Q4(1,3)*((0.5*x(8)*(x(12)^2)+0.5*x(12)*(x(8)^2))^(-4/3))...
*(x(8)*x(12)+0.5*(x(8)^2))-(C/(3*U))*Q5(1,3)*((0.5*x(12)*(x(22)^2)+0.5*x(22)*...
(x(12)^2))^(-4/3))*(x(12)*x(22)+0.5*(x(22)^2));

gf_nlp(13)=-(C/(3*U))*x(55)*((0.5*x(6)*(x(13)^2)+0.5*x(13)*(x(6)^2))^(-4/3))...
*(x(6)*x(13)+0.5*(x(6)^2));

gf_nlp(14)=-(C/(3*U))*Q4(2,3)*((0.5*x(9)*(x(14)^2)+0.5*x(14)*(x(9)^2))^(-4/3))...
*(x(9)*x(14)+0.5*(x(9)^2));

gf_nlp(15)=-(C/(3*U))*x(56)*((0.5*x(7)*(x(15)^2)+0.5*x(15)*(x(7)^2))^(-4/3))...
*(x(7)*x(15)+0.5*(x(7)^2));

gf_nlp(16)=-(C/(3*U))*Q4(3,3)*((0.5*x(10)*(x(16)^2)+0.5*x(16)*(x(10)^2))^(-4/3))...
*(x(10)*x(16)+0.5*(x(10)^2));

gf_nlp(17)=-(C/(3*U))*Q4(3,4)*((0.5*x(11)*(x(17)^2)+0.5*x(17)*(x(11)^2))^(-4/3))...
*(x(11)*x(17)+0.5*(x(11)^2));

gf_nlp(18)=-(C/(3*U))*Q5(1,2)*((0.5*x(18)*(x(21)^2)+0.5*x(21)*(x(18)^2))^(-4/3))...
*(x(18)*x(21)+0.5*(x(21)^2));

gf_nlp(19)=-(C/(3*U))*Q5(2,2)*((0.5*x(19)*(x(23)^2)+0.5*x(23)*(x(19)^2))^(-4/3))...
*(x(19)*x(23)+0.5*(x(23)^2));

gf_nlp(20)=-(C/(3*U))*Q5(3,2)*((0.5*x(20)*(x(24)^2)+0.5*x(24)*(x(20)^2))^(-4/3))...
*(x(20)*x(24)+0.5*(x(24)^2));

gf_nlp(21)=-(C/(3*U))*Q5(1,2)*((0.5*x(18)*(x(21)^2)+0.5*x(21)*(x(18)^2))^(-4/3))...
*(x(18)*x(21)+0.5*(x(18)^2))-(C/(3*U))*Q6(1,2)*((0.5*x(21)*(x(29)^2)+0.5*x(29)*...
(x(21)^2))^(-4/3))*(x(21)*x(29)+0.5*(x(29)^2));

gf_nlp(22)=-(C/(3*U))*Q5(1,3)*((0.5*x(12)*(x(22)^2)+0.5*x(22)*(x(12)^2))^(-4/3))...
*(x(12)*x(22)+0.5*(x(12)^2))-(C/(3*U))*x(57)*((0.5*x(22)*(x(30)^2)+0.5*x(30)*... 
(x(22)^2))^(-4/3))*(x(22)*x(30)+0.5*(x(30)^2));

gf_nlp(23)=-(C/(3*U))*Q5(2,2)*((0.5*x(19)*(x(23)^2)+0.5*x(23)*(x(19)^2))^(-4/3))...
*(x(19)*x(23)+0.5*(x(19)^2));

gf_nlp(24)=-(C/(3*U))*Q5(3,2)*((0.5*x(20)*(x(24)^2)+0.5*x(24)*(x(20)^2))^(-4/3))...
*(x(20)*x(24)+0.5*(x(20)^2));

gf_nlp(25)=-(C/(3*U))*Q6(4,2)*((0.5*x(25)*(x(31)^2)+0.5*x(31)*(x(25)^2))^(-4/3))...
*(x(25)*x(31)+0.5*(x(31)^2));

gf_nlp(26)=-(C/(3*U))*x(58)*((0.5*x(26)*(x(32)^2)+0.5*x(32)*(x(26)^2))^(-4/3))...
*(x(26)*x(32)+0.5*(x(32)^2));

gf_nlp(27)=-(C/(3*U))*Q6(5,2)*((0.5*x(27)*(x(33)^2)+0.5*x(33)*(x(27)^2))^(-4/3))...
*(x(27)*x(33)+0.5*(x(33)^2));

gf_nlp(28)=-(C/(3*U))*Q6(5,3)*((0.5*x(28)*(x(34)^2)+0.5*x(34)*(x(28)^2))^(-4/3))...
*(x(28)*x(34)+0.5*(x(34)^2));

gf_nlp(29)=-(C/(3*U))*Q6(1,2)*((0.5*x(21)*(x(29)^2)+0.5*x(29)*(x(21)^2))^(-4/3))...
*(x(21)*x(29)+0.5*(x(21)^2));

gf_nlp(30)=-(C/(3*U))*x(57)*((0.5*x(22)*(x(30)^2)+0.5*x(30)*(x(22)^2))^(-4/3))...
*(x(22)*x(30)+0.5*(x(22)^2));

gf_nlp(31)=-(C/(3*U))*Q6(4,2)*((0.5*x(25)*(x(31)^2)+0.5*x(31)*(x(25)^2))^(-4/3))...
*(x(25)*x(31)+0.5*(x(25)^2));

gf_nlp(32)=-(C/(3*U))*x(58)*((0.5*x(26)*(x(32)^2)+0.5*x(32)*(x(26)^2))^(-4/3))...
*(x(26)*x(32)+0.5*(x(26)^2));

gf_nlp(33)=-(C/(3*U))*Q6(5,2)*((0.5*x(27)*(x(33)^2)+0.5*x(33)*(x(27)^2))^(-4/3))...
*(x(27)*x(33)+0.5*(x(27)^2));

gf_nlp(34)=-(C/(3*U))*Q6(5,3)*((0.5*x(28)*(x(34)^2)+0.5*x(34)*(x(28)^2))^(-4/3))...
*(x(28)*x(34)+0.5*(x(28)^2));

gf_nlp(35)=-(C/(3*U))*Q7(1,1)*((0.5*x(35)*(x(40)^2)+0.5*x(40)*(x(35)^2))^(-4/3))...
*(x(35)*x(40)+0.5*(x(40)^2));

gf_nlp(36)=-(C/(3*U))*Q7(2,1)*((0.5*x(36)*(x(41)^2)+0.5*x(41)*(x(36)^2))^(-4/3))...
*(x(36)*x(41)+0.5*(x(41)^2));

gf_nlp(37)=-(C/(3*U))*Q7(3,1)*((0.5*x(37)*(x(42)^2)+0.5*x(42)*(x(37)^2))^(-4/3))...
*(x(37)*x(42)+0.5*(x(42)^2));

gf_nlp(38)=-(C/(3*U))*Q7(4,1)*((0.5*x(38)*(x(43)^2)+0.5*x(43)*(x(38)^2))^(-4/3))...
*(x(38)*x(43)+0.5*(x(43)^2));

gf_nlp(39)=-(C/(3*U))*Q7(5,1)*((0.5*x(39)*(x(44)^2)+0.5*x(44)*(x(39)^2))^(-4/3))...
*(x(39)*x(44)+0.5*(x(44)^2));

gf_nlp(40)=-(C/(3*U))*Q7(1,1)*((0.5*x(35)*(x(40)^2)+0.5*x(40)*(x(35)^2))^(-4/3))...
*(x(35)*x(40)+0.5*(x(35)^2));

gf_nlp(41)=-(C/(3*U))*Q7(2,1)*((0.5*x(36)*(x(41)^2)+0.5*x(41)*(x(36)^2))^(-4/3))...
*(x(36)*x(41)+0.5*(x(36)^2));

gf_nlp(42)=-(C/(3*U))*Q7(3,1)*((0.5*x(37)*(x(42)^2)+0.5*x(42)*(x(37)^2))^(-4/3))...
*(x(37)*x(42)+0.5*(x(37)^2));

gf_nlp(43)=-(C/(3*U))*Q7(4,1)*((0.5*x(38)*(x(43)^2)+0.5*x(43)*(x(38)^2))^(-4/3))...
*(x(38)*x(43)+0.5*(x(38)^2));

gf_nlp(44)=-(C/(3*U))*Q7(5,1)*((0.5*x(39)*(x(44)^2)+0.5*x(44)*(x(39)^2))^(-4/3))...
*(x(39)*x(44)+0.5*(x(39)^2));

gf_nlp(45)=-(C/(3*U))*Qh(1)*((0.5*(TINc(1)-TOUThu)*(x(45)^2+TINc(1)*x(45)-...
TOUThu*x(45)))^(-4/3))*(0.5*(TINc(1)-TOUThu)*(2*x(45)+TINc(1)-TOUThu));

gf_nlp(46)=-(C/(3*U))*Qh(2)*((0.5*(TINc(2)-TOUThu)*(x(46)^2+TINc(2)*x(46)-...
TOUThu*x(46)))^(-4/3))*(0.5*(TINc(2)-TOUThu)*(2*x(46)+TINc(2)-TOUThu));

gf_nlp(47)=-(C/(3*U))*Qh(3)*((0.5*(TINc(3)-TOUThu)*(x(47)^2+TINc(3)*x(47)-...
TOUThu*x(47)))^(-4/3))*(0.5*(TINc(3)-TOUThu)*(2*x(47)+TINc(3)-TOUThu));

gf_nlp(48)=-(C/(3*U))*Qh(4)*((0.5*(TINc(4)-TOUThu)*(x(48)^2+TINc(4)*x(48)-...
TOUThu*x(48)))^(-4/3))*(0.5*(TINc(4)-TOUThu)*(2*x(48)+TINc(4)-TOUThu));

gf_nlp(49)=-(C/(3*U))*Qh(5)*((0.5*(TINc(5)-TOUThu)*(x(49)^2+TINc(5)*x(49)-...
TOUThu*x(49)))^(-4/3))*(0.5*(TINc(5)-TOUThu)*(2*x(49)+TINc(5)-TOUThu));

gf_nlp(50)=-(C/(3*U))*Qc(1)*((0.5*(TOUTh(1)-TINcu)*(x(50)^2+TOUTh(1)*x(50)-...
TINcu*x(50)))^(-4/3))*(0.5*(TOUTh(1)-TINcu)*(2*x(50)+TOUTh(1)-TINcu));

gf_nlp(51)=-(C/(3*U))*Qc(2)*((0.5*(TOUTh(2)-TINcu)*(x(51)^2+TOUTh(2)*x(51)-...
TINcu*x(51)))^(-4/3))*(0.5*(TOUTh(2)-TINcu)*(2*x(51)+TOUTh(2)-TINcu));

gf_nlp(52)=-(C/(3*U))*Qc(3)*((0.5*(TOUTh(3)-TINcu)*(x(52)^2+TOUTh(3)*x(52)-...
TINcu*x(52)))^(-4/3))*(0.5*(TOUTh(3)-TINcu)*(2*x(52)+TOUTh(3)-TINcu));

gf_nlp(53)=-(C/(3*U))*Qc(4)*((0.5*(TOUTh(4)-TINcu)*(x(53)^2+TOUTh(4)*x(53)-...
TINcu*x(53)))^(-4/3))*(0.5*(TOUTh(4)-TINcu)*(2*x(53)+TOUTh(4)-TINcu));

gf_nlp(54)=-(C/(3*U))*Qc(5)*((0.5*(TOUTh(5)-TINcu)*(x(54)^2+TOUTh(5)*x(54)-...
TINcu*x(54)))^(-4/3))*(0.5*(TOUTh(5)-TINcu)*(2*x(54)+TOUTh(5)-TINcu));

gf_nlp(55)=-(C/U)*((0.5*(x(8)^2)*x(12)+0.5*(x(12)^2)*x(8))^(-1/3))+...
(C/U)*((0.5*(x(6)^2)*x(13)+0.5*(x(13)^2)*x(6))^(-1/3))+...
(C/U)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3));

gf_nlp(56)=-(C/U)*((0.5*(x(9)^2)*x(14)+0.5*(x(14)^2)*x(9))^(-1/3))+...
(C/U)*((0.5*(x(7)^2)*x(15)+0.5*(x(15)^2)*x(7))^(-1/3))+...
(C/U)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3)); 

gf_nlp(57)=-(C/U)*((0.5*(x(21)^2)*x(29)+0.5*(x(29)^2)*x(21))^(-1/3))+...
(C/U)*((0.5*(x(22)^2)*x(30)+0.5*(x(30)^2)*x(22))^(-1/3))+...
(C/U)*((0.5*(x(27)^2)*x(33)+0.5*(x(33)^2)*x(27))^(-1/3))-...
(C/U)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3)); 

gf_nlp(58)=-(C/U)*((0.5*(x(25)^2)*x(31)+0.5*(x(31)^2)*x(25))^(-1/3))+...
(C/U)*((0.5*(x(26)^2)*x(32)+0.5*(x(32)^2)*x(26))^(-1/3))+...
(C/U)*((0.5*(x(27)^2)*x(33)+0.5*(x(33)^2)*x(27))^(-1/3))-...
(C/U)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3)); 

gf_nlp(59)=-(C/U)*Fh(1)*((0.5*(x(1)^2)*x(2)+0.5*(x(2)^2)*x(1))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(3)^2)*x(4)+0.5*(x(4)^2)*x(3))^(-1/3))+CHU*Fh(1)+...
(C/U)*Fh(1)*((0.5*(TINc(5)-TOUThu)*(x(49)^2+TINc(5)*x(49)-TOUThu*x(49)))^(-1/3));

gf_nlp(60)=-(C/U)*Fh(1)*((0.5*(x(3)^2)*x(4)+0.5*(x(4)^2)*x(3))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(4)^2)*x(6)+0.5*(x(6)^2)*x(4))^(-1/3));

gf_nlp(61)=-(C/U)*Fh(1)*((0.5*(x(4)^2)*x(6)+0.5*(x(6)^2)*x(4))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(8)^2)*x(12)+0.5*(x(12)^2)*x(8))^(-1/3))-...
(C/U)*Fh(1)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3));

gf_nlp(62)=-(C/U)*Fh(1)*((0.5*(x(8)^2)*x(12)+0.5*(x(12)^2)*x(8))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(12)^2)*x(22)+0.5*(x(22)^2)*x(12))^(-1/3));

gf_nlp(63)=-(C/U)*Fh(1)*((0.5*(x(12)^2)*x(22)+0.5*(x(22)^2)*x(12))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(21)^2)*x(29)+0.5*(x(29)^2)*x(21))^(-1/3))-...
(C/U)*Fh(1)*((0.5*(x(27)^2)*x(33)+0.5*(x(33)^2)*x(27))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3));

gf_nlp(64)=(C/U)*Fh(1)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*Fh(1)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))-...
(C/U)*Fh(1)*((0.5*(x(21)^2)*x(29)+0.5*(x(29)^2)*x(21))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(27)^2)*x(33)+0.5*(x(33)^2)*x(27))^(-1/3))-...
(C/U)*Fh(1)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3))+...
(C/U)*Fh(1)*((0.5*(x(35)^2)*x(40)+0.5*(x(40)^2)*x(35))^(-1/3));

gf_nlp(65)=-(C/U)*Fh(1)*((0.5*(x(35)^2)*x(40)+0.5*(x(40)^2)*x(35))^(-1/3))+...
CCU*Fh(1)+(C/U)*Fh(1)*((0.5*(TOUTh(1)-TINcu)*(x(50)^2+TOUTh(1)*x(50)-TINcu*...
x(50)))^(-1/3));

gf_nlp(66)=(C/U)*Fh(2)*((0.5*(x(5)^2)*x(7)+0.5*(x(7)^2)*x(5))^(-1/3));

gf_nlp(67)=-(C/U)*Fh(2)*((0.5*(x(5)^2)*x(7)+0.5*(x(7)^2)*x(5))^(-1/3))+...
(C/U)*Fh(2)*((0.5*(x(9)^2)*x(14)+0.5*(x(14)^2)*x(9))^(-1/3))-...
(C/U)*Fh(2)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))+...
(C/U)*Fh(2)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3));

gf_nlp(68)=-(C/U)*Fh(2)*((0.5*(x(9)^2)*x(14)+0.5*(x(14)^2)*x(9))^(-1/3))-...
(C/U)*Fh(2)*((0.5*(x(18)^2)*x(21)+0.5*(x(21)^2)*x(18))^(-1/3))+...
(C/U)*Fh(2)*((0.5*(x(12)^2)*x(22)+0.5*(x(22)^2)*x(12))^(-1/3))+...
(C/U)*Fh(2)*((0.5*(x(19)^2)*x(23)+0.5*(x(23)^2)*x(19))^(-1/3));

gf_nlp(69)=(C/U)*Fh(2)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*Fh(2)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))+...
(C/U)*Fh(2)*((0.5*(x(18)^2)*x(21)+0.5*(x(21)^2)*x(18))^(-1/3))-...
(C/U)*Fh(2)*((0.5*(x(12)^2)*x(22)+0.5*(x(22)^2)*x(12))^(-1/3))-...
(C/U)*Fh(2)*((0.5*(x(19)^2)*x(23)+0.5*(x(23)^2)*x(19))^(-1/3));

gf_nlp(70)=(C/U)*Fh(2)*((0.5*(x(36)^2)*x(41)+0.5*(x(41)^2)*x(36))^(-1/3));

gf_nlp(71)=-(C/U)*Fh(2)*((0.5*(x(36)^2)*x(41)+0.5*(x(41)^2)*x(36))^(-1/3))+...
CCU*Fh(2)+(C/U)*Fh(2)*((0.5*(TOUTh(2)-TINcu)*(x(51)^2+TOUTh(2)*x(51)-TINcu*...
x(51)))^(-1/3));

gf_nlp(72)=(C/U)*Fh(3)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3));

gf_nlp(73)=-(C/U)*Fh(3)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*Fh(3)*((0.5*(x(18)^2)*x(21)+0.5*(x(21)^2)*x(18))^(-1/3))+...
(C/U)*Fh(3)*((0.5*(x(12)^2)*x(22)+0.5*(x(22)^2)*x(12))^(-1/3))+...
(C/U)*Fh(3)*((0.5*(x(20)^2)*x(24)+0.5*(x(24)^2)*x(20))^(-1/3));

gf_nlp(74)=(C/U)*Fh(3)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*Fh(3)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))+...
(C/U)*Fh(3)*((0.5*(x(18)^2)*x(21)+0.5*(x(21)^2)*x(18))^(-1/3))-...
(C/U)*Fh(3)*((0.5*(x(12)^2)*x(22)+0.5*(x(22)^2)*x(12))^(-1/3))-...
(C/U)*Fh(3)*((0.5*(x(20)^2)*x(24)+0.5*(x(24)^2)*x(20))^(-1/3));

gf_nlp(75)=(C/U)*Fh(3)*((0.5*(x(37)^2)*x(42)+0.5*(x(42)^2)*x(37))^(-1/3));

gf_nlp(76)=-(C/U)*Fh(3)*((0.5*(x(37)^2)*x(42)+0.5*(x(42)^2)*x(37))^(-1/3))+...
CCU*Fh(3)+(C/U)*Fh(3)*((0.5*(TOUTh(3)-TINcu)*(x(52)^2+TOUTh(3)*x(52)-TINcu*...
x(52)))^(-1/3));

gf_nlp(77)=-(C/U)*Fh(4)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))+...
(C/U)*Fh(4)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))+...
(C/U)*Fh(4)*((0.5*(x(25)^2)*x(31)+0.5*(x(31)^2)*x(25))^(-1/3))-...
(C/U)*Fh(4)*((0.5*(x(27)^2)*x(33)+0.5*(x(33)^2)*x(27))^(-1/3))+...
(C/U)*Fh(4)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3));

gf_nlp(78)=(C/U)*Fh(4)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*Fh(4)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))-...
(C/U)*Fh(4)*((0.5*(x(25)^2)*x(31)+0.5*(x(31)^2)*x(25))^(-1/3))+...
(C/U)*Fh(4)*((0.5*(x(27)^2)*x(33)+0.5*(x(33)^2)*x(27))^(-1/3))-...
(C/U)*Fh(4)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3))+...
(C/U)*Fh(4)*((0.5*(x(38)^2)*x(43)+0.5*(x(43)^2)*x(38))^(-1/3));

gf_nlp(79)=-(C/U)*Fh(4)*((0.5*(x(38)^2)*x(43)+0.5*(x(43)^2)*x(38))^(-1/3))+...
CCU*Fh(4)+(C/U)*Fh(4)*((0.5*(TOUTh(4)-TINcu)*(x(53)^2+TOUTh(4)*x(53)-TINcu*...
x(53)))^(-1/3));

gf_nlp(80)=-(C/U)*Fh(5)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))+...
(C/U)*Fh(5)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))+...
(C/U)*Fh(5)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3));

gf_nlp(81)=(C/U)*Fh(5)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*Fh(5)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))-...
(C/U)*Fh(5)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3))+...
(C/U)*Fh(5)*((0.5*(x(39)^2)*x(44)+0.5*(x(44)^2)*x(39))^(-1/3));

gf_nlp(82)=-(C/U)*Fh(5)*((0.5*(x(39)^2)*x(44)+0.5*(x(44)^2)*x(39))^(-1/3))+...
CCU*Fh(5)+(C/U)*Fh(5)*((0.5*(TOUTh(5)-TINcu)*(x(54)^2+TOUTh(5)*x(54)-TINcu*...
x(54)))^(-1/3));

gf_nlp(83)=-CHU*Fc(1)-(C/U)*Fc(1)*((0.5*(TINc(1)-TOUThu)*(x(45)^2+TINc(1)*...
x(45)-TOUThu*x(45)))^(-1/3));

gf_nlp(84)=-CHU*Fc(2)-(C/U)*Fc(2)*((0.5*(TINc(2)-TOUThu)*(x(46)^2+TINc(2)*...
x(46)-TOUThu*x(46)))^(-1/3));

gf_nlp(85)=-CHU*Fc(3)-(C/U)*Fc(3)*((0.5*(TINc(3)-TOUThu)*(x(47)^2+TINc(3)*...
x(47)-TOUThu*x(47)))^(-1/3));

gf_nlp(86)=-CHU*Fc(4)-(C/U)*Fc(4)*((0.5*(TINc(4)-TOUThu)*(x(48)^2+TINc(4)*...
x(48)-TOUThu*x(48)))^(-1/3));

gf_nlp(87)=-CHU*Fc(5)-(C/U)*Fc(5)*((0.5*(TINc(5)-TOUThu)*(x(49)^2+TINc(5)*...
x(49)-TOUThu*x(49)))^(-1/3));

gf_nlp(88)=(C/U)*Fc(3)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*Fc(3)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3));

gf_nlp(89)=0;

gf_nlp(90)=(C/U)*Fc(2)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))-...
(C/U)*Fc(2)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))+...
(C/U)*Fc(2)*((0.5*(x(18)^2)*x(21)+0.5*(x(21)^2)*x(18))^(-1/3))-...
(C/U)*Fc(2)*((0.5*(x(12)^2)*x(22)+0.5*(x(22)^2)*x(12))^(-1/3));

gf_nlp(91)=-(C/U)*Fc(2)*((0.5*(x(18)^2)*x(21)+0.5*(x(21)^2)*x(18))^(-1/3))+...
(C/U)*Fc(2)*((0.5*(x(12)^2)*x(22)+0.5*(x(22)^2)*x(12))^(-1/3))+...
(C/U)*Fc(2)*((0.5*(x(27)^2)*x(33)+0.5*(x(33)^2)*x(27))^(-1/3))-...
(C/U)*Fc(2)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3));

gf_nlp(92)=-(C/U)*Fc(2)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))+...
(C/U)*Fc(2)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3))-...
(C/U)*Fc(2)*((0.5*(x(27)^2)*x(33)+0.5*(x(33)^2)*x(27))^(-1/3))+...
(C/U)*Fc(2)*((0.5*(x(28)^2)*x(34)+0.5*(x(34)^2)*x(28))^(-1/3));

gf_nlp(93)=-(C/U)*Fc(3)*((0.5*(x(10)^2)*x(16)+0.5*(x(16)^2)*x(10))^(-1/3))+...
(C/U)*Fc(3)*((0.5*(x(11)^2)*x(17)+0.5*(x(17)^2)*x(11))^(-1/3));

gg_nlp=zeros(158,93);
gg_nlp(1,59)= 1;
gg_nlp(2,59)= -1;
gg_nlp(2,60)= 1;
gg_nlp(3,60)= -1;
gg_nlp(3,61)=1;
gg_nlp(4,66)=-1;
gg_nlp(4,67)=1;
gg_nlp(5,61)=-1;
gg_nlp(5,62)=1;
gg_nlp(6,67)=-1;
gg_nlp(6,68)=1;
gg_nlp(7,72)=-1;
gg_nlp(7,73)=1;
gg_nlp(8,62)=-1;
gg_nlp(8,63)=1;
gg_nlp(9,68)=-1;
gg_nlp(9,69)=1;
gg_nlp(10,73)=-1;
gg_nlp(10,74)=1;
gg_nlp(11,63)=-1;
gg_nlp(11,64)=1;
gg_nlp(12,77)=-1;
gg_nlp(12,78)=1;
gg_nlp(13,80)=-1;
gg_nlp(13,81)=1;
gg_nlp(14,64)=-1;
gg_nlp(14,65)=1;
gg_nlp(15,70)=-1;
gg_nlp(15,71)=1;
gg_nlp(16,75)=-1;
gg_nlp(16,76)=1;
gg_nlp(17,78)=-1;
gg_nlp(17,79)=1;
gg_nlp(18,81)=-1;
gg_nlp(18,82)=1;
gg_nlp(19,59)=Fh(1)/Fc(5);
gg_nlp(20,59)=-Fh(1)/Fc(4);
gg_nlp(20,60)=Fh(1)/Fc(4);
gg_nlp(21,60)=-Fh(1)/Fc(4);
gg_nlp(21,61)=Fh(1)/Fc(4);
gg_nlp(21,66)=-Fh(2)/Fc(4);
gg_nlp(21,67)=Fh(2)/Fc(4);
gg_nlp(22,62)=Fh(1)/Fc(3);
gg_nlp(22,64)=-Fh(1)/Fc(3);
gg_nlp(22,68)=Fh(2)/Fc(3);
gg_nlp(22,69)=-Fh(2)/Fc(3);
gg_nlp(22,73)=Fh(3)/Fc(3);
gg_nlp(22,74)=-Fh(3)/Fc(3);
gg_nlp(22,77)=Fh(4)/Fc(3);
gg_nlp(22,78)=-Fh(4)/Fc(3);
gg_nlp(22,80)=Fh(5)/Fc(3);
gg_nlp(22,81)=-Fh(5)/Fc(3);
gg_nlp(22,88)=-1;
gg_nlp(22,90)=-Fc(2)/Fc(3);
gg_nlp(22,92)=Fc(2)/Fc(3);
gg_nlp(22,93)=1;
gg_nlp(23,61)=-Fh(1)/Fc(4);
gg_nlp(23,64)=Fh(1)/Fc(4);
gg_nlp(23,67)=-Fh(2)/Fc(4);
gg_nlp(23,69)=Fh(2)/Fc(4);
gg_nlp(23,72)=-Fh(3)/Fc(4);
gg_nlp(23,74)=Fh(3)/Fc(4);
gg_nlp(23,77)=-Fh(4)/Fc(4);
gg_nlp(23,78)=Fh(4)/Fc(4);
gg_nlp(23,80)=-Fh(5)/Fc(4);
gg_nlp(23,81)=Fh(5)/Fc(4);
gg_nlp(23,88)=Fc(3)/Fc(4);
gg_nlp(23,90)=Fc(2)/Fc(4);
gg_nlp(23,92)=-Fc(2)/Fc(4);
gg_nlp(23,93)=-Fc(3)/Fc(4);
gg_nlp(24,90)=-1;
gg_nlp(24,91)=1;
gg_nlp(25,62)=-Fh(1)/Fc(3);
gg_nlp(25,63)=Fh(1)/Fc(3);
gg_nlp(25,68)=-Fh(2)/Fc(3);
gg_nlp(25,69)=Fh(2)/Fc(3);
gg_nlp(25,73)=-Fh(3)/Fc(3);
gg_nlp(25,74)=Fh(3)/Fc(3);
gg_nlp(25,90)=Fc(2)/Fc(3);
gg_nlp(25,91)=-Fc(2)/Fc(3);
gg_nlp(26,91)=-1;
gg_nlp(26,92)=1;
gg_nlp(27,63)=-Fh(1)/Fc(3);
gg_nlp(27,64)=Fh(1)/Fc(3);
gg_nlp(27,77)=-Fh(4)/Fc(3);
gg_nlp(27,78)=Fh(4)/Fc(3);
gg_nlp(27,80)=-Fh(5)/Fc(3);
gg_nlp(27,81)=Fh(5)/Fc(3);
gg_nlp(27,91)=Fc(2)/Fc(3);
gg_nlp(27,92)=-Fc(2)/Fc(3);
gg_nlp(28,64)=-Fh(1)/Fc(1);
gg_nlp(28,65)=Fh(1)/Fc(1);
gg_nlp(28,70)=-Fh(2)/Fc(1);
gg_nlp(28,71)=Fh(2)/Fc(1);
gg_nlp(28,75)=-Fh(3)/Fc(1);
gg_nlp(28,76)=Fh(3)/Fc(1);
gg_nlp(28,78)=-Fh(4)/Fc(1);
gg_nlp(28,79)=Fh(4)/Fc(1);
gg_nlp(28,81)=-Fh(5)/Fc(1);
gg_nlp(28,82)=Fh(5)/Fc(1);
gg_nlp(29,65)=-1;
gg_nlp(30,71)=-1;
gg_nlp(31,76)=-1;
gg_nlp(32,79)=-1;
gg_nlp(33,82)=-1;
gg_nlp(34,83)=1;
gg_nlp(35,84)=1;
gg_nlp(36,85)=1;
gg_nlp(37,86)=1;
gg_nlp(38,87)=1;
gg_nlp(38,59)=-Fh(1)/Fc(5);
gg_nlp(39,59)=-Fh(1);
gg_nlp(40,59)=Fh(1);
gg_nlp(40,60)=-Fh(1);
gg_nlp(41,60)=Fh(1);
gg_nlp(41,61)=-Fh(1);
gg_nlp(42,66)=Fh(2);
gg_nlp(42,67)=-Fh(2);
gg_nlp(43,55)=-1;
gg_nlp(43,61)=Fh(1);
gg_nlp(43,62)=-Fh(1);
gg_nlp(44,55)=1;
gg_nlp(45,56)=-1;
gg_nlp(45,67)=Fh(2);
gg_nlp(45,68)=-Fh(2);
gg_nlp(46,56)=1;
gg_nlp(47,55)=1;
gg_nlp(47,56)=1;
gg_nlp(47,61)=-Fh(1);
gg_nlp(47,64)=Fh(1);
gg_nlp(47,67)=-Fh(2);
gg_nlp(47,69)=Fh(2);
gg_nlp(47,73)=-Fh(3);
gg_nlp(47,74)=Fh(3);
gg_nlp(47,77)=-Fh(4);
gg_nlp(47,78)=Fh(4);
gg_nlp(47,80)=-Fh(5);
gg_nlp(47,81)=Fh(5);
gg_nlp(47,88)=Fc(3);
gg_nlp(47,90)=Fc(2);
gg_nlp(47,92)=-Fc(2);
gg_nlp(47,93)=-Fc(3);
gg_nlp(48,55)=-1;
gg_nlp(48,56)=-1;
gg_nlp(48,61)=Fh(1);
gg_nlp(48,64)=-Fh(1);
gg_nlp(48,67)=Fh(2);
gg_nlp(48,69)=-Fh(2);
gg_nlp(48,72)=Fh(3);
gg_nlp(48,74)=-Fh(3);
gg_nlp(48,77)=Fh(4);
gg_nlp(48,78)=-Fh(4);
gg_nlp(48,80)=Fh(5);
gg_nlp(48,81)=-Fh(5);
gg_nlp(48,88)=-Fc(3);
gg_nlp(48,90)=-Fc(2);
gg_nlp(48,92)=Fc(2);
gg_nlp(48,93)=Fc(3);
gg_nlp(49,68)=-Fh(2);
gg_nlp(49,69)=Fh(2);
gg_nlp(49,73)=-Fh(3);
gg_nlp(49,74)=Fh(3);
gg_nlp(49,90)=Fc(2);
gg_nlp(49,91)=-Fc(2);
gg_nlp(50,62)=Fh(1);
gg_nlp(50,63)=-Fh(1);
gg_nlp(50,68)=Fh(2);
gg_nlp(50,69)=-Fh(2);
gg_nlp(50,73)=Fh(3);
gg_nlp(50,74)=-Fh(3);
gg_nlp(50,90)=-Fc(2);
gg_nlp(50,91)=Fc(2);
gg_nlp(51,68)=Fh(2);
gg_nlp(51,69)=-Fh(2);
gg_nlp(52,73)=Fh(3);
gg_nlp(52,74)=-Fh(3);
gg_nlp(53,57)=-1;
gg_nlp(53,63)=Fh(1);
gg_nlp(53,64)=-Fh(1);
gg_nlp(54,57)=1;
gg_nlp(55,58)=-1;
gg_nlp(55,77)=Fh(4);
gg_nlp(55,78)=-Fh(4);
gg_nlp(56,58)=1;
gg_nlp(57,57)=1;
gg_nlp(57,58)=1;
gg_nlp(57,63)=-Fh(1);
gg_nlp(57,64)=Fh(1);
gg_nlp(57,77)=-Fh(4);
gg_nlp(57,78)=Fh(4);
gg_nlp(57,91)=Fc(2);
gg_nlp(57,92)=-Fc(2);
gg_nlp(58,57)=-1;
gg_nlp(58,58)=-1;
gg_nlp(58,63)=Fh(1);
gg_nlp(58,64)=-Fh(1);
gg_nlp(58,77)=Fh(4);
gg_nlp(58,78)=-Fh(4);
gg_nlp(58,80)=Fh(5);
gg_nlp(58,81)=-Fh(5);
gg_nlp(58,91)=-Fc(2);
gg_nlp(58,92)=Fc(2);
gg_nlp(59,64)=Fh(1);
gg_nlp(59,65)=-Fh(1);
gg_nlp(60,70)=Fh(2);
gg_nlp(60,71)=-Fh(2);
gg_nlp(61,75)=Fh(3);
gg_nlp(61,76)=-Fh(3);
gg_nlp(62,78)=Fh(4);
gg_nlp(62,79)=-Fh(4);
gg_nlp(63,81)=Fh(5);
gg_nlp(63,82)=-Fh(5);
gg_nlp(64,65)=Fh(1);
gg_nlp(65,71)=Fh(2);
gg_nlp(66,76)=Fh(3);
gg_nlp(67,79)=Fh(4);
gg_nlp(68,82)=Fh(5);
gg_nlp(69,83)=-Fc(1);
gg_nlp(70,84)=-Fc(2);
gg_nlp(71,85)=-Fc(3);
gg_nlp(72,86)=-Fc(4);
gg_nlp(73,59)=Fh(1);
gg_nlp(73,87)=-Fc(5);
gg_nlp(74,1)=1;
gg_nlp(74,59)=-Fh(1)/Fc(5);
gg_nlp(74,87)=1;
gg_nlp(75,2)=1;
gg_nlp(75,59)=-1;
gg_nlp(75,87)=1;
gg_nlp(76,3)=1;
gg_nlp(76,59)=-1+Fh(1)/Fc(4);
gg_nlp(76,64)=-Fh(1)/Fc(4);
gg_nlp(76,66)=Fh(2)/Fc(4);
gg_nlp(76,69)=-Fh(2)/Fc(4);
gg_nlp(76,72)=Fh(3)/Fc(4);
gg_nlp(76,74)=-Fh(3)/Fc(4);
gg_nlp(76,77)=Fh(4)/Fc(4);
gg_nlp(76,78)=-Fh(4)/Fc(4);
gg_nlp(76,80)=Fh(5)/Fc(4);
gg_nlp(76,81)=-Fh(5)/Fc(4);
gg_nlp(76,88)=-Fc(3)/Fc(4);
gg_nlp(76,89)=1;
gg_nlp(76,90)=-Fc(2)/Fc(4);
gg_nlp(76,92)=Fc(2)/Fc(4);
gg_nlp(76,93)=Fc(3)/Fc(4);
gg_nlp(77,4)=1;
gg_nlp(77,60)=-1+Fh(1)/Fc(4);
gg_nlp(77,64)=-Fh(1)/Fc(4);
gg_nlp(77,66)=Fh(2)/Fc(4);
gg_nlp(77,69)=-Fh(2)/Fc(4);
gg_nlp(77,72)=Fh(3)/Fc(4);
gg_nlp(77,74)=-Fh(3)/Fc(4);
gg_nlp(77,77)=Fh(4)/Fc(4);
gg_nlp(77,78)=-Fh(4)/Fc(4);
gg_nlp(77,80)=Fh(5)/Fc(4);
gg_nlp(77,81)=-Fh(5)/Fc(4);
gg_nlp(77,88)=-Fc(3)/Fc(4);
gg_nlp(77,89)=1;
gg_nlp(77,90)=-Fc(2)/Fc(4);
gg_nlp(77,92)=Fc(2)/Fc(4);
gg_nlp(77,93)=Fc(3)/Fc(4);
gg_nlp(78,4)=1;
gg_nlp(78,60)=-1+Fh(1)/Fc(4);
gg_nlp(78,64)=-Fh(1)/Fc(4);
gg_nlp(78,66)=Fh(2)/Fc(4);
gg_nlp(78,69)=-Fh(2)/Fc(4);
gg_nlp(78,72)=Fh(3)/Fc(4);
gg_nlp(78,74)=-Fh(3)/Fc(4);
gg_nlp(78,77)=Fh(4)/Fc(4);
gg_nlp(78,78)=-Fh(4)/Fc(4);
gg_nlp(78,80)=Fh(5)/Fc(4);
gg_nlp(78,81)=-Fh(5)/Fc(4);
gg_nlp(78,88)=-Fc(3)/Fc(4);
gg_nlp(78,89)=1;
gg_nlp(78,90)=-Fc(2)/Fc(4);
gg_nlp(78,92)=Fc(2)/Fc(4);
gg_nlp(78,93)=Fc(3)/Fc(4);
gg_nlp(79,6)=1;
gg_nlp(79,61)=-1+Fh(1)/Fc(4);
gg_nlp(79,64)=-Fh(1)/Fc(4);
gg_nlp(79,67)=Fh(2)/Fc(4);
gg_nlp(79,69)=-Fh(2)/Fc(4);
gg_nlp(79,72)=Fh(3)/Fc(4);
gg_nlp(79,74)=-Fh(3)/Fc(4);
gg_nlp(79,77)=Fh(4)/Fc(4);
gg_nlp(79,78)=-Fh(4)/Fc(4);
gg_nlp(79,80)=Fh(5)/Fc(4);
gg_nlp(79,81)=-Fh(5)/Fc(4);
gg_nlp(79,88)=-Fc(3)/Fc(4);
gg_nlp(79,89)=1;
gg_nlp(79,90)=-Fc(2)/Fc(4);
gg_nlp(79,92)=Fc(2)/Fc(4);
gg_nlp(79,93)=Fc(3)/Fc(4);
gg_nlp(80,5)=1;
gg_nlp(80,60)=Fh(1)/Fc(4);
gg_nlp(80,64)=-Fh(1)/Fc(4);
gg_nlp(80,66)=-1+Fh(2)/Fc(4);
gg_nlp(80,69)=-Fh(2)/Fc(4);
gg_nlp(80,72)=Fh(3)/Fc(4);
gg_nlp(80,74)=-Fh(3)/Fc(4);
gg_nlp(80,77)=Fh(4)/Fc(4);
gg_nlp(80,78)=-Fh(4)/Fc(4);
gg_nlp(80,80)=Fh(5)/Fc(4);
gg_nlp(80,81)=-Fh(5)/Fc(4);
gg_nlp(80,88)=-Fc(3)/Fc(4);
gg_nlp(80,89)=1;
gg_nlp(80,90)=-Fc(2)/Fc(4);
gg_nlp(80,92)=Fc(2)/Fc(4);
gg_nlp(80,93)=Fc(3)/Fc(4);
gg_nlp(81,7)=1;
gg_nlp(81,61)=Fh(1)/Fc(4);
gg_nlp(81,64)=-Fh(1)/Fc(4);
gg_nlp(81,67)=-1+Fh(2)/Fc(4);
gg_nlp(81,69)=-Fh(2)/Fc(4);
gg_nlp(81,72)=Fh(3)/Fc(4);
gg_nlp(81,74)=-Fh(3)/Fc(4);
gg_nlp(81,77)=Fh(4)/Fc(4);
gg_nlp(81,78)=-Fh(4)/Fc(4);
gg_nlp(81,80)=Fh(5)/Fc(4);
gg_nlp(81,81)=-Fh(5)/Fc(4);
gg_nlp(81,88)=-Fc(3)/Fc(4);
gg_nlp(81,89)=1;
gg_nlp(81,90)=-Fc(2)/Fc(4);
gg_nlp(81,92)=Fc(2)/Fc(4);
gg_nlp(81,93)=Fc(3)/Fc(4);
gg_nlp(82,8)=1;
gg_nlp(82,61)=-1;
gg_nlp(82,88)=1;
gg_nlp(83,12)=1;
gg_nlp(83,62)=-1+Fh(1)/Fc(3);
gg_nlp(83,64)=-Fh(1)/Fc(3);
gg_nlp(83,68)=Fh(2)/Fc(3);
gg_nlp(83,69)=-Fh(2)/Fc(3);
gg_nlp(83,73)=Fh(3)/Fc(3);
gg_nlp(83,74)=-Fh(3)/Fc(3);
gg_nlp(83,77)=Fh(4)/Fc(3);
gg_nlp(83,78)=-Fh(4)/Fc(3);
gg_nlp(83,80)=Fh(5)/Fc(3);
gg_nlp(83,81)=-Fh(5)/Fc(3);
gg_nlp(83,90)=-Fc(2)/Fc(3);
gg_nlp(83,92)=Fc(2)/Fc(3);
gg_nlp(83,93)=1;
gg_nlp(84,6)=1;
gg_nlp(84,61)=-1+Fh(1)/Fc(4);
gg_nlp(84,64)=-Fh(1)/Fc(4);
gg_nlp(84,67)=Fh(2)/Fc(4);
gg_nlp(84,69)=-Fh(2)/Fc(4);
gg_nlp(84,72)=Fh(3)/Fc(4);
gg_nlp(84,74)=-Fh(3)/Fc(4);
gg_nlp(84,77)=Fh(4)/Fc(4);
gg_nlp(84,78)=-Fh(4)/Fc(4);
gg_nlp(84,80)=Fh(5)/Fc(4);
gg_nlp(84,81)=-Fh(5)/Fc(4);
gg_nlp(84,88)=-Fc(3)/Fc(4);
gg_nlp(84,89)=1;
gg_nlp(84,90)=-Fc(2)/Fc(4);
gg_nlp(84,92)=Fc(2)/Fc(4);
gg_nlp(84,93)=Fc(3)/Fc(4);
gg_nlp(85,13)=1;
gg_nlp(85,62)=-1;
gg_nlp(85,89)=1;
gg_nlp(86,9)=1;
gg_nlp(86,67)=-1;
gg_nlp(86,88)=1;
gg_nlp(87,14)=1;
gg_nlp(87,62)=Fh(1)/Fc(3);
gg_nlp(87,64)=-Fh(1)/Fc(3);
gg_nlp(87,68)=-1+Fh(2)/Fc(3);
gg_nlp(87,69)=-Fh(2)/Fc(3);
gg_nlp(87,73)=Fh(3)/Fc(3);
gg_nlp(87,74)=-Fh(3)/Fc(3);
gg_nlp(87,77)=Fh(4)/Fc(3);
gg_nlp(87,78)=-Fh(4)/Fc(3);
gg_nlp(87,80)=Fh(5)/Fc(3);
gg_nlp(87,81)=-Fh(5)/Fc(3);
gg_nlp(87,90)=-Fc(2)/Fc(3);
gg_nlp(87,92)=Fc(2)/Fc(3);
gg_nlp(87,93)=1;
gg_nlp(88,7)=1;
gg_nlp(88,61)=Fh(1)/Fc(4);
gg_nlp(88,64)=-Fh(1)/Fc(4);
gg_nlp(88,67)=-1+Fh(2)/Fc(4);
gg_nlp(88,69)=-Fh(2)/Fc(4);
gg_nlp(88,72)=Fh(3)/Fc(4);
gg_nlp(88,74)=-Fh(3)/Fc(4);
gg_nlp(88,77)=Fh(4)/Fc(4);
gg_nlp(88,78)=-Fh(4)/Fc(4);
gg_nlp(88,80)=Fh(5)/Fc(4);
gg_nlp(88,81)=-Fh(5)/Fc(4);
gg_nlp(88,88)=-Fc(3)/Fc(4);
gg_nlp(88,89)=1;
gg_nlp(88,90)=-Fc(2)/Fc(4);
gg_nlp(88,92)=Fc(2)/Fc(4);
gg_nlp(88,93)=Fc(3)/Fc(4);
gg_nlp(89,15)=1;
gg_nlp(89,68)=-1;
gg_nlp(89,89)=1;
gg_nlp(90,10)=1;
gg_nlp(90,72)=-1;
gg_nlp(90,88)=1;
gg_nlp(91,16)=1;
gg_nlp(91,62)=Fh(1)/Fc(3);
gg_nlp(91,64)=-Fh(1)/Fc(3);
gg_nlp(91,68)=Fh(2)/Fc(3);
gg_nlp(91,69)=-Fh(2)/Fc(3);
gg_nlp(91,73)=-1+Fh(3)/Fc(3);
gg_nlp(91,74)=-Fh(3)/Fc(3);
gg_nlp(91,77)=Fh(4)/Fc(3);
gg_nlp(91,78)=-Fh(4)/Fc(3);
gg_nlp(91,80)=Fh(5)/Fc(3);
gg_nlp(91,81)=-Fh(5)/Fc(3);
gg_nlp(91,90)=-Fc(2)/Fc(3);
gg_nlp(91,92)=Fc(2)/Fc(3);
gg_nlp(91,93)=1;
gg_nlp(92,11)=1;
gg_nlp(92,61)=Fh(1)/Fc(4);
gg_nlp(92,64)=-Fh(1)/Fc(4);
gg_nlp(92,67)=Fh(2)/Fc(4);
gg_nlp(92,69)=-Fh(2)/Fc(4);
gg_nlp(92,72)=-1+Fh(3)/Fc(4);
gg_nlp(92,74)=-Fh(3)/Fc(4);
gg_nlp(92,77)=Fh(4)/Fc(4);
gg_nlp(92,78)=-Fh(4)/Fc(4);
gg_nlp(92,80)=Fh(5)/Fc(4);
gg_nlp(92,81)=-Fh(5)/Fc(4);
gg_nlp(92,88)=-Fc(3)/Fc(4);
gg_nlp(92,89)=1;
gg_nlp(92,90)=-Fc(2)/Fc(4);
gg_nlp(92,92)=Fc(2)/Fc(4);
gg_nlp(92,93)=Fc(3)/Fc(4);
gg_nlp(93,17)=1;
gg_nlp(93,73)=-1;
gg_nlp(93,89)=1;
gg_nlp(94,18)=1;
gg_nlp(94,62)=-1;
gg_nlp(94,90)=1;
gg_nlp(95,21)=1;
gg_nlp(95,63)=-1;
gg_nlp(95,91)=1;
gg_nlp(96,12)=1;
gg_nlp(96,62)=-1+Fh(1)/Fc(3);
gg_nlp(96,64)=-Fh(1)/Fc(3);
gg_nlp(96,68)=Fh(2)/Fc(3);
gg_nlp(96,69)=-Fh(2)/Fc(3);
gg_nlp(96,73)=Fh(3)/Fc(3);
gg_nlp(96,74)=-Fh(3)/Fc(3);
gg_nlp(96,77)=Fh(4)/Fc(3);
gg_nlp(96,78)=-Fh(4)/Fc(3);
gg_nlp(96,80)=Fh(5)/Fc(3);
gg_nlp(96,81)=-Fh(5)/Fc(3);
gg_nlp(96,90)=-Fc(2)/Fc(3);
gg_nlp(96,92)=Fc(2)/Fc(3);
gg_nlp(96,93)=1;
gg_nlp(97,22)=1;
gg_nlp(97,63)=-1+Fh(1)/Fc(3);
gg_nlp(97,64)=-Fh(1)/Fc(3);
gg_nlp(97,77)=Fh(4)/Fc(3);
gg_nlp(97,78)=-Fh(4)/Fc(3);
gg_nlp(97,80)=Fh(5)/Fc(3);
gg_nlp(97,81)=-Fh(5)/Fc(3);
gg_nlp(97,91)=-Fc(2)/Fc(3);
gg_nlp(97,92)=Fc(2)/Fc(3);
gg_nlp(97,93)=1;
gg_nlp(98,19)=1;
gg_nlp(98,68)=-1;
gg_nlp(98,90)=1;
gg_nlp(99,23)=1;
gg_nlp(99,69)=-1;
gg_nlp(99,91)=1;
gg_nlp(100,20)=1;
gg_nlp(100,73)=-1;
gg_nlp(100,90)=1;
gg_nlp(101,24)=1;
gg_nlp(101,74)=-1;
gg_nlp(101,91)=1;
gg_nlp(102,21)=1;
gg_nlp(102,63)=-1;
gg_nlp(102,91)=1;
gg_nlp(103,29)=1;
gg_nlp(103,64)=-1;
gg_nlp(103,92)=1;
gg_nlp(104,22)=1;
gg_nlp(104,63)=-1+Fh(1)/Fc(3);
gg_nlp(104,64)=-Fh(1)/Fc(3);
gg_nlp(104,77)=Fh(4)/Fc(3);
gg_nlp(104,78)=-Fh(4)/Fc(3);
gg_nlp(104,80)=Fh(5)/Fc(3);
gg_nlp(104,81)=-Fh(5)/Fc(3);
gg_nlp(104,91)=-Fc(2)/Fc(3);
gg_nlp(104,92)=Fc(2)/Fc(3);
gg_nlp(104,93)=1;
gg_nlp(105,30)=1;
gg_nlp(105,64)=-1;
gg_nlp(105,93)=1;
gg_nlp(106,25)=1;
gg_nlp(106,77)=-1;
gg_nlp(106,91)=1;
gg_nlp(107,31)=1;
gg_nlp(107,78)=-1;
gg_nlp(107,92)=1;
gg_nlp(108,26)=1;
gg_nlp(108,63)=Fh(1)/Fc(3);
gg_nlp(108,64)=-Fh(1)/Fc(3);
gg_nlp(108,77)=-1+Fh(4)/Fc(3);
gg_nlp(108,78)=-Fh(4)/Fc(3);
gg_nlp(108,80)=Fh(5)/Fc(3);
gg_nlp(108,81)=-Fh(5)/Fc(3);
gg_nlp(108,91)=-Fc(2)/Fc(3);
gg_nlp(108,92)=Fc(2)/Fc(3);
gg_nlp(108,93)=1;
gg_nlp(109,32)=1;
gg_nlp(109,78)=-1;
gg_nlp(109,93)=1;
gg_nlp(110,27)=1;
gg_nlp(110,80)=-1;
gg_nlp(110,91)=1;
gg_nlp(111,33)=1;
gg_nlp(111,81)=-1;
gg_nlp(111,92)=1;
gg_nlp(112,28)=1;
gg_nlp(112,63)=Fh(1)/Fc(3);
gg_nlp(112,64)=-Fh(1)/Fc(3);
gg_nlp(112,77)=Fh(4)/Fc(3);
gg_nlp(112,78)=-Fh(4)/Fc(3);
gg_nlp(112,80)=-1+Fh(5)/Fc(3);
gg_nlp(112,81)=-Fh(5)/Fc(3);
gg_nlp(112,91)=-Fc(2)/Fc(3);
gg_nlp(112,92)=Fc(2)/Fc(3);
gg_nlp(112,93)=1;
gg_nlp(113,34)=1;
gg_nlp(113,81)=-1;
gg_nlp(113,93)=1;
gg_nlp(114,35)=1;
gg_nlp(114,64)=-1+Fh(1)/Fc(1);
gg_nlp(114,65)=-Fh(1)/Fc(1);
gg_nlp(114,70)=Fh(2)/Fc(1);
gg_nlp(114,71)=-Fh(2)/Fc(1);
gg_nlp(114,75)=Fh(3)/Fc(1);
gg_nlp(114,76)=-Fh(3)/Fc(1);
gg_nlp(114,78)=Fh(4)/Fc(1);
gg_nlp(114,79)=-Fh(4)/Fc(1);
gg_nlp(114,81)=Fh(5)/Fc(1);
gg_nlp(114,82)=-Fh(5)/Fc(1);
gg_nlp(115,40)=1;
gg_nlp(115,65)=-1;
gg_nlp(116,36)=1;
gg_nlp(116,64)=Fh(1)/Fc(1);
gg_nlp(116,65)=-Fh(1)/Fc(1);
gg_nlp(116,70)=-1+Fh(2)/Fc(1);
gg_nlp(116,71)=-Fh(2)/Fc(1);
gg_nlp(116,75)=Fh(3)/Fc(1);
gg_nlp(116,76)=-Fh(3)/Fc(1);
gg_nlp(116,78)=Fh(4)/Fc(1);
gg_nlp(116,79)=-Fh(4)/Fc(1);
gg_nlp(116,81)=Fh(5)/Fc(1);
gg_nlp(116,82)=-Fh(5)/Fc(1);
gg_nlp(117,41)=1;
gg_nlp(117,71)=-1;
gg_nlp(118,37)=1;
gg_nlp(118,64)=Fh(1)/Fc(1);
gg_nlp(118,65)=-Fh(1)/Fc(1);
gg_nlp(118,70)=Fh(2)/Fc(1);
gg_nlp(118,71)=-Fh(2)/Fc(1);
gg_nlp(118,75)=-1+Fh(3)/Fc(1);
gg_nlp(118,76)=-Fh(3)/Fc(1);
gg_nlp(118,78)=Fh(4)/Fc(1);
gg_nlp(118,79)=-Fh(4)/Fc(1);
gg_nlp(118,81)=Fh(5)/Fc(1);
gg_nlp(118,82)=-Fh(5)/Fc(1);
gg_nlp(119,42)=1;
gg_nlp(119,76)=-1;
gg_nlp(120,38)=1;
gg_nlp(120,64)=Fh(1)/Fc(1);
gg_nlp(120,65)=-Fh(1)/Fc(1);
gg_nlp(120,70)=Fh(2)/Fc(1);
gg_nlp(120,71)=-Fh(2)/Fc(1);
gg_nlp(120,75)=Fh(3)/Fc(1);
gg_nlp(120,76)=-Fh(3)/Fc(1);
gg_nlp(120,78)=-1+Fh(4)/Fc(1);
gg_nlp(120,79)=-Fh(4)/Fc(1);
gg_nlp(120,81)=Fh(5)/Fc(1);
gg_nlp(120,82)=-Fh(5)/Fc(1);
gg_nlp(121,43)=1;
gg_nlp(121,79)=-1;
gg_nlp(122,39)=1;
gg_nlp(122,64)=Fh(1)/Fc(1);
gg_nlp(122,65)=-Fh(1)/Fc(1);
gg_nlp(122,70)=Fh(2)/Fc(1);
gg_nlp(122,71)=-Fh(2)/Fc(1);
gg_nlp(122,75)=Fh(3)/Fc(1);
gg_nlp(122,76)=-Fh(3)/Fc(1);
gg_nlp(122,78)=Fh(4)/Fc(1);
gg_nlp(122,79)=-Fh(4)/Fc(1);
gg_nlp(122,81)=-1+Fh(5)/Fc(1);
gg_nlp(122,82)=-Fh(5)/Fc(1);
gg_nlp(123,44)=1;
gg_nlp(123,82)=-1;
gg_nlp(124,50)=1;
gg_nlp(124,65)=-1;
gg_nlp(125,51)=1;
gg_nlp(125,71)=-1;
gg_nlp(126,52)=1;
gg_nlp(126,76)=-1;
gg_nlp(127,53)=1;
gg_nlp(127,79)=-1;
gg_nlp(128,54)=1;
gg_nlp(128,82)=-1;
gg_nlp(129,45)=1;
gg_nlp(129,83)=1;
gg_nlp(130,46)=1;
gg_nlp(130,84)=1;
gg_nlp(131,47)=1;
gg_nlp(131,85)=1;
gg_nlp(132,48)=1;
gg_nlp(132,86)=1;
gg_nlp(133,49)=1;
gg_nlp(133,59)=-Fh(1)/Fc(5);
gg_nlp(133,87)=1;
gg_nlp(134,1)=-1;
gg_nlp(135,3)=-1;
gg_nlp(136,4)=-1;
gg_nlp(137,5)=-1;
gg_nlp(138,8)=-1;
gg_nlp(139,6)=-1;
gg_nlp(140,9)=-1;
gg_nlp(141,7)=-1;
gg_nlp(142,10)=-1;
gg_nlp(143,11)=-1;
gg_nlp(144,18)=-1;
gg_nlp(145,12)=-1;
gg_nlp(146,19)=-1;
gg_nlp(147,20)=-1;
gg_nlp(148,21)=-1;
gg_nlp(149,22)=-1;
gg_nlp(150,25)=-1;
gg_nlp(151,26)=-1;
gg_nlp(152,27)=-1;
gg_nlp(153,28)=-1;
gg_nlp(154,35)=-1;
gg_nlp(155,36)=-1;
gg_nlp(156,37)=-1;
gg_nlp(157,38)=-1;
gg_nlp(158,39)=-1;

if nargout > 2 || size(x,2) > 93

C = CF*ones(1,35);
    
B = zeros(158,35);
B(39,1)=-omega(1);
B(40,2)=-omega(1);
B(41,3)=-omega(1);
B(42,4)=-omega(2);
B(43,5)=-omega(1);
B(44,6)=-omega(1);
B(45,7)=-omega(2);
B(46,8)=-omega(2);
B(47,9)=-omega(3);
B(48,10)=-omega(3);
B(49,11)=-omega(1);
B(50,12)=-omega(1);
B(51,13)=-omega(2);
B(52,14)=-omega(3);
B(53,15)=-omega(1);
B(54,16)=-omega(1);
B(55,17)=-omega(4);
B(56,18)=-omega(4);
B(57,19)=-omega(5);
B(58,20)=-omega(5);
B(59,21)=-omega(1);
B(60,22)=-omega(2);
B(61,23)=-omega(3);
B(62,24)=-omega(4);
B(63,25)=-omega(5);
B(64,31)=-omega(1);
B(65,32)=-omega(2);
B(66,33)=-omega(3);
B(67,34)=-omega(4);
B(68,35)=-omega(5);
B(69,26)=-omega(1);
B(70,27)=-omega(2);
B(71,28)=-omega(3);
B(72,29)=-omega(4);
B(73,30)=-omega(5);
B(74,1)=alfa;
B(75,1)=alfa;
B(76,2)=alfa;
B(77,2)=alfa;
B(78,3)=alfa;
B(79,3)=alfa;
B(80,4)=alfa;
B(81,4)=alfa;
B(82,5)=alfa;
B(83,5)=alfa;
B(84,6)=alfa;
B(85,6)=alfa;
B(86,7)=alfa;
B(87,7)=alfa;
B(88,8)=alfa;
B(89,8)=alfa;
B(90,9)=alfa;
B(91,9)=alfa;
B(92,10)=alfa;
B(93,10)=alfa;
B(94,11)=alfa;
B(95,11)=alfa;
B(96,12)=alfa;
B(97,12)=alfa;
B(98,13)=alfa;
B(99,13)=alfa;
B(100,14)=alfa;
B(101,14)=alfa;
B(102,15)=alfa;
B(103,15)=alfa;
B(104,16)=alfa;
B(105,16)=alfa;
B(106,17)=alfa;
B(107,17)=alfa;
B(108,18)=alfa;
B(109,18)=alfa;
B(110,19)=alfa;
B(111,19)=alfa;
B(112,20)=alfa;
B(113,20)=alfa;
B(114,21)=alfa;
B(115,21)=alfa;
B(116,22)=alfa;
B(117,22)=alfa;
B(118,23)=alfa;
B(119,23)=alfa;
B(120,24)=alfa;
B(121,24)=alfa;
B(122,25)=alfa;
B(123,25)=alfa; 
B(124,31)=alfa;
B(125,32)=alfa;
B(126,33)=alfa;
B(127,34)=alfa;
B(128,35)=alfa;
B(129,26)=alfa;
B(130,27)=alfa;
B(131,28)=alfa;
B(132,29)=alfa;
B(133,30)=alfa;
end

if size(x,2) > 93         % NLP com integralidade relaxada
  gf_nlp = [gf_nlp C];
  gg_nlp = [gg_nlp B];
end
gg_nlp=gg_nlp';
