function [ys,check]=model_steadystate(junk,ys)
check=0;

%X_X_X=ys(5:81,1);
%X_X_X=zeros(81,1);
X_X_X=[0;0;0;0;  1.123387605065928e+000;  5.179507036413639e-001;  -2.091143519676516e+000;  -2.091084680291438e-001;  9.677460374771604e-001;  -3.088418593154907e-001;  2.782058229172606e+000;  1.665316422901766e+000;  -4.844888835320180e+001;  -6.131125740934371e+002;  -2.197185350660502e+001;  1.672907280742844e-001;  1.210987514646240e-001;  -1.626858994754924e+000;  -1.690148684296642e+000;  -1.690643976475363e+000;  -2.431232565387712e+003;  -9.591511060966237e+000;  -5.231510834957766e+001;  5.385062329029318e+000;  -3.744171599743758e+001;  -2.478202334226309e+000;  -1.352446500025474e+001;  -9.195697761996188e-002;  1.115472680678842e-002;  9.999826043494361e-001;  3.120941294610447e+000;  9.879541705828602e-001;  7.142586966255084e-004;  7.600931510321339e-003;  9.998812553819759e-001;  8.366902164083436e-004;  8.322484416469577e-003;  9.981415303129025e-001;  6.935155508087580e-003;  6.211311642759715e-002;  8.881874238471024e+000;  9.169345900818962e-001;  -7.169972136887629e+000;  -6.398110056927791e+000;  -3.980820414258209e+000;  4.809127251183493e-003;  -9.377800190031154e-001;  1.105602725121468e-001;  1.138780245221088e+000;  -4.929560420175079e+000;  -4.783352147569248e+000;  -2.666706119444438e+000;  9.923030234072799e-003;  1.450048601531899e-002;  6.260435988025287e+000;  -1.097171882853586e+001;  -8.566816034897027e+000;  -8.315121923669379e+000;  4.762406702957432e+000;  -1.661115398052833e+000;  5.008163477379949e-001;  8.503981005525445e-001;  -6.645533937660539e-002;  -5.378159309961514e+000;  -5.481819323520875e+000;  -3.374684649387517e+000;  4.941943129354257e-002;  4.999539609798843e-002;  5.000089961812475e-002;  5.000174665402529e-002;  -7.081361982070511e+000;  9.511673657695712e-001;  9.500068331153461e-001;  9.499986646450127e-001;  9.499974075593776e-001;  8.081362004176363e+000;  -4.166861500711246e-003;  -3.595318443931717e+001;  -2.364098094490368e+003;  -1.343869753255159e+001;  5.776950702345137e-003];
X_X_X=X_X_X(5:end,1);

[X_X_X,fval]=fminunc(@sub_ono_0,X_X_X,optimset('MaxFunEvals',40000,'MaxIter',10^9,'GradObj','on','TolFun',1.0e-7,'Display','iter'));
%[X_X_X,fval]=fminsearch(@sub_ono,X_X_X,optimset('MaxFunEvals',4000,'MaxIter',10^9,'Display','off'));
%[X_X_X,fval]=fminsearch(@sub_ono,X_X_X,optimset('MaxFunEvals',40000,'MaxIter',10^9,'Display','off'));
if fval>10^-6
  %[X_X_X,fval]=fminsearch(@sub_ono,X_X_X,optimset('MaxFunEvals',40000,'MaxIter',10^9,'Display','off'));
  if fval>10^-6
    disp(fval);
    %[X_X_X,fval]=fminsearch(@sub_ono,X_X_X,optimset('MaxFunEvals',40000,'MaxIter',10^9,'Display','off'));
    %[X_X_X,fval]=fminsearch(@sub_ono,X_X_X,optimset('MaxFunEvals',40000,'MaxIter',10^9,'Display','off'));
    %disp(fval);
  end
end

global alfa_L1 alfa_L2 alfa_L3 alfa_L4 alfa_O1 alfa_O2 alfa_O3 alfa_O4 bbb fi_B fi_D fi_HB fi_P fi_S fi_Y gama_B gama_G gama_P gama_Y nu_1Y_0 nu_2Y_0 nu_3Y_0 nu_4Y_0 nu_B_0 nu_F1_0 nu_F2_0 nu_F3_0 nu_F4_0 nu_FB_0 nu_FD_0 nu_FO_0 nu_FP_0 nu_FR_0 nu_FS_0 nu_FY_0 nu_G_0 nu_H_0 nu_L_0 nu_M_0 nu_O_0 nu_R_0 nu_T_0 nu_TO_0 omega omega_2 omega_H psi sigma teta; 
z_1Y=nu_1Y_0;
z_2Y=nu_2Y_0;
z_3Y=nu_3Y_0;
z_4Y=nu_4Y_0;
z_B=nu_B_0;
z_F1=nu_F1_0;
z_F2=nu_F2_0;
z_F3=nu_F3_0;
z_F4=nu_F4_0;
z_FB=nu_FB_0;
z_FD=nu_FD_0;
z_FO=nu_FO_0;
z_FP=nu_FP_0;
z_FR=nu_FR_0;
z_FS=nu_FS_0;
z_FY=nu_FY_0;
z_G=nu_G_0;
z_H=nu_H_0;
z_L=nu_L_0;
z_M=nu_M_0;
z_O=nu_O_0;
z_R=nu_R_0;
z_T=nu_T_0;
z_TO=nu_TO_0;


X_X_X=[zeros(4,1);X_X_X];
b_1=X_X_X(1,1);
b_2=X_X_X(2,1);
b_3=X_X_X(3,1);
b_4=X_X_X(4,1);
b_FD=X_X_X(5,1);
b_FF=X_X_X(6,1);
b_GD=X_X_X(7,1);
b_GF=X_X_X(8,1);
b_HD=X_X_X(9,1);
b_HF=X_X_X(10,1);
c=X_X_X(11,1);
d_1=X_X_X(12,1);
d_2=X_X_X(13,1);
d_3=X_X_X(14,1);
d_4=X_X_X(15,1);
d_F=X_X_X(16,1);
f=X_X_X(17,1);
g=X_X_X(18,1);
i=X_X_X(19,1);
i_1=X_X_X(20,1);
i_2=X_X_X(21,1);
i_3=X_X_X(22,1);
i_4=X_X_X(23,1);
k_1=X_X_X(24,1);
k_2=X_X_X(25,1);
k_3=X_X_X(26,1);
k_4=X_X_X(27,1);
l=X_X_X(28,1);
limda=X_X_X(29,1);
limda_1b=X_X_X(30,1);
limda_1p=X_X_X(31,1);
limda_2b=X_X_X(32,1);
limda_2d=X_X_X(33,1);
limda_2p=X_X_X(34,1);
limda_3b=X_X_X(35,1);
limda_3d=X_X_X(36,1);
limda_3p=X_X_X(37,1);
limda_4b=X_X_X(38,1);
limda_4d=X_X_X(39,1);
limda_4p=X_X_X(40,1);
m=X_X_X(41,1);
o_1=X_X_X(42,1);
o_2=X_X_X(43,1);
o_3=X_X_X(44,1);
o_4=X_X_X(45,1);
p=X_X_X(46,1);
p_D=X_X_X(47,1);
p_F=X_X_X(48,1);
p_O=X_X_X(49,1);
p_2=X_X_X(50,1);
p_3=X_X_X(51,1);
p_4=X_X_X(52,1);
r_D=X_X_X(53,1);
r_F=X_X_X(54,1);
s_1=X_X_X(55,1);
s_2=X_X_X(56,1);
s_3=X_X_X(57,1);
s_4=X_X_X(58,1);
s_F=X_X_X(59,1);
t_L=X_X_X(60,1);
t_O=X_X_X(61,1);
w=X_X_X(62,1);
w_1=X_X_X(63,1);
w_2=X_X_X(64,1);
w_3=X_X_X(65,1);
w_4=X_X_X(66,1);
x_F1=X_X_X(67,1);
x_F2=X_X_X(68,1);
x_F3=X_X_X(69,1);
x_F4=X_X_X(70,1);
x_FF=X_X_X(71,1);
x_H1=X_X_X(72,1);
x_H2=X_X_X(73,1);
x_H3=X_X_X(74,1);
x_H4=X_X_X(75,1);
x_HF=X_X_X(76,1);
y=X_X_X(77,1);
y_2=X_X_X(78,1);
y_3=X_X_X(79,1);
y_4=X_X_X(80,1);
y_F=X_X_X(81,1);



ys=[b_1;b_2;b_3;b_4;b_FD;b_FF;b_GD;b_GF;b_HD;b_HF;c;d_1;d_2;d_3;d_4;d_F;f;g;i;i_1;i_2;i_3;i_4;k_1;k_2;k_3;k_4;l;limda;limda_1b;limda_1p;limda_2b;limda_2d;limda_2p;limda_3b;limda_3d;limda_3p;limda_4b;limda_4d;limda_4p;m;o_1;o_2;o_3;o_4;p;p_D;p_F;p_O;p_2;p_3;p_4;r_D;r_F;s_1;s_2;s_3;s_4;s_F;t_L;t_O;w;w_1;w_2;w_3;w_4;x_F1;x_F2;x_F3;x_F4;x_FF;x_H1;x_H2;x_H3;x_H4;x_HF;y;y_2;y_3;y_4;y_F;z_1Y;z_2Y;z_3Y;z_4Y;z_B;z_F1;z_F2;z_F3;z_F4;z_FB;z_FD;z_FO;z_FP;z_FR;z_FS;z_FY;z_G;z_H;z_L;z_M;z_O;z_R;z_T;z_TO];
return



function [UU,DDD]=sub_ono_0(X_X_X)

global alfa_L1 alfa_L2 alfa_L3 alfa_L4 alfa_O1 alfa_O2 alfa_O3 alfa_O4 bbb fi_B fi_D fi_HB fi_P fi_S fi_Y gama_B gama_G gama_P gama_Y nu_1Y_0 nu_2Y_0 nu_3Y_0 nu_4Y_0 nu_B_0 nu_F1_0 nu_F2_0 nu_F3_0 nu_F4_0 nu_FB_0 nu_FD_0 nu_FO_0 nu_FP_0 nu_FR_0 nu_FS_0 nu_FY_0 nu_G_0 nu_H_0 nu_L_0 nu_M_0 nu_O_0 nu_R_0 nu_T_0 nu_TO_0 omega omega_2 omega_H psi sigma teta; 
z_1Y=nu_1Y_0;
z_2Y=nu_2Y_0;
z_3Y=nu_3Y_0;
z_4Y=nu_4Y_0;
z_B=nu_B_0;
z_F1=nu_F1_0;
z_F2=nu_F2_0;
z_F3=nu_F3_0;
z_F4=nu_F4_0;
z_FB=nu_FB_0;
z_FD=nu_FD_0;
z_FO=nu_FO_0;
z_FP=nu_FP_0;
z_FR=nu_FR_0;
z_FS=nu_FS_0;
z_FY=nu_FY_0;
z_G=nu_G_0;
z_H=nu_H_0;
z_L=nu_L_0;
z_M=nu_M_0;
z_O=nu_O_0;
z_R=nu_R_0;
z_T=nu_T_0;
z_TO=nu_TO_0;

X_X_X=[zeros(4,1);X_X_X];
b_1=X_X_X(1,1);
b_2=X_X_X(2,1);
b_3=X_X_X(3,1);
b_4=X_X_X(4,1);
b_FD=X_X_X(5,1);
b_FF=X_X_X(6,1);
b_GD=X_X_X(7,1);
b_GF=X_X_X(8,1);
b_HD=X_X_X(9,1);
b_HF=X_X_X(10,1);
c=X_X_X(11,1);
d_1=X_X_X(12,1);
d_2=X_X_X(13,1);
d_3=X_X_X(14,1);
d_4=X_X_X(15,1);
d_F=X_X_X(16,1);
f=X_X_X(17,1);
g=X_X_X(18,1);
i=X_X_X(19,1);
i_1=X_X_X(20,1);
i_2=X_X_X(21,1);
i_3=X_X_X(22,1);
i_4=X_X_X(23,1);
k_1=X_X_X(24,1);
k_2=X_X_X(25,1);
k_3=X_X_X(26,1);
k_4=X_X_X(27,1);
l=X_X_X(28,1);
limda=X_X_X(29,1);
limda_1b=X_X_X(30,1);
limda_1p=X_X_X(31,1);
limda_2b=X_X_X(32,1);
limda_2d=X_X_X(33,1);
limda_2p=X_X_X(34,1);
limda_3b=X_X_X(35,1);
limda_3d=X_X_X(36,1);
limda_3p=X_X_X(37,1);
limda_4b=X_X_X(38,1);
limda_4d=X_X_X(39,1);
limda_4p=X_X_X(40,1);
m=X_X_X(41,1);
o_1=X_X_X(42,1);
o_2=X_X_X(43,1);
o_3=X_X_X(44,1);
o_4=X_X_X(45,1);
p=X_X_X(46,1);
p_D=X_X_X(47,1);
p_F=X_X_X(48,1);
p_O=X_X_X(49,1);
p_2=X_X_X(50,1);
p_3=X_X_X(51,1);
p_4=X_X_X(52,1);
r_D=X_X_X(53,1);
r_F=X_X_X(54,1);
s_1=X_X_X(55,1);
s_2=X_X_X(56,1);
s_3=X_X_X(57,1);
s_4=X_X_X(58,1);
s_F=X_X_X(59,1);
t_L=X_X_X(60,1);
t_O=X_X_X(61,1);
w=X_X_X(62,1);
w_1=X_X_X(63,1);
w_2=X_X_X(64,1);
w_3=X_X_X(65,1);
w_4=X_X_X(66,1);
x_F1=X_X_X(67,1);
x_F2=X_X_X(68,1);
x_F3=X_X_X(69,1);
x_F4=X_X_X(70,1);
x_FF=X_X_X(71,1);
x_H1=X_X_X(72,1);
x_H2=X_X_X(73,1);
x_H3=X_X_X(74,1);
x_H4=X_X_X(75,1);
x_HF=X_X_X(76,1);
y=X_X_X(77,1);
y_2=X_X_X(78,1);
y_3=X_X_X(79,1);
y_4=X_X_X(80,1);
y_F=X_X_X(81,1);


U(1,1)=exp(-c)*(1+exp(z_H)*y^psi-bbb*exp(z_H)*y^psi)-limda;
U(2,1)=exp(z_M-m)+fi_HB*b_HF^2*exp(-3*m)-limda+bbb*limda*exp(p+y);
U(3,1)=(1-bbb*exp(-p-y+r_D))*limda;
U(4,1)=fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f);
U(5,1)=limda*(-exp(s_1)+bbb*exp(s_1)+bbb*exp(d_1));
U(6,1)=limda*(-exp(s_2)+bbb*exp(s_2)+bbb*exp(d_2));
U(7,1)=limda*(-exp(s_3)+bbb*exp(s_3)+bbb*exp(d_3));
U(8,1)=limda*(-exp(s_4)+bbb*exp(s_4)+bbb*exp(d_4));
U(9,1)=-fi_HB*x_HF-limda*exp(s_F)+bbb*limda*(exp(s_F)+exp(d_F));
U(10,1)=-exp(l*(omega+1))+limda*(1-t_L)*exp(w);
U(11,1)=exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p);
U(12,1)=1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D);
U(13,1)=-limda_4b+limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(14,1)=-2*limda_4p*fi_B*b_4/exp(p_4);
U(15,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d;
U(16,1)=-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(17,1)=-limda_4b+limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(18,1)=limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p);
U(19,1)=exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p);
U(20,1)=exp(k_4)-(1-sigma)*exp(k_4-y)-exp(i_4);
U(21,1)=exp(y_4)-exp(-teta*(p_4-p_F)+y_F);
U(22,1)=1-exp(z_4Y-alfa_L4*z_L-alfa_O4*p-y*alfa_L4);
U(23,1)=1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D);
U(24,1)=-limda_3b+limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(25,1)=-2*limda_3p*fi_B*b_3/exp(p_3);
U(26,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d;
U(27,1)=-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(28,1)=-limda_3b+limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(29,1)=limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p);
U(30,1)=exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p);
U(31,1)=exp(k_3)-(1-sigma)*exp(k_3-y)-exp(i_3);
U(32,1)=exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g);
U(33,1)=1-exp(z_3Y-alfa_L3*z_L-alfa_O3*p-y*alfa_L3);
U(34,1)=1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D);
U(35,1)=-limda_2b+limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(36,1)=-2*limda_2p*fi_B*b_2/exp(p_2);
U(37,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d;
U(38,1)=-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(39,1)=-limda_2b+limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(40,1)=limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p);
U(41,1)=exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p);
U(42,1)=exp(k_2)-(1-sigma)*exp(k_2-y)-exp(i_2);
U(43,1)=exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g);
U(44,1)=1-exp(z_2Y-alfa_L2*z_L-alfa_O2*p-y*alfa_L2);
U(45,1)=1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D);
U(46,1)=-limda_1b+limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1);
U(47,1)=-2*limda_1p*fi_B*b_1/exp(p_O);
U(48,1)=-limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+limda_1b;
U(49,1)=-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1);
U(50,1)=exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p);
U(51,1)=exp(k_1)-(1-sigma)*exp(k_1-y)-exp(i_1);
U(52,1)=1-exp(z_1Y-alfa_L1*z_L+alfa_O1*z_O-y*alfa_L1-y*alfa_O1);
U(53,1)=exp(2*p_F+p-f)-exp(z_FP);
U(54,1)=exp(2*y_F+y)-exp(z_FY);
U(55,1)=b_FD*exp(-p_F-y_F)-exp(z_FB);
U(56,1)=x_F1-z_F1;
U(57,1)=x_F2-z_F2;
U(58,1)=x_F3-z_F3;
U(59,1)=x_F4-z_F4;
U(60,1)=exp(d_F-p_F-y_F)-exp(z_FD);
U(61,1)=exp(y+p+2*y_F)-exp(z_FS);
U(62,1)=exp(p-f)-exp(z_FO);
U(63,1)=exp(r_F)-exp(z_FR);
U(64,1)=exp(i)-exp(i_1)-exp(i_2)-exp(i_3)-exp(i_4);
U(65,1)=exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4);
U(66,1)=1-exp(c)-exp(i)-exp(g)-exp(y_F)+exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i));
U(67,1)=1-omega_H*exp((1-teta)*p_D)-(1-omega_H)*exp((1-teta)*p_F);
U(68,1)=exp((1-teta)*p_D)-omega_2*exp((1-teta)*p_2)-(1-omega_2)*exp((1-teta)*p_3);
U(69,1)=x_HF+x_FF-1;
U(70,1)=x_H1+x_F1-1;
U(71,1)=x_H2+x_F2-1;
U(72,1)=x_H3+x_F3-1;
U(73,1)=x_H4+x_F4-1;
U(74,1)=b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4);
U(75,1)=b_GF+b_FF+b_HF;
U(76,1)=exp(p_D*g)+b_GD+b_GF-(t_L*exp(w)+t_O*exp(p_O+o_1)+exp(r_D-p-y)*b_GD+exp(r_F+f-p-y)*b_GF+exp(m)-exp(m-p-y));
U(77,1)=r_D-gama_P*p-gama_Y*y-z_R;
U(78,1)=t_O-z_TO;
U(79,1)=exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T;
U(80,1)=z_B*b_GD-b_GF;
U(81,1)=exp(g)-exp(z_G);


UU=sum(real(U).^2)+10000*sum(imag(U).^2);
if 1==1
  DDD(1,1)=2*(exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))*(exp(p_O+o_1)-exp(r_D+o_1+p_O-y-p))+4*(-limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+limda_1b)*limda_1p*fi_B*b_1*exp(-p_O)+4*(-limda_1b+limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*limda_1p*alfa_L1*exp(-w_1)*fi_B*b_1+4*(-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*fi_B*b_1+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*exp(p_O+o_1)+8*limda_1p^2*fi_B^2*b_1/exp(p_O)^2;
DDD(2,1)=2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*(exp(p_2+y_2)-exp(r_D+y_2+p_2-y-p))+8*(limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p))*limda_2p*exp(-p_2)*fi_B*b_2+4*(-limda_2b+limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*alfa_O2*exp(-o_2)*fi_B*b_2+4*(-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*fi_B*b_2+4*((-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d)*fi_B*b_2*limda_2p+4*(-limda_2b+limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*alfa_L2*exp(-w_2)*fi_B*b_2+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*exp(p_2+y_2)+8*limda_2p^2*fi_B^2*b_2/exp(p_2)^2;
DDD(3,1)=2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*(exp(p_3+y_3)-exp(r_D+y_3+p_3-y-p))+8*(limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p))*limda_3p*exp(-p_3)*fi_B*b_3+4*(-limda_3b+limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*alfa_O3*exp(-o_3)*fi_B*b_3+4*(-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*fi_B*b_3+4*((-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d)*fi_B*b_3*limda_3p+4*(-limda_3b+limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*alfa_L3*exp(-w_3)*fi_B*b_3+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*exp(p_3+y_3)+8*limda_3p^2*fi_B^2*b_3/exp(p_3)^2;
DDD(4,1)=4*((-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d)*fi_B*b_4*limda_4p+4*(-limda_4b+limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*alfa_L4*exp(-w_4)*fi_B*b_4+2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*(exp(p_4+y_4)-exp(r_D+y_4+p_4-y-p))+8*(limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p))*limda_4p*exp(-p_4)*fi_B*b_4+4*(-limda_4b+limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*alfa_O4*exp(-o_4)*fi_B*b_4+4*(-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*fi_B*b_4+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*exp(p_4+y_4)+8*limda_4p^2*fi_B^2*b_4/exp(p_4)^2;
DDD(5,1)=2*b_FD*exp(-2*p_F-2*y_F)-2*exp(-p_F-y_F+z_FB)+2*b_GD+2*b_FD+2*b_HD+2*b_1*exp(p_O+o_1)+2*b_2*exp(p_2+y_2)+2*b_3*exp(p_3+y_3)+2*b_4*exp(p_4+y_4);
DDD(6,1)=2*b_GF+2*b_FF+2*b_HF;
DDD(7,1)=2*(exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T)*gama_B+2*(z_B*b_GD-b_GF)*z_B+2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*(1-exp(-p-y+r_D))+2*b_GD+2*b_FD+2*b_HD+2*b_1*exp(p_O+o_1)+2*b_2*exp(p_2+y_2)+2*b_3*exp(p_3+y_3)+2*b_4*exp(p_4+y_4);
DDD(8,1)=2*(exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T)*gama_B-2*z_B*b_GD+4*b_GF+2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*(1-exp(-p-y+r_F+f))+2*b_FF+2*b_HF;
DDD(9,1)=2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*(1-exp(-p-y+r_D))+2*b_GD+2*b_FD+2*b_HD+2*b_1*exp(p_O+o_1)+2*b_2*exp(p_2+y_2)+2*b_3*exp(p_3+y_3)+2*b_4*exp(p_4+y_4);
DDD(10,1)=4*(exp(z_M-m)+fi_HB*b_HF^2*exp(-3*m)-limda+bbb*limda*exp(p+y))*fi_HB*b_HF*exp(-3*m)+2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*(1-exp(-p-y+r_F+f))+2*(fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f))*fi_HB*exp(-2*m)+2*b_GF+2*b_FF+2*b_HF;
DDD(11,1)=-2*(exp(-c)*(1+exp(z_H)*y^psi-bbb*exp(z_H)*y^psi)-limda)*exp(-c)*(1+exp(z_H)*y^psi-bbb*exp(z_H)*y^psi)+2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*exp(c)-2*(exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g))*exp(-teta*p_2)*omega_H*exp(c)+2*(1-exp(c)-exp(i)-exp(g)-exp(y_F)+exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i)))*(-exp(c)+exp((1-teta)*p_F)*(1-omega_H)*exp(c))-2*(exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g))*exp(-teta*p_3)*omega_H*exp(c);
DDD(12,1)=2*(exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))*exp(d_1)-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*x_H1*exp(d_1)+2*(1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))*(2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p-2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))+2*limda^2*(-exp(s_1)+bbb*exp(s_1)+bbb*exp(d_1))*bbb*exp(d_1);
DDD(13,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*x_H2*exp(d_2)+2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*exp(d_2)+2*(1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))*(2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p-2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))+2*limda^2*(-exp(s_2)+bbb*exp(s_2)+bbb*exp(d_2))*bbb*exp(d_2);
DDD(14,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*x_H3*exp(d_3)+2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*exp(d_3)+2*(1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))*(2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p-2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))+2*limda^2*(-exp(s_3)+bbb*exp(s_3)+bbb*exp(d_3))*bbb*exp(d_3);
DDD(15,1)=2*(1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))*(2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p-2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*x_H4*exp(d_4)+2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*exp(d_4)+2*limda^2*(-exp(s_4)+bbb*exp(s_4)+bbb*exp(d_4))*bbb*exp(d_4);
DDD(16,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*x_HF*exp(d_F)+2*(-fi_HB*x_HF-limda*exp(s_F)+bbb*limda*(exp(s_F)+exp(d_F)))*bbb*limda*exp(d_F)+2*(exp(d_F-p_F-y_F)-exp(z_FD))*exp(d_F-p_F-y_F);
DDD(17,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*b_HF*exp(-p-y+r_F+f)-2*(fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f))*bbb*limda*exp(-p-y+r_F+f)-2*(exp(p-f)-exp(z_FO))*exp(p-f)-2*(exp(2*p_F+p-f)-exp(z_FP))*exp(2*p_F+p-f)-2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*exp(-p-y+r_F+f)*b_GF;
DDD(18,1)=-2*(exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T)*gama_G*exp(p_D+g)-2*(exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g))*exp(-teta*(p_2-p_D)+g)+2*(exp(g)-exp(z_G))*exp(g)+2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*p_D*exp(p_D*g)-2*(1-exp(c)-exp(i)-exp(g)-exp(y_F)+exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i)))*exp(g)-2*(exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g))*exp(-teta*(p_3-p_D)+g);
DDD(19,1)=-2*(exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g))*exp(-teta*p_2)*omega_H*exp(i)+2*(exp(i)-exp(i_1)-exp(i_2)-exp(i_3)-exp(i_4))*exp(i)+2*(1-exp(c)-exp(i)-exp(g)-exp(y_F)+exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i)))*(-exp(i)+exp((1-teta)*p_F)*(1-omega_H)*exp(i))-2*(exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g))*exp(-teta*p_3)*omega_H*exp(i);
DDD(20,1)=2*(exp(d_1)+3*exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p)-exp(k_1)+exp(k_1-y)-exp(k_1-y)*sigma-exp(i)+exp(i_2)+exp(i_3)+exp(i_4))*exp(i_1);
DDD(21,1)=2*exp(i_2)*(-exp(k_2)+exp(k_2-y)-exp(k_2-y)*sigma+3*exp(i_2)+exp(d_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p)-exp(i)+exp(i_1)+exp(i_3)+exp(i_4));
DDD(22,1)=2*exp(i_3)*(-exp(k_3)+exp(k_3-y)-exp(k_3-y)*sigma+3*exp(i_3)+exp(d_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p)-exp(i)+exp(i_1)+exp(i_2)+exp(i_4));
DDD(23,1)=2*exp(i_4)*(-exp(k_4)+exp(k_4-y)-exp(k_4-y)*sigma+3*exp(i_4)+exp(d_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p)-exp(i)+exp(i_1)+exp(i_2)+exp(i_3));
DDD(24,1)=2*(exp(k_1)-(1-sigma)*exp(k_1-y)-exp(i_1))*(exp(k_1)-(1-sigma)*exp(k_1-y))-2*(-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1);
DDD(25,1)=2*(exp(k_2)-(1-sigma)*exp(k_2-y)-exp(i_2))*(exp(k_2)-(1-sigma)*exp(k_2-y))-2*(-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
DDD(26,1)=2*(exp(k_3)-(1-sigma)*exp(k_3-y)-exp(i_3))*(exp(k_3)-(1-sigma)*exp(k_3-y))-2*(-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
DDD(27,1)=2*(exp(k_4)-(1-sigma)*exp(k_4-y)-exp(i_4))*(exp(k_4)-(1-sigma)*exp(k_4-y))-2*(-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
DDD(28,1)=-2*(-exp(l*(omega+1))+limda*(1-t_L)*exp(w))*(omega+1)*exp(l*(omega+1));
DDD(29,1)=2*(exp(z_M-m)+fi_HB*b_HF^2*exp(-3*m)-limda+bbb*limda*exp(p+y))*(-1+bbb*exp(p+y))-2*exp(-c)*(1+exp(z_H)*y^psi-bbb*exp(z_H)*y^psi)+2*limda+2*(-exp(l*(omega+1))+limda*(1-t_L)*exp(w))*(1-t_L)*exp(w)+2*(-fi_HB*x_HF-limda*exp(s_F)+bbb*limda*(exp(s_F)+exp(d_F)))*(-exp(s_F)+bbb*(exp(s_F)+exp(d_F)))+2*(fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f))*(-1-bbb*exp(-p-y+r_F+f))+2*limda*(-exp(s_4)+bbb*exp(s_4)+bbb*exp(d_4))^2+2*limda*(-exp(s_3)+bbb*exp(s_3)+bbb*exp(d_3))^2+2*limda*(-exp(s_2)+bbb*exp(s_2)+bbb*exp(d_2))^2+2*limda*(-exp(s_1)+bbb*exp(s_1)+bbb*exp(d_1))^2+2*(1-bbb*exp(-p-y+r_D))^2*limda;
DDD(30,1)=-2*limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+6*limda_1b-2*limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1)-2+4*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p-4*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D)+2*(-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*(-1+exp(p-r_D)*(1-sigma));
DDD(31,1)=-2*(-limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+limda_1b)*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+2*(-limda_1b+limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1)+2*(1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))*(-2*fi_D*(p+y)*exp(-d_1+o_1)+2*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))+2*(-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1)+8*limda_1p*fi_B^2*b_1^2/exp(p_O)^2;
DDD(32,1)=8*limda_2b+2*limda_2d*teta*exp(-p_2)-2*limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-2*limda_2p*fi_P*p*exp(-r_D+y-p_2+p)-2*limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1)+2*(-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*(-1+exp(p-r_D)*(1-sigma))+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d)*exp(p_2)-2*limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1)-2+4*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p-4*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D);
DDD(33,1)=2*(limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p))*teta*exp(-p_2)-2*(-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p-4*limda_2p*fi_Y*y*exp(y+p-r_D)-2*limda_2b*exp(p_2)+2*limda_2d;
DDD(34,1)=2*(limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p))*(-exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-fi_P*p*exp(-r_D+y-p_2+p))+2*(-limda_2b+limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1)+2*(-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1)+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d)*(-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y+2*fi_Y*y*exp(y+p-r_D))+2*(-limda_2b+limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1)+2*(1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))*(-2*fi_D*(p+y)*exp(-d_2+y_2)+2*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))+8*limda_2p*fi_B^2*b_2^2/exp(p_2)^2;
DDD(35,1)=8*limda_3b+2*limda_3d*teta*exp(-p_3)-2*limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-2*limda_3p*fi_P*p*exp(-r_D+y-p_3+p)-2*limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1)+2*(-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*(-1+exp(p-r_D)*(1-sigma))+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d)*exp(p_3)-2*limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1)-2+4*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p-4*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D);
DDD(36,1)=2*(limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p))*teta*exp(-p_3)-2*(-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p-4*limda_3p*fi_Y*y*exp(y+p-r_D)-2*limda_3b*exp(p_3)+2*limda_3d;
DDD(37,1)=2*(limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p))*(-exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-fi_P*p*exp(-r_D+y-p_3+p))+2*(-limda_3b+limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1)+2*(-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1)+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d)*(-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y+2*fi_Y*y*exp(y+p-r_D))+2*(-limda_3b+limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1)+2*(1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))*(-2*fi_D*(p+y)*exp(-d_3+y_3)+2*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))+8*limda_3p*fi_B^2*b_3^2/exp(p_3)^2;
DDD(38,1)=2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d)*exp(p_4)+8*limda_4b-2*limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)-2+4*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p-4*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D)+2*limda_4d*teta*exp(-p_4)-2*limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-2*limda_4p*fi_P*p*exp(-r_D+y-p_4+p)-2*limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)+2*(-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*(-1+exp(p-r_D)*(1-sigma));
DDD(39,1)=-2*(-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p-4*limda_4p*fi_Y*y*exp(y+p-r_D)-2*limda_4b*exp(p_4)+2*limda_4d+2*(limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p))*teta*exp(-p_4);
DDD(40,1)=2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d)*(-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y+2*fi_Y*y*exp(y+p-r_D))+2*(-limda_4b+limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)+2*(1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))*(-2*fi_D*(p+y)*exp(-d_4+y_4)+2*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))+2*(limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p))*(-exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-fi_P*p*exp(-r_D+y-p_4+p))+2*(-limda_4b+limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)+2*(-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)+8*limda_4p*fi_B^2*b_4^2/exp(p_4)^2;
DDD(41,1)=2*(exp(z_M-m)+fi_HB*b_HF^2*exp(-3*m)-limda+bbb*limda*exp(p+y))*(-exp(z_M-m)-3*fi_HB*b_HF^2*exp(-3*m))+2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*(exp(m)-exp(m-y-p))-4*(fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f))*fi_HB*b_HF*exp(-2*m)+2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*(-exp(m)+exp(m-y-p));
DDD(42,1)=2*(exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))*(b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))+2*(exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T)*t_O*exp(p_O+o_1)-2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w_1+o_1)+2*(1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))*(-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))+2*(-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1)-2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*t_O*exp(p_O+o_1)+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*b_1*exp(p_O+o_1);
DDD(43,1)=2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*exp(o_2+y_2)-2*(-limda_2b+limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
DDD(44,1)=2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*exp(o_3+y_3)-2*(-limda_3b+limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
DDD(45,1)=2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*exp(o_4+y_4)-2*(-limda_4b+limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
DDD(46,1)=2*(limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p))*(-2*limda_3p*exp(-p_3)*fi_P-limda_3p*fi_P*exp(-r_D+y-p_3+p)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p))+2*(-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*(limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(2*fi_P*p+2*fi_D*(p+y)))+2*(-limda_3b+limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*alfa_O3*exp(-o_3)*(2*fi_P*p+2*fi_D*(p+y))+2*(exp(p-f)-exp(z_FO))*exp(p-f)+2*(1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))*(-2*fi_D*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*exp(-d_3+y_3+y+p-r_D)+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))+2*(fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f))*bbb*limda*exp(-p-y+r_F+f)+2*(1-exp(z_4Y-alfa_L4*z_L-alfa_O4*p-y*alfa_L4))*alfa_O4*exp(z_4Y-alfa_L4*z_L-alfa_O4*p-y*alfa_L4)+2*(exp(2*p_F+p-f)-exp(z_FP))*exp(2*p_F+p-f)+2*(-limda_3b+limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*alfa_L3*exp(-w_3)*(2*fi_P*p+2*fi_D*(p+y))+2*(-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*(limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(2*fi_P*p+2*fi_D*(p+y)))+2*(limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p))*(-2*limda_4p*exp(-p_4)*fi_P-limda_4p*fi_P*exp(-r_D+y-p_4+p)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p))+2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*b_4*exp(r_D+y_4+p_4-y-p)+2*(1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))*(-2*fi_D*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*exp(-d_1+o_1+y+p-r_D)+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))+2*(-limda_4b+limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*alfa_O4*exp(-o_4)*(2*fi_P*p+2*fi_D*(p+y))-4*(-limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+limda_1b)*limda_1p*fi_D*(p+y)*exp(-p_O)+4*(-limda_1b+limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*limda_1p*alfa_L1*exp(-w_1)*fi_D*(p+y)+2*(1-exp(z_2Y-alfa_L2*z_L-alfa_O2*p-y*alfa_L2))*alfa_O2*exp(z_2Y-alfa_L2*z_L-alfa_O2*p-y*alfa_L2)+2*(limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p))*(-2*limda_2p*exp(-p_2)*fi_P-limda_2p*fi_P*exp(-r_D+y-p_2+p)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p))+2*(-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*(limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(2*fi_P*p+2*fi_D*(p+y)))+2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*(exp(-p-y+r_D)*b_GD+exp(-p-y+r_F+f)*b_GF-exp(m-y-p))+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d)*((-2*fi_P*p-2*fi_D*(p+y))*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D))+2*(-limda_2b+limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*alfa_O2*exp(-o_2)*(2*fi_P*p+2*fi_D*(p+y))+2*(exp(y+p+2*y_F)-exp(z_FS))*exp(y+p+2*y_F)+2*(-limda_4b+limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*alfa_L4*exp(-w_4)*(2*fi_P*p+2*fi_D*(p+y))+2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*(b_HD*exp(-p-y+r_D)+b_HF*exp(-p-y+r_F+f)+exp(m-y-p))-2*(r_D-gama_P*p-gama_Y*y-z_R)*gama_P+2*(1-bbb*exp(-p-y+r_D))*limda^2*bbb*exp(-p-y+r_D)+2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*b_3*exp(r_D+y_3+p_3-y-p)+2*(-limda_2b+limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*alfa_L2*exp(-w_2)*(2*fi_P*p+2*fi_D*(p+y))+2*(1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))*(-2*fi_D*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*exp(-d_2+y_2+y+p-r_D)+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))+2*(1-exp(z_3Y-alfa_L3*z_L-alfa_O3*p-y*alfa_L3))*alfa_O3*exp(z_3Y-alfa_L3*z_L-alfa_O3*p-y*alfa_L3)+2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*b_2*exp(r_D+y_2+p_2-y-p)+2*(exp(z_M-m)+fi_HB*b_HF^2*exp(-3*m)-limda+bbb*limda*exp(p+y))*bbb*limda*exp(p+y)+2*(exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))*b_1*exp(r_D+o_1+p_O-y-p)+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d)*((-2*fi_P*p-2*fi_D*(p+y))*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D))+2*(-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*(limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1)+2*limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*fi_D*(p+y))+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d)*((-2*fi_P*p-2*fi_D*(p+y))*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D))+2*(1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))*(-2*fi_D*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*exp(-d_4+y_4+y+p-r_D)+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D));
DDD(47,1)=-2*(exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T)*gama_G*exp(p_D+g)-2*(exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g))*teta*exp(-teta*(p_2-p_D)+g)+2*(exp((1-teta)*p_D)-omega_2*exp((1-teta)*p_2)-(1-omega_2)*exp((1-teta)*p_3))*(1-teta)*exp((1-teta)*p_D)-2*(1-omega_H*exp((1-teta)*p_D)-(1-omega_H)*exp((1-teta)*p_F))*omega_H*(1-teta)*exp((1-teta)*p_D)+2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*g*exp(p_D*g)-2*(exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g))*teta*exp(-teta*(p_3-p_D)+g);
DDD(48,1)=-2*(exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g))*teta*exp(-teta*(p_2-p_F)+y_F)-2*(exp(y_4)-exp(-teta*(p_4-p_F)+y_F))*teta*exp(-teta*(p_4-p_F)+y_F)-2*(1-omega_H*exp((1-teta)*p_D)-(1-omega_H)*exp((1-teta)*p_F))*(1-omega_H)*(1-teta)*exp((1-teta)*p_F)-2*(exp(d_F-p_F-y_F)-exp(z_FD))*exp(d_F-p_F-y_F)-2*(b_FD*exp(-p_F-y_F)-exp(z_FB))*b_FD*exp(-p_F-y_F)+4*(exp(2*p_F+p-f)-exp(z_FP))*exp(2*p_F+p-f)+2*(1-exp(c)-exp(i)-exp(g)-exp(y_F)+exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i)))*(1-teta)*exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i))-2*(exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g))*teta*exp(-teta*(p_3-p_F)+y_F);
DDD(49,1)=2*(exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))*(b_1*exp(p_O+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))+2*(exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T)*t_O*exp(p_O+o_1)+2*(-limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+limda_1b)*limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)-2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*t_O*exp(p_O+o_1)+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*b_1*exp(p_O+o_1)-8*limda_1p^2*fi_B^2*b_1^2/exp(p_O)^2;
DDD(50,1)=2*(exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g))*(teta*exp(-teta*p_2)*omega_H*(exp(c)+exp(i))+teta*exp(-teta*(p_2-p_F)+y_F)+teta*exp(-teta*(p_2-p_D)+g))+2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*(b_2*exp(p_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))+2*(limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p))*(-limda_2d*teta*exp(-p_2)+limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)+limda_2p*fi_P*p*exp(-r_D+y-p_2+p))+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d)*limda_2b*exp(p_2)-2*(exp((1-teta)*p_D)-omega_2*exp((1-teta)*p_2)-(1-omega_2)*exp((1-teta)*p_3))*omega_2*(1-teta)*exp((1-teta)*p_2)+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*b_2*exp(p_2+y_2)-8*limda_2p^2*fi_B^2*b_2^2/exp(p_2)^2;
DDD(51,1)=2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*(b_3*exp(p_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))+2*(limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p))*(-limda_3d*teta*exp(-p_3)+limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)+limda_3p*fi_P*p*exp(-r_D+y-p_3+p))+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d)*limda_3b*exp(p_3)-2*(exp((1-teta)*p_D)-omega_2*exp((1-teta)*p_2)-(1-omega_2)*exp((1-teta)*p_3))*(1-omega_2)*(1-teta)*exp((1-teta)*p_3)+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*b_3*exp(p_3+y_3)-8*limda_3p^2*fi_B^2*b_3^2/exp(p_3)^2+2*(exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g))*(teta*exp(-teta*p_3)*omega_H*(exp(c)+exp(i))+teta*exp(-teta*(p_3-p_F)+y_F)+teta*exp(-teta*(p_3-p_D)+g));
DDD(52,1)=2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d)*limda_4b*exp(p_4)+2*(exp(y_4)-exp(-teta*(p_4-p_F)+y_F))*teta*exp(-teta*(p_4-p_F)+y_F)+2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*(b_4*exp(p_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))+2*(limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p))*(-limda_4d*teta*exp(-p_4)+limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)+limda_4p*fi_P*p*exp(-r_D+y-p_4+p))+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*b_4*exp(p_4+y_4)-8*limda_4p^2*fi_B^2*b_4^2/exp(p_4)^2;
DDD(53,1)=2*(limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p))*limda_4p*fi_P*p*exp(-r_D+y-p_4+p)-2*(exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))*b_1*exp(r_D+o_1+p_O-y-p)-2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*exp(-p-y+r_D)*b_GD-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*b_HD*exp(-p-y+r_D)-2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*b_2*exp(r_D+y_2+p_2-y-p)-2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*b_3*exp(r_D+y_3+p_3-y-p)-2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*b_4*exp(r_D+y_4+p_4-y-p)-4*((-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d)*limda_2p*fi_Y*y*exp(y+p-r_D)+2*(limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p))*limda_2p*fi_P*p*exp(-r_D+y-p_2+p)-4*(1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D)-4*((-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d)*limda_4p*fi_Y*y*exp(y+p-r_D)-4*(1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D)-4*(1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D)-4*((-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d)*limda_3p*fi_Y*y*exp(y+p-r_D)+2*(limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p))*limda_3p*fi_P*p*exp(-r_D+y-p_3+p)-4*(1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D)+2*(-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*(-limda_4b*exp(p-r_D)*(1-sigma)-limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))+2*(-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*(-limda_3b*exp(p-r_D)*(1-sigma)-limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))+2*(-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*(-limda_1b*exp(p-r_D)*(1-sigma)-limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))+2*(-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*(-limda_2b*exp(p-r_D)*(1-sigma)-limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))-2*gama_Y*y-2*gama_P*p+2*r_D-2*z_R-2*(1-bbb*exp(-p-y+r_D))*limda^2*bbb*exp(-p-y+r_D);
DDD(54,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*b_HF*exp(-p-y+r_F+f)-2*(fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f))*bbb*limda*exp(-p-y+r_F+f)+2*(exp(r_F)-exp(z_FR))*exp(r_F)-2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*exp(-p-y+r_F+f)*b_GF;
DDD(55,1)=2*limda^2*(-exp(s_1)+bbb*exp(s_1)+bbb*exp(d_1))*exp(s_1)*(-1+bbb);
DDD(56,1)=2*limda^2*(-exp(s_2)+bbb*exp(s_2)+bbb*exp(d_2))*exp(s_2)*(-1+bbb);
DDD(57,1)=2*limda^2*(-exp(s_3)+bbb*exp(s_3)+bbb*exp(d_3))*exp(s_3)*(-1+bbb);
DDD(58,1)=2*limda^2*(-exp(s_4)+bbb*exp(s_4)+bbb*exp(d_4))*exp(s_4)*(-1+bbb);
DDD(59,1)=2*(-fi_HB*x_HF-limda*exp(s_F)+bbb*limda*exp(s_F)+bbb*limda*exp(d_F))*limda*exp(s_F)*(-1+bbb);
DDD(60,1)=2*exp(w)*(limda^2*exp(w)*t_L-x_HF*exp(d_F)-x_H4*exp(d_4)-x_H3*exp(d_3)-x_H2*exp(d_2)-x_H1*exp(d_1)-b_HF*exp(-p-y+r_F+f)-b_HD*exp(-p-y+r_D)+3*exp(w)*t_L-limda^2*exp(w)+limda*exp(l*(omega+1))+exp(-p-y+r_D)*b_GD+exp(-p-y+r_F+f)*b_GF+2*t_O*exp(p_O+o_1)-gama_G*exp(p_D+g)+b_HF+b_HD-z_T-exp(w)+exp(c)+2*exp(m)-2*exp(m-y-p)+gama_B*b_GF+gama_B*b_GD-exp(p_D*g)-b_GD-b_GF);
DDD(61,1)=2*(exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T)*exp(p_O+o_1)+2*t_O-2*z_TO-2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*exp(p_O+o_1);
DDD(62,1)=2*(exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T)*exp(w)*t_L+2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w)+2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*(-exp(w)+exp(w)*t_L)+2*(-exp(l*(omega+1))+limda*(1-t_L)*exp(w))*limda*(1-t_L)*exp(w)-2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*exp(w)*t_L;
DDD(63,1)=2*(exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))*exp(w_1+o_1)-2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w_1+o_1)-2*(-limda_1b+limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1);
DDD(64,1)=-2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w_2+y_2)+2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*exp(w_2+y_2)-2*(-limda_2b+limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
DDD(65,1)=-2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w_3+y_3)+2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*exp(w_3+y_3)-2*(-limda_3b+limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
DDD(66,1)=-2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w_4+y_4)-2*(-limda_4b+limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)+2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*exp(w_4+y_4);
DDD(67,1)=2*x_H1+4*x_F1-2-2*z_F1;
DDD(68,1)=2*x_H2+4*x_F2-2-2*z_F2;
DDD(69,1)=2*x_H3+4*x_F3-2-2*z_F3;
DDD(70,1)=4*x_F4-2*z_F4+2*x_H4-2;
DDD(71,1)=2*x_HF+2*x_FF-2;
DDD(72,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*exp(d_1)+2*x_H1+2*x_F1-2;
DDD(73,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*exp(d_2)+2*x_H2+2*x_F2-2;
DDD(74,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*exp(d_3)+2*x_H3+2*x_F3-2;
DDD(75,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*exp(d_4)+2*x_H4+2*x_F4-2;
DDD(76,1)=-2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*exp(d_F)-2*(-fi_HB*x_HF-limda*exp(s_F)+bbb*limda*(exp(s_F)+exp(d_F)))*fi_HB+2*x_HF+2*x_FF-2;
DDD(77,1)=2*(-limda_3b+limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*alfa_L3*exp(-w_3)*(2*fi_Y*y+2*fi_D*(p+y))+2*(1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))*(-2*fi_D*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*exp(-d_3+y_3+y+p-r_D)+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))+2*(fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f))*bbb*limda*exp(-p-y+r_F+f)+2*(-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*(limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(2*fi_Y*y+2*fi_D*(p+y)))+2*(-limda_4b+limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*alfa_O4*exp(-o_4)*(2*fi_Y*y+2*fi_D*(p+y))+2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*b_4*exp(r_D+y_4+p_4-y-p)+2*(1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))*(-2*fi_D*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*exp(-d_1+o_1+y+p-r_D)+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D))-4*(-limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+limda_1b)*limda_1p*fi_D*(p+y)*exp(-p_O)+4*(-limda_1b+limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*limda_1p*alfa_L1*exp(-w_1)*fi_D*(p+y)-2*(limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p))*limda_4p*fi_P*p*exp(-r_D+y-p_4+p)+2*(exp(p_D*g)+b_GD+b_GF-exp(w)*t_L-t_O*exp(p_O+o_1)-exp(-p-y+r_D)*b_GD-exp(-p-y+r_F+f)*b_GF-exp(m)+exp(m-y-p))*(exp(-p-y+r_D)*b_GD+exp(-p-y+r_F+f)*b_GF-exp(m-y-p))+2*(exp(y+p+2*y_F)-exp(z_FS))*exp(y+p+2*y_F)+2*(-limda_2b+limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*alfa_O2*exp(-o_2)*(2*fi_Y*y+2*fi_D*(p+y))+2*(-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*(limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(2*fi_Y*y+2*fi_D*(p+y)))+2*(exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p))*(b_HD*exp(-p-y+r_D)+b_HF*exp(-p-y+r_F+f)+exp(m-y-p))+2*(exp(2*y_F+y)-exp(z_FY))*exp(2*y_F+y)+2*(-limda_2b+limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*alfa_L2*exp(-w_2)*(2*fi_Y*y+2*fi_D*(p+y))+2*(1-bbb*exp(-p-y+r_D))*limda^2*bbb*exp(-p-y+r_D)-2*(limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p))*limda_2p*fi_P*p*exp(-r_D+y-p_2+p)+2*(-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1))*(limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1)+2*limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*fi_D*(p+y))-2*(r_D-gama_P*p-gama_Y*y-z_R)*gama_Y+2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*b_3*exp(r_D+y_3+p_3-y-p)+2*(-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*(limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(2*fi_Y*y+2*fi_D*(p+y)))+2*(-limda_3b+limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*alfa_O3*exp(-o_3)*(2*fi_Y*y+2*fi_D*(p+y))+2*(1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))*(-2*fi_D*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*exp(-d_2+y_2+y+p-r_D)+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d)*((-2*fi_Y*y-2*fi_D*(p+y)-2*fi_Y)*limda_2p+2*limda_2p*fi_Y*exp(y+p-r_D)+2*limda_2p*fi_Y*y*exp(y+p-r_D))+2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*b_2*exp(r_D+y_2+p_2-y-p)+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d)*((-2*fi_Y*y-2*fi_D*(p+y)-2*fi_Y)*limda_4p+2*limda_4p*fi_Y*exp(y+p-r_D)+2*limda_4p*fi_Y*y*exp(y+p-r_D))+2*(1-exp(z_4Y-alfa_L4*z_L-alfa_O4*p-y*alfa_L4))*alfa_L4*exp(z_4Y-alfa_L4*z_L-alfa_O4*p-y*alfa_L4)+2*(exp(z_M-m)+fi_HB*b_HF^2*exp(-3*m)-limda+bbb*limda*exp(p+y))*bbb*limda*exp(p+y)+2*(exp(-c)*(1+exp(z_H)*y^psi-bbb*exp(z_H)*y^psi)-limda)*exp(-c)*(exp(z_H)*y^psi*psi/y-bbb*exp(z_H)*y^psi*psi/y)+2*(exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p))*b_1*exp(r_D+o_1+p_O-y-p)+2*(-limda_4b+limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*alfa_L4*exp(-w_4)*(2*fi_Y*y+2*fi_D*(p+y))+2*(exp(k_3)-(1-sigma)*exp(k_3-y)-exp(i_3))*(1-sigma)*exp(k_3-y)+2*(exp(k_4)-(1-sigma)*exp(k_4-y)-exp(i_4))*(1-sigma)*exp(k_4-y)+2*(1-exp(z_2Y-alfa_L2*z_L-alfa_O2*p-y*alfa_L2))*alfa_L2*exp(z_2Y-alfa_L2*z_L-alfa_O2*p-y*alfa_L2)+2*(exp(k_2)-(1-sigma)*exp(k_2-y)-exp(i_2))*(1-sigma)*exp(k_2-y)+2*(1-exp(z_3Y-alfa_L3*z_L-alfa_O3*p-y*alfa_L3))*alfa_L3*exp(z_3Y-alfa_L3*z_L-alfa_O3*p-y*alfa_L3)+2*(exp(k_1)-(1-sigma)*exp(k_1-y)-exp(i_1))*(1-sigma)*exp(k_1-y)-2*(1-exp(z_1Y-alfa_L1*z_L+alfa_O1*z_O-y*alfa_L1-y*alfa_O1))*(-alfa_L1-alfa_O1)*exp(z_1Y-alfa_L1*z_L+alfa_O1*z_O-y*alfa_L1-y*alfa_O1)+2*((-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d)*((-2*fi_Y*y-2*fi_D*(p+y)-2*fi_Y)*limda_3p+2*limda_3p*fi_Y*exp(y+p-r_D)+2*limda_3p*fi_Y*y*exp(y+p-r_D))-2*(limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p))*limda_3p*fi_P*p*exp(-r_D+y-p_3+p)+2*(1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))*(-2*fi_D*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*exp(-d_4+y_4+y+p-r_D)+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D));
DDD(78,1)=-2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w_2+y_2)+2*(exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g))*exp(y_2)+2*(exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))*(exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p))+2*(-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1))*limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1)+2*(1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))*(-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D))+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*b_2*exp(p_2+y_2);
DDD(79,1)=-2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w_3+y_3)+2*(exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))*(exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p))+2*(-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1))*limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1)+2*(1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))*(-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D))+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*b_3*exp(p_3+y_3)+2*(exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g))*exp(y_3);
DDD(80,1)=-2*(exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4))*exp(w_4+y_4)+2*(1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))*(-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D))+2*(exp(y_4)-exp(-teta*(p_4-p_F)+y_F))*exp(y_4)+2*(exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))*(exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p))+2*(-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1))*limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1)+2*(b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4))*b_4*exp(p_4+y_4);
DDD(81,1)=-2*(exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g))*exp(-teta*(p_2-p_F)+y_F)-2*(exp(y_4)-exp(-teta*(p_4-p_F)+y_F))*exp(-teta*(p_4-p_F)+y_F)+4*(exp(y+p+2*y_F)-exp(z_FS))*exp(y+p+2*y_F)-2*(exp(d_F-p_F-y_F)-exp(z_FD))*exp(d_F-p_F-y_F)-2*(b_FD*exp(-p_F-y_F)-exp(z_FB))*b_FD*exp(-p_F-y_F)+4*(exp(2*y_F+y)-exp(z_FY))*exp(2*y_F+y)-2*(1-exp(c)-exp(i)-exp(g)-exp(y_F)+exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i)))*exp(y_F)-2*(exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g))*exp(-teta*(p_3-p_F)+y_F);
end
DDD=DDD(5:end,:);
return








function UU=sub_ono(X_X_X)

global alfa_L1 alfa_L2 alfa_L3 alfa_L4 alfa_O1 alfa_O2 alfa_O3 alfa_O4 bbb fi_B fi_D fi_HB fi_P fi_S fi_Y gama_B gama_G gama_P gama_Y nu_1Y_0 nu_2Y_0 nu_3Y_0 nu_4Y_0 nu_B_0 nu_F1_0 nu_F2_0 nu_F3_0 nu_F4_0 nu_FB_0 nu_FD_0 nu_FO_0 nu_FP_0 nu_FR_0 nu_FS_0 nu_FY_0 nu_G_0 nu_H_0 nu_L_0 nu_M_0 nu_O_0 nu_R_0 nu_T_0 nu_TO_0 omega omega_2 omega_H psi sigma teta; 
z_1Y=nu_1Y_0;
z_2Y=nu_2Y_0;
z_3Y=nu_3Y_0;
z_4Y=nu_4Y_0;
z_B=nu_B_0;
z_F1=nu_F1_0;
z_F2=nu_F2_0;
z_F3=nu_F3_0;
z_F4=nu_F4_0;
z_FB=nu_FB_0;
z_FD=nu_FD_0;
z_FO=nu_FO_0;
z_FP=nu_FP_0;
z_FR=nu_FR_0;
z_FS=nu_FS_0;
z_FY=nu_FY_0;
z_G=nu_G_0;
z_H=nu_H_0;
z_L=nu_L_0;
z_M=nu_M_0;
z_O=nu_O_0;
z_R=nu_R_0;
z_T=nu_T_0;
z_TO=nu_TO_0;


b_1=X_X_X(1,1);
b_2=X_X_X(2,1);
b_3=X_X_X(3,1);
b_4=X_X_X(4,1);
b_FD=X_X_X(5,1);
b_FF=X_X_X(6,1);
b_GD=X_X_X(7,1);
b_GF=X_X_X(8,1);
b_HD=X_X_X(9,1);
b_HF=X_X_X(10,1);
c=X_X_X(11,1);
d_1=X_X_X(12,1);
d_2=X_X_X(13,1);
d_3=X_X_X(14,1);
d_4=X_X_X(15,1);
d_F=X_X_X(16,1);
f=X_X_X(17,1);
g=X_X_X(18,1);
i=X_X_X(19,1);
i_1=X_X_X(20,1);
i_2=X_X_X(21,1);
i_3=X_X_X(22,1);
i_4=X_X_X(23,1);
k_1=X_X_X(24,1);
k_2=X_X_X(25,1);
k_3=X_X_X(26,1);
k_4=X_X_X(27,1);
l=X_X_X(28,1);
limda=X_X_X(29,1);
limda_1b=X_X_X(30,1);
limda_1p=X_X_X(31,1);
limda_2b=X_X_X(32,1);
limda_2d=X_X_X(33,1);
limda_2p=X_X_X(34,1);
limda_3b=X_X_X(35,1);
limda_3d=X_X_X(36,1);
limda_3p=X_X_X(37,1);
limda_4b=X_X_X(38,1);
limda_4d=X_X_X(39,1);
limda_4p=X_X_X(40,1);
m=X_X_X(41,1);
o_1=X_X_X(42,1);
o_2=X_X_X(43,1);
o_3=X_X_X(44,1);
o_4=X_X_X(45,1);
p=X_X_X(46,1);
p_D=X_X_X(47,1);
p_F=X_X_X(48,1);
p_O=X_X_X(49,1);
p_2=X_X_X(50,1);
p_3=X_X_X(51,1);
p_4=X_X_X(52,1);
r_D=X_X_X(53,1);
r_F=X_X_X(54,1);
s_1=X_X_X(55,1);
s_2=X_X_X(56,1);
s_3=X_X_X(57,1);
s_4=X_X_X(58,1);
s_F=X_X_X(59,1);
t_L=X_X_X(60,1);
t_O=X_X_X(61,1);
w=X_X_X(62,1);
w_1=X_X_X(63,1);
w_2=X_X_X(64,1);
w_3=X_X_X(65,1);
w_4=X_X_X(66,1);
x_F1=X_X_X(67,1);
x_F2=X_X_X(68,1);
x_F3=X_X_X(69,1);
x_F4=X_X_X(70,1);
x_FF=X_X_X(71,1);
x_H1=X_X_X(72,1);
x_H2=X_X_X(73,1);
x_H3=X_X_X(74,1);
x_H4=X_X_X(75,1);
x_HF=X_X_X(76,1);
y=X_X_X(77,1);
y_2=X_X_X(78,1);
y_3=X_X_X(79,1);
y_4=X_X_X(80,1);
y_F=X_X_X(81,1);


U(1,1)=exp(-c)*(1+exp(z_H)*y^psi-bbb*exp(z_H)*y^psi)-limda;
U(2,1)=exp(z_M-m)+fi_HB*b_HF^2*exp(-3*m)-limda+bbb*limda*exp(p+y);
U(3,1)=(1-bbb*exp(-p-y+r_D))*limda;
U(4,1)=fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f);
U(5,1)=limda*(-exp(s_1)+bbb*exp(s_1)+bbb*exp(d_1));
U(6,1)=limda*(-exp(s_2)+bbb*exp(s_2)+bbb*exp(d_2));
U(7,1)=limda*(-exp(s_3)+bbb*exp(s_3)+bbb*exp(d_3));
U(8,1)=limda*(-exp(s_4)+bbb*exp(s_4)+bbb*exp(d_4));
U(9,1)=-fi_HB*x_HF-limda*exp(s_F)+bbb*limda*(exp(s_F)+exp(d_F));
U(10,1)=-exp(l*(omega+1))+limda*(1-t_L)*exp(w);
U(11,1)=exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p);
U(12,1)=1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D);
U(13,1)=-limda_4b+limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(14,1)=-2*limda_4p*fi_B*b_4/exp(p_4);
U(15,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d;
U(16,1)=-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(17,1)=-limda_4b+limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(18,1)=limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p);
U(19,1)=exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p);
U(20,1)=exp(k_4)-(1-sigma)*exp(k_4-y)-exp(i_4);
U(21,1)=exp(y_4)-exp(-teta*(p_4-p_F)+y_F);
U(22,1)=1-exp(z_4Y-alfa_L4*z_L-alfa_O4*p-y*alfa_L4);
U(23,1)=1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D);
U(24,1)=-limda_3b+limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(25,1)=-2*limda_3p*fi_B*b_3/exp(p_3);
U(26,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d;
U(27,1)=-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(28,1)=-limda_3b+limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(29,1)=limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p);
U(30,1)=exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p);
U(31,1)=exp(k_3)-(1-sigma)*exp(k_3-y)-exp(i_3);
U(32,1)=exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g);
U(33,1)=1-exp(z_3Y-alfa_L3*z_L-alfa_O3*p-y*alfa_L3);
U(34,1)=1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D);
U(35,1)=-limda_2b+limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(36,1)=-2*limda_2p*fi_B*b_2/exp(p_2);
U(37,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d;
U(38,1)=-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(39,1)=-limda_2b+limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(40,1)=limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p);
U(41,1)=exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p);
U(42,1)=exp(k_2)-(1-sigma)*exp(k_2-y)-exp(i_2);
U(43,1)=exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g);
U(44,1)=1-exp(z_2Y-alfa_L2*z_L-alfa_O2*p-y*alfa_L2);
U(45,1)=1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D);
U(46,1)=-limda_1b+limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1);
U(47,1)=-2*limda_1p*fi_B*b_1/exp(p_O);
U(48,1)=-limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+limda_1b;
U(49,1)=-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1);
U(50,1)=exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p);
U(51,1)=exp(k_1)-(1-sigma)*exp(k_1-y)-exp(i_1);
U(52,1)=1-exp(z_1Y-alfa_L1*z_L+alfa_O1*z_O-y*alfa_L1-y*alfa_O1);
U(53,1)=exp(2*p_F+p-f)-exp(z_FP);
U(54,1)=exp(2*y_F+y)-exp(z_FY);
U(55,1)=b_FD*exp(-p_F-y_F)-exp(z_FB);
U(56,1)=x_F1-z_F1;
U(57,1)=x_F2-z_F2;
U(58,1)=x_F3-z_F3;
U(59,1)=x_F4-z_F4;
U(60,1)=exp(d_F-p_F-y_F)-exp(z_FD);
U(61,1)=exp(y+p+2*y_F)-exp(z_FS);
U(62,1)=exp(p-f)-exp(z_FO);
U(63,1)=exp(r_F)-exp(z_FR);
U(64,1)=exp(i)-exp(i_1)-exp(i_2)-exp(i_3)-exp(i_4);
U(65,1)=exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4);
U(66,1)=1-exp(c)-exp(i)-exp(g)-exp(y_F)+exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i));
U(67,1)=1-omega_H*exp((1-teta)*p_D)-(1-omega_H)*exp((1-teta)*p_F);
U(68,1)=exp((1-teta)*p_D)-omega_2*exp((1-teta)*p_2)-(1-omega_2)*exp((1-teta)*p_3);
U(69,1)=x_HF+x_FF-1;
U(70,1)=x_H1+x_F1-1;
U(71,1)=x_H2+x_F2-1;
U(72,1)=x_H3+x_F3-1;
U(73,1)=x_H4+x_F4-1;
U(74,1)=b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4);
U(75,1)=b_GF+b_FF+b_HF;
U(76,1)=exp(p_D*g)+b_GD+b_GF-(t_L*exp(w)+t_O*exp(p_O+o_1)+exp(r_D-p-y)*b_GD+exp(r_F+f-p-y)*b_GF+exp(m)-exp(m-p-y));
U(77,1)=r_D-gama_P*p-gama_Y*y-z_R;
U(78,1)=t_O-z_TO;
U(79,1)=exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T;
U(80,1)=z_B*b_GD-b_GF;
U(81,1)=exp(g)-exp(z_G);


UU=sum(real(U).^2)+10000*sum(imag(U).^2);
return





function UU=sub_ono2(X_X_X)

global alfa_L1 alfa_L2 alfa_L3 alfa_L4 alfa_O1 alfa_O2 alfa_O3 alfa_O4 bbb fi_B fi_D fi_HB fi_P fi_S fi_Y gama_B gama_G gama_P gama_Y nu_1Y_0 nu_2Y_0 nu_3Y_0 nu_4Y_0 nu_B_0 nu_F1_0 nu_F2_0 nu_F3_0 nu_F4_0 nu_FB_0 nu_FD_0 nu_FO_0 nu_FP_0 nu_FR_0 nu_FS_0 nu_FY_0 nu_G_0 nu_H_0 nu_L_0 nu_M_0 nu_O_0 nu_R_0 nu_T_0 nu_TO_0 omega omega_2 omega_H psi sigma teta; 
z_1Y=nu_1Y_0;
z_2Y=nu_2Y_0;
z_3Y=nu_3Y_0;
z_4Y=nu_4Y_0;
z_B=nu_B_0;
z_F1=nu_F1_0;
z_F2=nu_F2_0;
z_F3=nu_F3_0;
z_F4=nu_F4_0;
z_FB=nu_FB_0;
z_FD=nu_FD_0;
z_FO=nu_FO_0;
z_FP=nu_FP_0;
z_FR=nu_FR_0;
z_FS=nu_FS_0;
z_FY=nu_FY_0;
z_G=nu_G_0;
z_H=nu_H_0;
z_L=nu_L_0;
z_M=nu_M_0;
z_O=nu_O_0;
z_R=nu_R_0;
z_T=nu_T_0;
z_TO=nu_TO_0;


b_1=X_X_X(1,1);
b_2=X_X_X(2,1);
b_3=X_X_X(3,1);
b_4=X_X_X(4,1);
b_FD=X_X_X(5,1);
b_GD=X_X_X(6,1);
b_GF=X_X_X(7,1);
b_HF=X_X_X(8,1);
c=X_X_X(9,1);
d_F=X_X_X(10,1);
f=X_X_X(11,1);
g=X_X_X(12,1);
i=X_X_X(13,1);
i_1=X_X_X(14,1);
i_4=X_X_X(15,1);
k_1=X_X_X(16,1);
k_2=X_X_X(17,1);
k_3=X_X_X(18,1);
k_4=X_X_X(19,1);
limda=X_X_X(20,1);
limda_1b=X_X_X(21,1);
limda_1p=X_X_X(22,1);
limda_2b=X_X_X(23,1);
limda_2d=X_X_X(24,1);
limda_2p=X_X_X(25,1);
limda_3b=X_X_X(26,1);
limda_3d=X_X_X(27,1);
limda_3p=X_X_X(28,1);
limda_4b=X_X_X(29,1);
limda_4d=X_X_X(30,1);
limda_4p=X_X_X(31,1);
m=X_X_X(32,1);
o_1=X_X_X(33,1);
o_2=X_X_X(34,1);
o_3=X_X_X(35,1);
o_4=X_X_X(36,1);
p=X_X_X(37,1);
p_D=X_X_X(38,1);
p_F=X_X_X(39,1);
p_O=X_X_X(40,1);
p_2=X_X_X(41,1);
p_3=X_X_X(42,1);
p_4=X_X_X(43,1);
r_F=X_X_X(44,1);
s_1=X_X_X(45,1);
s_2=X_X_X(46,1);
s_3=X_X_X(47,1);
s_4=X_X_X(48,1);
s_F=X_X_X(49,1);
t_L=X_X_X(50,1);
t_O=X_X_X(51,1);
w=X_X_X(52,1);
w_2=X_X_X(53,1);
w_4=X_X_X(54,1);
x_F1=X_X_X(55,1);
x_H2=X_X_X(56,1);
x_H3=X_X_X(57,1);
x_H4=X_X_X(58,1);
y=X_X_X(59,1);
y_2=X_X_X(60,1);
y_3=X_X_X(61,1);
y_4=X_X_X(62,1);
y_F=X_X_X(63,1);


b_FF=-b_GF-b_HF;
b_HD=-b_FD*exp(-2*p_F-2*y_F)+exp(-p_F-y_F+z_FB)-b_GD-b_FD-b_1*exp(p_O+o_1)-b_2*exp(p_2+y_2)-b_3*exp(p_3+y_3)-b_4*exp(p_4+y_4);
d_1=log((1-bbb)/bbb)+s_1;
d_2=log((1-bbb)/bbb)+s_2;
d_3=log((1-bbb)/bbb)+s_3;
d_4=log((1-bbb)/bbb)+s_4;
i_2=log(-1/2*(bbb*b_4*exp(-y-p+m+p_4+y_4)+2*bbb*b_2*t_L*exp(w+p_2+y_2)-bbb*b_2*z_TO*exp(-p_O-o_1+p_2+y_2)-bbb*b_4*t_O*exp(-p_O-o_1+p_4+y_4)-bbb*sigma*b_GD*exp(k_2-y)+bbb*b_2*gama_B*b_GF*exp(p_2+y_2)-bbb*b_2*z_T*exp(p_2+y_2)+bbb*b_2*b_GF*exp(-y-p+r_F+f+p_2+y_2)+exp(s_4)*b_GD*bbb+exp(p_4+y_4)*bbb*b_GD-exp(w_4+y_4)*bbb*b_GD+exp(k_4)*bbb*b_GD-exp(k_4-y)*bbb*b_GD-2*exp(i_4)*bbb*b_GD-exp(o_4+y_4)*bbb*b_GD-bbb*b_GD*exp(s_2)+exp(k_4-y)*sigma*bbb*b_GD-2*bbb*b_4*t_L*exp(w+p_4+y_4)-bbb*b_GD*exp(p_2+y_2)-bbb*b_2*exp(p_2+y_2+p_D*g)+bbb*b_GD*exp(o_2+y_2)+bbb*b_2*exp(m+p_2+y_2)-bbb*b_2*b_GF*exp(p_2+y_2)+bbb*b_4*b_GF*exp(p_4+y_4)-bbb*b_4*b_GF*exp(y_4+p_4-y-p+r_F+f)-bbb*b_4*gama_B*b_GD*exp(p_4+y_4)+2*bbb*b_2*t_O*exp(p_2+y_2+p_O+o_1)-bbb*b_4*exp(y_4+p_4+m)-bbb*b_2*exp(-y-p+m+p_2+y_2)-bbb*b_4*gama_B*b_GF*exp(p_4+y_4)+bbb*b_2*gama_B*b_GD*exp(p_2+y_2)+b_GD*exp(s_2)+bbb*b_GD*exp(k_2-y)+bbb*b_GD*exp(w_2+y_2)-2*bbb*b_4*t_O*exp(p_4+y_4+p_O+o_1)-bbb*b_2*gama_G*exp(p_2+y_2+p_D+g)+bbb*b_4*z_TO*exp(-p_O-o_1+p_4+y_4)+bbb*b_4*exp(p_4+y_4+p_D*g)+bbb*b_4*z_T*exp(p_4+y_4)-exp(s_4)*b_GD-bbb*b_GD*exp(k_2)+bbb*b_2*t_O*exp(-p_O-o_1+p_2+y_2)+bbb*b_4*gama_G*exp(p_4+y_4+p_D+g))/bbb/b_GD);
i_3=log(1/2*(3*bbb*b_GD*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_2*exp(m+p_2+y_2+p_O+o_1)+bbb*b_2*exp(m+y+p+p_2+y_2+p_O+o_1)-6*bbb*b_4*t_L*exp(w+y+p+p_4+y_4+p_O+o_1)+3*bbb*sigma*b_GD*exp(k_4+p+p_O+o_1)+b_GD*exp(s_2+y+p+p_O+o_1)-3*b_GD*exp(s_4+y+p+p_O+o_1)-8*bbb*b_GD*exp(i_4+y+p+p_O+o_1)-bbb*b_GD*exp(k_2+y+p+p_O+o_1)-3*bbb*b_GD*exp(k_4+p+p_O+o_1)+2*bbb*b_2*t_O*exp(y+p+p_2+y_2+2*p_O+2*o_1)+2*bbb*b_2*t_L*exp(w+y+p+p_2+y_2+p_O+o_1)+bbb*b_2*b_GF*exp(r_F+f+p_2+y_2+p_O+o_1)-bbb*b_2*gama_G*exp(y+p+p_2+y_2+p_O+o_1+p_D+g)+3*bbb*b_4*z_TO*exp(p_4+y_4+y+p)-6*bbb*b_4*t_O*exp(y+p+p_4+y_4+2*p_O+2*o_1)-3*bbb*b_4*gama_B*b_GF*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_2*z_T*exp(y+p+p_2+y_2+p_O+o_1)-2*bbb*b_GD*exp(i_1+y+p+p_O+o_1)-3*bbb*b_4*t_O*exp(p_4+y_4+y+p)+3*bbb*b_4*gama_G*exp(y+p+p_4+y_4+p_O+o_1+p_D+g)-bbb*b_2*exp(y+p+p_2+y_2+p_O+o_1+p_D*g)+3*bbb*b_GD*exp(s_4+y+p+p_O+o_1)+2*bbb*b_GD*exp(i+y+p+p_O+o_1)-3*bbb*b_4*b_GF*exp(r_F+f+p_4+y_4+p_O+o_1)+bbb*b_2*t_O*exp(p_2+y_2+y+p)+bbb*b_GD*exp(y+p+o_2+y_2+p_O+o_1)+bbb*b_2*gama_B*b_GD*exp(y+p+p_2+y_2+p_O+o_1)+3*bbb*b_4*b_GF*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_2*b_GF*exp(y+p+p_2+y_2+p_O+o_1)-bbb*b_GD*exp(y+p+p_2+y_2+p_O+o_1)+3*bbb*b_4*exp(m+p_4+y_4+p_O+o_1)+3*bbb*b_4*z_T*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_2*z_TO*exp(p_2+y_2+y+p)-3*bbb*b_4*exp(m+y+p+p_4+y_4+p_O+o_1)+3*bbb*b_GD*exp(k_4+y+p+p_O+o_1)-3*bbb*b_4*gama_B*b_GD*exp(y+p+p_4+y_4+p_O+o_1)+bbb*b_GD*exp(k_2+p+p_O+o_1)-bbb*sigma*b_GD*exp(k_2+p+p_O+o_1)-3*bbb*b_GD*exp(y+p+w_4+y_4+p_O+o_1)-3*bbb*b_GD*exp(y+p+o_4+y_4+p_O+o_1)+bbb*b_2*gama_B*b_GF*exp(y+p+p_2+y_2+p_O+o_1)+bbb*b_GD*exp(y+p+w_2+y_2+p_O+o_1)-bbb*b_GD*exp(s_2+y+p+p_O+o_1)+3*bbb*b_4*exp(y+p+p_4+y_4+p_O+o_1+p_D*g))/bbb/b_GD)-y-p-p_O-o_1;
l=log(limda*(1-t_L)*exp(w))/(omega+1);
r_D=log((exp(y+p+p_O+o_1+p_D*g)-2*t_L*exp(w+y+p+p_O+o_1)-2*t_O*exp(2*p_O+2*o_1+y+p)-gama_B*b_GD*exp(y+p+p_O+o_1)-gama_B*b_GF*exp(y+p+p_O+o_1)+gama_G*exp(y+p+p_O+o_1+p_D+g)+z_T*exp(y+p+p_O+o_1)-t_O*exp(y+p)+z_TO*exp(y+p)-b_GF*exp(p_O+o_1+r_F+f)+b_GD*exp(y+p+p_O+o_1)+b_GF*exp(y+p+p_O+o_1)-exp(m+y+p_O+o_1+p)+exp(p_O+o_1+m))/b_GD)-p_O-o_1;
w_1=-2*o_1+log((bbb*b_GD*exp(k_1+y+p+p_O+o_1)-bbb*b_GD*exp(y+p+p_4+y_4+p_O+o_1)+bbb*b_1*exp(y+p+2*p_O+2*o_1+p_D*g)-b_GD*exp(s_1+y+p+p_O+o_1)+bbb*b_1*exp(2*p_O+2*o_1+m)-bbb*b_1*b_GF*exp(2*p_O+2*o_1+r_F+f)+2*bbb*b_4*t_L*exp(w+y+p+p_4+y_4+p_O+o_1)+bbb*b_1*b_GF*exp(2*p_O+2*o_1+y+p)-bbb*sigma*b_GD*exp(k_4+p+p_O+o_1)+bbb*sigma*b_GD*exp(k_1+p+p_O+o_1)+bbb*b_GD*exp(s_1+y+p+p_O+o_1)+b_GD*exp(s_4+y+p+p_O+o_1)+bbb*b_GD*exp(2*p_O+2*o_1+y+p)+2*bbb*b_GD*exp(i_4+y+p+p_O+o_1)+bbb*b_GD*exp(k_4+p+p_O+o_1)-2*bbb*b_1*t_L*exp(w+2*p_O+2*o_1+y+p)+bbb*b_1*z_TO*exp(y+p+p_O+o_1)-bbb*b_4*z_TO*exp(p_4+y_4+y+p)+2*bbb*b_4*t_O*exp(y+p+p_4+y_4+2*p_O+2*o_1)+bbb*b_4*gama_B*b_GF*exp(y+p+p_4+y_4+p_O+o_1)+bbb*b_1*gama_G*exp(y+p+2*p_O+2*o_1+p_D+g)-2*bbb*b_GD*exp(i_1+y+p+p_O+o_1)+bbb*b_4*t_O*exp(p_4+y_4+y+p)-bbb*b_4*gama_G*exp(y+p+p_4+y_4+p_O+o_1+p_D+g)-bbb*b_GD*exp(s_4+y+p+p_O+o_1)+bbb*b_4*b_GF*exp(r_F+f+p_4+y_4+p_O+o_1)-bbb*b_GD*exp(k_1+p+p_O+o_1)-bbb*b_4*b_GF*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_1*gama_B*b_GD*exp(2*p_O+2*o_1+y+p)-2*bbb*b_1*t_O*exp(3*p_O+3*o_1+y+p)-bbb*b_4*exp(m+p_4+y_4+p_O+o_1)-bbb*b_4*z_T*exp(y+p+p_4+y_4+p_O+o_1)+bbb*b_1*z_T*exp(2*p_O+2*o_1+y+p)+bbb*b_4*exp(m+y+p+p_4+y_4+p_O+o_1)-bbb*b_GD*exp(k_4+y+p+p_O+o_1)+bbb*b_4*gama_B*b_GD*exp(y+p+p_4+y_4+p_O+o_1)+bbb*b_GD*exp(y+p+w_4+y_4+p_O+o_1)-bbb*b_1*exp(2*p_O+2*o_1+m+y+p)+bbb*b_GD*exp(y+p+o_4+y_4+p_O+o_1)-bbb*b_4*exp(y+p+p_4+y_4+p_O+o_1+p_D*g)-bbb*b_1*gama_B*b_GF*exp(2*p_O+2*o_1+y+p)-bbb*b_1*t_O*exp(y+p+p_O+o_1))/bbb/b_GD)-y-p-p_O;
w_3=log(-(4*bbb*b_GD*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_3*exp(m+p_3+y_3+p_O+o_1)-bbb*b_2*exp(m+p_2+y_2+p_O+o_1)+bbb*b_2*exp(m+y+p+p_2+y_2+p_O+o_1)-bbb*b_GD*exp(s_3+y+p+p_O+o_1)-bbb*sigma*b_GD*exp(k_3+p+p_O+o_1)-bbb*b_3*z_TO*exp(p_3+y_3+y+p)-8*bbb*b_4*t_L*exp(w+y+p+p_4+y_4+p_O+o_1)+bbb*b_3*gama_B*b_GF*exp(y+p+p_3+y_3+p_O+o_1)+4*bbb*sigma*b_GD*exp(k_4+p+p_O+o_1)+b_GD*exp(s_2+y+p+p_O+o_1)+bbb*b_3*exp(m+y+p+p_3+y_3+p_O+o_1)+bbb*b_3*gama_B*b_GD*exp(y+p+p_3+y_3+p_O+o_1)+bbb*b_3*t_O*exp(p_3+y_3+y+p)-4*b_GD*exp(s_4+y+p+p_O+o_1)-10*bbb*b_GD*exp(i_4+y+p+p_O+o_1)-bbb*b_GD*exp(k_2+y+p+p_O+o_1)-4*bbb*b_GD*exp(k_4+p+p_O+o_1)+2*bbb*b_2*t_O*exp(y+p+p_2+y_2+2*p_O+2*o_1)+2*bbb*b_2*t_L*exp(w+y+p+p_2+y_2+p_O+o_1)+2*bbb*b_3*t_L*exp(w+y+p+p_3+y_3+p_O+o_1)+bbb*b_2*b_GF*exp(r_F+f+p_2+y_2+p_O+o_1)-bbb*b_2*gama_G*exp(y+p+p_2+y_2+p_O+o_1+p_D+g)+bbb*b_GD*exp(y+p+o_3+y_3+p_O+o_1)+4*bbb*b_4*z_TO*exp(p_4+y_4+y+p)-8*bbb*b_4*t_O*exp(y+p+p_4+y_4+2*p_O+2*o_1)-4*bbb*b_4*gama_B*b_GF*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_3*z_T*exp(y+p+p_3+y_3+p_O+o_1)-bbb*b_2*z_T*exp(y+p+p_2+y_2+p_O+o_1)+b_GD*exp(s_3+y+p+p_O+o_1)-2*bbb*b_GD*exp(i_1+y+p+p_O+o_1)-4*bbb*b_4*t_O*exp(p_4+y_4+y+p)+4*bbb*b_4*gama_G*exp(y+p+p_4+y_4+p_O+o_1+p_D+g)-bbb*b_2*exp(y+p+p_2+y_2+p_O+o_1+p_D*g)+4*bbb*b_GD*exp(s_4+y+p+p_O+o_1)+2*bbb*b_GD*exp(i+y+p+p_O+o_1)-4*bbb*b_4*b_GF*exp(r_F+f+p_4+y_4+p_O+o_1)+bbb*b_2*t_O*exp(p_2+y_2+y+p)-bbb*b_GD*exp(k_3+y+p+p_O+o_1)+bbb*b_GD*exp(y+p+o_2+y_2+p_O+o_1)-bbb*b_GD*exp(y+p+p_3+y_3+p_O+o_1)+bbb*b_2*gama_B*b_GD*exp(y+p+p_2+y_2+p_O+o_1)+4*bbb*b_4*b_GF*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_2*b_GF*exp(y+p+p_2+y_2+p_O+o_1)-bbb*b_3*b_GF*exp(y+p+p_3+y_3+p_O+o_1)-bbb*b_GD*exp(y+p+p_2+y_2+p_O+o_1)+bbb*b_GD*exp(k_3+p+p_O+o_1)+4*bbb*b_4*exp(m+p_4+y_4+p_O+o_1)+4*bbb*b_4*z_T*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_2*z_TO*exp(p_2+y_2+y+p)-4*bbb*b_4*exp(m+y+p+p_4+y_4+p_O+o_1)+4*bbb*b_GD*exp(k_4+y+p+p_O+o_1)-4*bbb*b_4*gama_B*b_GD*exp(y+p+p_4+y_4+p_O+o_1)-bbb*b_3*gama_G*exp(y+p+p_3+y_3+p_O+o_1+p_D+g)+bbb*b_GD*exp(k_2+p+p_O+o_1)-bbb*sigma*b_GD*exp(k_2+p+p_O+o_1)-4*bbb*b_GD*exp(y+p+w_4+y_4+p_O+o_1)-4*bbb*b_GD*exp(y+p+o_4+y_4+p_O+o_1)+bbb*b_2*gama_B*b_GF*exp(y+p+p_2+y_2+p_O+o_1)+bbb*b_GD*exp(y+p+w_2+y_2+p_O+o_1)+2*bbb*b_3*t_O*exp(y+p+p_3+y_3+2*p_O+2*o_1)-bbb*b_GD*exp(s_2+y+p+p_O+o_1)+bbb*b_3*b_GF*exp(r_F+f+p_3+y_3+p_O+o_1)+4*bbb*b_4*exp(y+p+p_4+y_4+p_O+o_1+p_D*g)-bbb*b_3*exp(y+p+p_3+y_3+p_O+o_1+p_D*g))/bbb/b_GD)-y-p-p_O-o_1-y_3;
x_F2=1/2*z_F2-1/2*x_H2+1/2;
x_F3=1/2*z_F3-1/2*x_H3+1/2;
x_F4=1/2*z_F4-1/2*x_H4+1/2;
x_FF=(limda*exp(s_F)-bbb*limda*exp(s_F)-bbb*limda*exp(d_F)+fi_HB)/fi_HB;
x_H1=-2*x_F1+z_F1+1;
x_HF=limda*(-exp(s_F)+bbb*exp(s_F)+bbb*exp(d_F))/fi_HB;


U(1,1)=exp(-c)*(1+exp(z_H)*y^psi-bbb*exp(z_H)*y^psi)-limda;
U(2,1)=exp(z_M-m)+fi_HB*b_HF^2*exp(-3*m)-limda+bbb*limda*exp(p+y);
U(3,1)=(1-bbb*exp(-p-y+r_D))*limda;
U(4,1)=fi_HB*b_HF*exp(-2*m)-limda-bbb*limda*exp(-p-y+r_F+f);
U(5,1)=limda*(-exp(s_1)+bbb*exp(s_1)+bbb*exp(d_1));
U(6,1)=limda*(-exp(s_2)+bbb*exp(s_2)+bbb*exp(d_2));
U(7,1)=limda*(-exp(s_3)+bbb*exp(s_3)+bbb*exp(d_3));
U(8,1)=limda*(-exp(s_4)+bbb*exp(s_4)+bbb*exp(d_4));
U(9,1)=-fi_HB*x_HF-limda*exp(s_F)+bbb*limda*(exp(s_F)+exp(d_F));
U(10,1)=-exp(l*(omega+1))+limda*(1-t_L)*exp(w);
U(11,1)=exp(c)+b_HD+b_HF+exp(m)-exp(w)+exp(w)*t_L-b_HD*exp(-p-y+r_D)-b_HF*exp(-p-y+r_F+f)-x_H1*exp(d_1)-x_H2*exp(d_2)-x_H3*exp(d_3)-x_H4*exp(d_4)-x_HF*exp(d_F)-exp(m-y-p);
U(12,1)=1-limda_4b-2*fi_D*(p+y)*exp(-d_4+y_4)*limda_4p+2*limda_4p*fi_D*(p+y)*exp(-d_4+y_4+y+p-r_D);
U(13,1)=-limda_4b+limda_4p*alfa_L4*exp(-w_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(14,1)=-2*limda_4p*fi_B*b_4/exp(p_4);
U(15,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_4^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_4p+2*limda_4p*fi_Y*y*exp(y+p-r_D)+limda_4b*exp(p_4)-limda_4d;
U(16,1)=-limda_4b+limda_4b*exp(p-r_D)*(1-sigma)+limda_4p*(1-alfa_L4-alfa_O4)*exp(p-r_D+y+y_4-k_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(17,1)=-limda_4b+limda_4p*alfa_O4*exp(-o_4)*(fi_P*p^2+fi_Y*y^2+fi_B*b_4^2+fi_D*(p+y)^2+1);
U(18,1)=limda_4b+limda_4d*teta*exp(-p_4)-limda_4p*exp(-p_4)*(2*fi_P*p-2*fi_B*b_4^2)-limda_4p*fi_P*p*exp(-r_D+y-p_4+p);
U(19,1)=exp(d_4)+exp(i_4)+exp(o_4+y_4)+b_4*exp(p_4+y_4)+exp(w_4+y_4)-exp(p_4+y_4)-b_4*exp(r_D+y_4+p_4-y-p);
U(20,1)=exp(k_4)-(1-sigma)*exp(k_4-y)-exp(i_4);
U(21,1)=exp(y_4)-exp(-teta*(p_4-p_F)+y_F);
U(22,1)=1-exp(z_4Y-alfa_L4*z_L-alfa_O4*p-y*alfa_L4);
U(23,1)=1-limda_3b-2*fi_D*(p+y)*exp(-d_3+y_3)*limda_3p+2*limda_3p*fi_D*(p+y)*exp(-d_3+y_3+y+p-r_D);
U(24,1)=-limda_3b+limda_3p*alfa_L3*exp(-w_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(25,1)=-2*limda_3p*fi_B*b_3/exp(p_3);
U(26,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_3^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_3p+2*limda_3p*fi_Y*y*exp(y+p-r_D)+limda_3b*exp(p_3)-limda_3d;
U(27,1)=-limda_3b+limda_3b*exp(p-r_D)*(1-sigma)+limda_3p*(1-alfa_L3-alfa_O3)*exp(p-r_D+y+y_3-k_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(28,1)=-limda_3b+limda_3p*alfa_O3*exp(-o_3)*(fi_P*p^2+fi_Y*y^2+fi_B*b_3^2+fi_D*(p+y)^2+1);
U(29,1)=limda_3b+limda_3d*teta*exp(-p_3)-limda_3p*exp(-p_3)*(2*fi_P*p-2*fi_B*b_3^2)-limda_3p*fi_P*p*exp(-r_D+y-p_3+p);
U(30,1)=exp(d_3)+exp(i_3)+exp(o_3+y_3)+b_3*exp(p_3+y_3)+exp(w_3+y_3)-exp(p_3+y_3)-b_3*exp(r_D+y_3+p_3-y-p);
U(31,1)=exp(k_3)-(1-sigma)*exp(k_3-y)-exp(i_3);
U(32,1)=exp(y_3)-exp(-teta*p_3)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_3-p_F)+y_F)-exp(-teta*(p_3-p_D)+g);
U(33,1)=1-exp(z_3Y-alfa_L3*z_L-alfa_O3*p-y*alfa_L3);
U(34,1)=1-limda_2b-2*fi_D*(p+y)*exp(-d_2+y_2)*limda_2p+2*limda_2p*fi_D*(p+y)*exp(-d_2+y_2+y+p-r_D);
U(35,1)=-limda_2b+limda_2p*alfa_L2*exp(-w_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(36,1)=-2*limda_2p*fi_B*b_2/exp(p_2);
U(37,1)=(-fi_P*p^2-fi_Y*y^2+fi_B*b_2^2-fi_D*(p+y)^2-1-2*fi_Y*y)*limda_2p+2*limda_2p*fi_Y*y*exp(y+p-r_D)+limda_2b*exp(p_2)-limda_2d;
U(38,1)=-limda_2b+limda_2b*exp(p-r_D)*(1-sigma)+limda_2p*(1-alfa_L2-alfa_O2)*exp(p-r_D+y+y_2-k_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(39,1)=-limda_2b+limda_2p*alfa_O2*exp(-o_2)*(fi_P*p^2+fi_Y*y^2+fi_B*b_2^2+fi_D*(p+y)^2+1);
U(40,1)=limda_2b+limda_2d*teta*exp(-p_2)-limda_2p*exp(-p_2)*(2*fi_P*p-2*fi_B*b_2^2)-limda_2p*fi_P*p*exp(-r_D+y-p_2+p);
U(41,1)=exp(d_2)+exp(i_2)+exp(o_2+y_2)+b_2*exp(p_2+y_2)+exp(w_2+y_2)-exp(p_2+y_2)-b_2*exp(r_D+y_2+p_2-y-p);
U(42,1)=exp(k_2)-(1-sigma)*exp(k_2-y)-exp(i_2);
U(43,1)=exp(y_2)-exp(-teta*p_2)*omega_H*(exp(c)+exp(i))-exp(-teta*(p_2-p_F)+y_F)-exp(-teta*(p_2-p_D)+g);
U(44,1)=1-exp(z_2Y-alfa_L2*z_L-alfa_O2*p-y*alfa_L2);
U(45,1)=1-limda_1b-2*fi_D*(p+y)*exp(-d_1+o_1)*limda_1p+2*limda_1p*fi_D*(p+y)*exp(-d_1+o_1+y+p-r_D);
U(46,1)=-limda_1b+limda_1p*alfa_L1*exp(-w_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1);
U(47,1)=-2*limda_1p*fi_B*b_1/exp(p_O);
U(48,1)=-limda_1p*(-fi_B*b_1^2+fi_D*(p+y)^2+1)*exp(-p_O)+limda_1b;
U(49,1)=-limda_1b+limda_1b*exp(p-r_D)*(1-sigma)+limda_1p*(1-alfa_L1-alfa_O1)*exp(p-r_D+y+o_1-k_1)*(fi_B*b_1^2+fi_D*(p+y)^2+1);
U(50,1)=exp(d_1)+exp(i_1)+b_1*exp(p_O+o_1)+exp(w_1+o_1)-exp(p_O+o_1)-b_1*exp(r_D+o_1+p_O-y-p);
U(51,1)=exp(k_1)-(1-sigma)*exp(k_1-y)-exp(i_1);
U(52,1)=1-exp(z_1Y-alfa_L1*z_L+alfa_O1*z_O-y*alfa_L1-y*alfa_O1);
U(53,1)=exp(2*p_F+p-f)-exp(z_FP);
U(54,1)=exp(2*y_F+y)-exp(z_FY);
U(55,1)=b_FD*exp(-p_F-y_F)-exp(z_FB);
U(56,1)=x_F1-z_F1;
U(57,1)=x_F2-z_F2;
U(58,1)=x_F3-z_F3;
U(59,1)=x_F4-z_F4;
U(60,1)=exp(d_F-p_F-y_F)-exp(z_FD);
U(61,1)=exp(y+p+2*y_F)-exp(z_FS);
U(62,1)=exp(p-f)-exp(z_FO);
U(63,1)=exp(r_F)-exp(z_FR);
U(64,1)=exp(i)-exp(i_1)-exp(i_2)-exp(i_3)-exp(i_4);
U(65,1)=exp(w)-exp(w_1+o_1)-exp(w_2+y_2)-exp(w_3+y_3)-exp(w_4+y_4);
U(66,1)=1-exp(c)-exp(i)-exp(g)-exp(y_F)+exp((1-teta)*p_F)*(1-omega_H)*(exp(c)+exp(i));
U(67,1)=1-omega_H*exp((1-teta)*p_D)-(1-omega_H)*exp((1-teta)*p_F);
U(68,1)=exp((1-teta)*p_D)-omega_2*exp((1-teta)*p_2)-(1-omega_2)*exp((1-teta)*p_3);
U(69,1)=x_HF+x_FF-1;
U(70,1)=x_H1+x_F1-1;
U(71,1)=x_H2+x_F2-1;
U(72,1)=x_H3+x_F3-1;
U(73,1)=x_H4+x_F4-1;
U(74,1)=b_GD+b_FD+b_HD+b_1*exp(p_O+o_1)+b_2*exp(p_2+y_2)+b_3*exp(p_3+y_3)+b_4*exp(p_4+y_4);
U(75,1)=b_GF+b_FF+b_HF;
U(76,1)=exp(p_D*g)+b_GD+b_GF-(t_L*exp(w)+t_O*exp(p_O+o_1)+exp(r_D-p-y)*b_GD+exp(r_F+f-p-y)*b_GF+exp(m)-exp(m-p-y));
U(77,1)=r_D-gama_P*p-gama_Y*y-z_R;
U(78,1)=t_O-z_TO;
U(79,1)=exp(w)*t_L+t_O*exp(p_O+o_1)+gama_B*(b_GD+b_GF)-gama_G*exp(p_D+g)-z_T;
U(80,1)=z_B*b_GD-b_GF;
U(81,1)=exp(g)-exp(z_G);


UU=sum(real(U).^2)+10000*sum(imag(U).^2);
return