function [ys,check,fval]=model_D7_steadystate(junk,ys)
check=0;
fval=0;
global M_
%keyboard;
if 1==0
    ys=real(ys);
    ys=fminsearch(@ono_ER,ys,optimset('Display','off','MaxFunEvals',350,'MaxIter',20000),[M_.fname '_static'],M_.params,[zeros(34,1)]);
    %X_X=[b_1_F,b_1_G,b_1_H_1,b_1_H_2,b_1_H_3,b_2_F,b_2_G,b_2_H_1,b_2_H_2,b_2_H_3,b_3_F,b_3_G,b_3_H_1,b_3_H_2,b_3_H_3,c_1,c_2,c_3,d_1,d_2,d_3,f_1,f_2,f_3,g_1,g_2,g_3,i_1,i_2,i_3,k_1,k_2,k_3,l_1,l_2,l_3,limda_1,limda_1_b,limda_1_d,limda_1_k,limda_1_p,limda_2,limda_2_b,limda_2_d,limda_2_k,limda_2_p,limda_3,limda_3_b,limda_3_d,limda_3_k,limda_3_p,m_1,m_1_F,m_1_H,m_2,m_2_F,m_2_H,m_3,m_3_F,m_3_H,p_1,p_1_P,p_2,p_2_P,p_3,p_3_P,r_1,r_2,r_3,s_1,s_2,s_3,tax_1_L,tax_1_Y,tax_2_L,tax_2_Y,tax_3_L,tax_3_Y,w_1,w_2,w_3,x_1_H_1,x_1_H_2,x_1_H_3,x_2_H_1,x_2_H_2,x_2_H_3,x_3_H_1,x_3_H_2,x_3_H_3,y_1,y_1_D,y_1_P,y_2,y_2_D,y_2_P,y_3,y_3_D,y_3_P].';
    %X_X=real(X_X);
    %X_X=ys0(1:99,1);
    X_X=ones(68,1);
    
    %X_X=fminsearch(@ono_eq,X_X,optimset('Display','final','MaxFunEvals',7350,'MaxIter',20000));

    %ys=ono_chast(X_X);
    %X_X=fminunc(@ono_chast,X_X,optimset('Display','off','MaxFunEvals',13,'MaxIter',20000,'GradObj','on','Hessian','off'));
    X_X=fsolve(@ono_eq,X_X,optimset('Display','off','MaxFunEvals',20000,'MaxIter',125,'Jacobian','on','NonlEqnAlgorithm','lm'));
    %ys=[X_X.',z_1_bbb,z_1_BF,z_1_EC,z_1_G,z_1_L,z_1_MF,z_1_MH,z_1_R,z_1_tL,z_1_tY,z_1_Y,z_2_bbb,z_2_BF,z_2_EC,z_2_G,z_2_L,z_2_MF,z_2_MH,z_2_R,z_2_tL,z_2_tY,z_2_Y,z_3_bbb,z_3_BF,z_3_G,z_3_L,z_3_MF,z_3_MH,z_3_R,z_3_tL,z_3_tY,z_3_Y,z_EP,z_Y,zzobs_1_C,zzobs_1_Ex,zzobs_1_G,zzobs_1_I,zzobs_1_Im,zzobs_1_M,zzobs_1_P,zzobs_1_R,zzobs_1_S,zzobs_1_Y,zzobs_2_C,zzobs_2_Ex,zzobs_2_G,zzobs_2_I,zzobs_2_Im,zzobs_2_M,zzobs_2_P,zzobs_2_R,zzobs_2_S,zzobs_2_W,zzobs_2_WL,zzobs_2_Y,zzobs_oil,zzobs_1_BF,zzobs_1_BH,zzobs_1_fr,zzobs_2_BF,zzobs_2_BH,zzobs_2_fr,zzobs_2_fn,zzobs_12_fn,zzobs_12_fr].';
    ys=ono_ys_chast(X_X);
end

for i=1:M_.param_nbr
    eval([deblank(M_.param_names(i,:)),'=M_.params(',num2str(i),');']);   
end

omega_1_1=1-omega_2_1;
omega_2_2=1-omega_1_2;
nu0_1_bbb=nu0_bbb;
nu0_2_bbb=nu0_bbb;
teta_1=teta;
teta_2=teta;

alfa_1=nu0_1_alfa;
alfa_2=nu0_2_alfa;

z_1_alfa=nu0_1_alfa;
z_2_alfa=nu0_2_alfa;

z_1_bbb=nu0_1_bbb;
z_1_BF=nu0_1_BF;
z_1_BK=nu0_1_BK;
z_1_DB=nu0_1_DB;
z_1_DBM=nu0_1_DBM;
z_1_G=nu0_1_G;
z_1_L=nu0_1_L;
z_1_MH=nu0_1_MH;
z_1_R=nu0_1_R;
z_1_SIG=nu0_1_SIG;
z_1_tL=nu0_1_tL;
z_1_tr=nu0_1_tr;
z_1_tY=nu0_1_tY;
z_1_Y=nu0_1_Y;
z_2_bbb=nu0_2_bbb;
z_2_BF=nu0_2_BF;
z_2_BK=nu0_2_BK;
z_2_DB=nu0_2_DB;
z_2_DBM=nu0_2_DBM;
z_2_G=nu0_2_G;
z_2_L=nu0_2_L;
z_2_MH=nu0_2_MH;
z_2_R=nu0_2_R;
z_2_SIG=nu0_2_SIG;
z_2_tL=nu0_2_tL;
z_2_tr=nu0_2_tr;
z_2_tY=nu0_2_tY;
z_2_Y=nu0_2_Y;
z_Y=nu0_Y;

k_1=k_1_SS;
k_2=k_2_SS;
r_2_F=r_2_F_SS;
p_1=p_1_SS;
c_1=c_1_SS;
b_2_B_1=b_2_B_1_SS;
r_1_F=r_1_F_SS;
c_2=c_2_SS;
f_1=f_1_SS;

y_1=nu0_Y;
y_2=y_1;


k_1=k_1_SS;
k_2=k_2_SS;
r_2_F=r_2_F_SS;
p_1=p_1_SS;
c_1=c_1_SS;
b_2_B_1=b_2_B_1_SS;
r_1_F=r_1_F_SS;
c_2=c_2_SS;
f_1=f_1_SS;

if 1==2
p_2=-nu0_1_R+p_1+nu0_2_R;
y_2_D=log(exp(c_2)+exp(k_2)-exp(k_2-y_2)+exp(k_2-y_2)*nu0_2_SIG+exp(nu0_2_G));
limda_2_p=-(exp(-2*p_2-2*y_2-2*nu0_2_R+4*r_2_F-k_2)*nu0_2_BF^2-exp(-2*p_2-2*y_2-nu0_2_R+3*r_2_F-k_2)*nu0_2_BF^2-exp(p_2-2*nu0_2_R+k_2+r_2_F)+exp(p_2-3*nu0_2_R+k_2+2*r_2_F)*nu0_2_DB-2*exp(-3*nu0_2_R+3*r_2_F-y_2)*nu0_2_BF+exp(p_2-2*nu0_2_R+k_2+r_2_F)*nu0_2_SIG-exp(p_2-3*nu0_2_R+k_2+2*r_2_F)*nu0_2_SIG*nu0_2_DB+2*exp(-3*nu0_2_R+3*r_2_F-y_2)*nu0_2_SIG*nu0_2_BF+exp(k_2-nu0_2_R+r_2_F)-exp(k_2-2*nu0_2_R+2*r_2_F)*nu0_2_DB+2*exp(-p_2-y_2-2*nu0_2_R+3*r_2_F)*nu0_2_BF)*exp(nu0_2_R-y_2-p_2)/(exp(-nu0_2_R+r_2_F)-exp(-2*nu0_2_R+2*r_2_F)*nu0_2_DB+2*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)*nu0_2_BF)/(-1+alfa_2);
limda_2=-exp(c_2*(1+omega_2_C))/(exp(c_2+nu0_2_MH)-exp(c_2+nu0_2_MH-nu0_bbb+p_2+y_2)-exp(c_2));
y_1_D=log(exp(c_1)+exp(k_1)-exp(k_1-y_2)+exp(k_1-y_2)*nu0_1_SIG+exp(nu0_1_G));
limda_1_p=(nu0_1_BF^2*exp(-2*nu0_1_R-2*y_2-2*p_1+4*r_1_F-k_1)-exp(-nu0_1_R-2*y_2-2*p_1+3*r_1_F-k_1)*nu0_1_BF^2-exp(-2*nu0_1_R+p_1+k_1+r_1_F)+exp(-3*nu0_1_R+p_1+k_1+2*r_1_F)*nu0_1_DB-2*exp(-3*nu0_1_R+3*r_1_F-y_2)*nu0_1_BF+exp(-2*nu0_1_R+p_1+k_1+r_1_F)*nu0_1_SIG-exp(-3*nu0_1_R+p_1+k_1+2*r_1_F)*nu0_1_SIG*nu0_1_DB+2*exp(-3*nu0_1_R+3*r_1_F-y_2)*nu0_1_SIG*nu0_1_BF+exp(k_1-nu0_1_R+r_1_F)-exp(k_1-2*nu0_1_R+2*r_1_F)*nu0_1_DB+2*exp(-2*nu0_1_R-y_2-p_1+3*r_1_F)*nu0_1_BF)*exp(nu0_1_R-y_2-p_1)/(-1+alfa_1)/(-exp(-nu0_1_R+r_1_F)+exp(-2*nu0_1_R+2*r_1_F)*nu0_1_DB-2*exp(-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*nu0_1_BF);
limda_1=exp(c_1*(1+omega_1_C))/(-exp(c_1+nu0_1_MH)+exp(c_1+nu0_1_MH-nu0_bbb+p_1+y_2)+exp(c_1));
b_1_B_2=.50000000000000000000000000000000*(nu0_1_BK*fi_1_HB*limda_1*nu0_1_tY*exp(p_1+y_2)-2.*exp(-1.*f_1)*nu0_1_BK*fi_1_HB*b_2_B_1*fi_1_HM*nu0_1_tY*teta+2.*nu0_1_BK*fi_1_HB*nu0_1_tr*fi_1_HM*teta*nu0_1_tY+nu0_1_BK*fi_1_HB*limda_1*exp(-1.*nu0_bbb+2.*p_1+2.*y_2)-1.*nu0_1_BK*fi_1_HB*limda_1*nu0_1_tY+2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_BK+2.*nu0_1_BF*fi_1_HM*fi_1_HB*teta+2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_tY+2.*nu0_1_BK*fi_1_HB*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*teta-2.*nu0_1_BK*fi_1_HB*fi_1_HM*teta*exp(nu0_1_G)-1.*nu0_1_BK*fi_1_HB*limda_1*nu0_1_tY*exp(-1.*nu0_bbb+2.*p_1+2.*y_2)-1.*nu0_1_BK*fi_1_HB*limda_1*teta-2.*nu0_1_BK*fi_1_HB*nu0_1_tr*fi_1_HM*teta+2.*nu0_1_BK*fi_1_HM*limda_1*exp(nu0_bbb+nu0_1_R-1.*p_1-1.*y_2)*teta+2.*nu0_1_BK*fi_1_HM*limda_1*exp(nu0_bbb+nu0_1_R-1.*p_1-1.*y_2)*nu0_1_tY-2.*nu0_1_BK*fi_1_HM*limda_1*exp(nu0_bbb+nu0_1_R-1.*p_1-1.*y_2)*nu0_1_tY*teta-1.*nu0_1_BK*fi_1_HB*limda_1*teta*exp(-1.*nu0_bbb+2.*p_1+2.*y_2)-2.*nu0_1_BK*fi_1_HB*nu0_1_tY*teta*fi_1_HM*exp(p_1+y_2+c_1+nu0_1_MH)-1.*nu0_1_BK*fi_1_HB*limda_1*teta*nu0_1_tY*exp(p_1+y_2)+nu0_1_BK*fi_1_HB*limda_1*teta*nu0_1_tY*exp(-1.*nu0_bbb+2.*p_1+2.*y_2)+nu0_1_BK*fi_1_HB*limda_1-2.*nu0_1_BK*fi_1_HB*nu0_1_tr*fi_1_HM*nu0_1_tY+2.*nu0_1_BK*fi_1_HB*fi_1_HM*exp(c_1+nu0_1_MH)-2.*nu0_1_BK*fi_1_HB*b_2_B_1*fi_1_HM*teta*exp(-1.*f_1+nu0_1_R-1.*y_2-1.*p_1)+nu0_1_BK*fi_1_HB*limda_1*teta*nu0_1_tY+2.*exp(-1.*f_1)*nu0_1_BK*fi_1_HB*b_2_B_1*fi_1_HM*nu0_1_tY+2.*nu0_1_BK*fi_1_HB*nu0_1_tY*teta*fi_1_HM*exp(nu0_1_G)+2.*exp(-1.*f_1)*nu0_1_BK*fi_1_HB*b_2_B_1*fi_1_HM*teta+2.*nu0_1_BK*fi_1_HB*nu0_1_tY*teta*fi_1_HM*exp(c_1+nu0_1_MH)+2.*nu0_1_BK*fi_1_HB*fi_1_HM*exp(nu0_1_G)-2.*nu0_1_BK*fi_1_HM*fi_1_HB*mu_B_1-2.*nu0_1_BK*fi_1_HM*limda_1*teta-2.*nu0_1_BK*fi_1_HM*limda_1*nu0_1_tY-1.*nu0_1_BK*fi_1_HB*limda_1*exp(p_1+y_2)+nu0_1_BK*fi_1_HB*limda_1*teta*exp(-1.*nu0_bbb+p_1+y_2)-2.*nu0_1_BK*fi_1_HB*fi_1_HM*exp(p_1+y_2+c_1+nu0_1_MH)-2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_BK*exp(nu0_1_R-1.*y_2-1.*p_1)+2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_BK*exp(nu0_1_R-1.*y_2-1.*p_1)*teta-2.*nu0_1_BK*fi_1_HM*fi_1_HB*mu_B_1*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY-2.*nu0_1_BK*fi_1_HB*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*teta*nu0_1_tY+2.*nu0_1_BK*fi_1_HM*fi_1_HB*mu_B_1*exp(nu0_1_R-1.*y_2-1.*p_1)-2.*nu0_1_BK*fi_1_HM*fi_1_HB*mu_B_1*exp(nu0_1_R-1.*y_2-1.*p_1)*teta+2.*nu0_1_BK*fi_1_HM*fi_1_HB*mu_B_1*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY*teta-2.*nu0_1_BK*fi_1_HB*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM+2.*nu0_1_BK*fi_1_HB*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*nu0_1_tY-2.*nu0_1_BK*fi_1_HB*fi_1_HM*nu0_1_tY*exp(c_1+nu0_1_MH)-2.*nu0_1_BF*fi_1_HM*fi_1_HB+2.*nu0_1_BK*fi_1_HB*fi_1_HM*nu0_1_tY*exp(p_1+y_2+c_1+nu0_1_MH)-2.*nu0_1_BK*fi_1_HM*limda_1*exp(nu0_bbb+nu0_1_R-1.*p_1-1.*y_2)+2.*nu0_1_BF*fi_1_HM*fi_1_HB*exp(nu0_1_R-1.*y_2-1.*p_1)-2.*nu0_1_BK*fi_1_HM*limda_1*exp(nu0_1_R-1.*y_2-1.*p_1)+2.*nu0_1_BK*fi_1_HM*limda_1-2.*nu0_1_BK*fi_1_HB*fi_1_HM*nu0_1_tY*exp(nu0_1_G)+nu0_1_BK*fi_1_HB*limda_1*nu0_1_tY*exp(-1.*nu0_bbb+p_1+y_2)-2.*nu0_1_BK*fi_1_HB*b_2_B_1*fi_1_HM*nu0_1_tY*exp(-1.*f_1+nu0_1_R-1.*y_2-1.*p_1)+2.*nu0_1_BK*fi_1_HB*b_2_B_1*fi_1_HM*exp(-1.*f_1+nu0_1_R-1.*y_2-1.*p_1)-2.*nu0_1_BK*fi_1_HB*fi_1_HM*teta*exp(c_1+nu0_1_MH)-2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_BK*teta-2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_BK*nu0_1_tY+2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_BK*nu0_1_tY*teta-2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_tY*teta+2.*nu0_1_BK*fi_1_HB*nu0_1_tr*fi_1_HM-1.*nu0_1_BK*fi_1_HB*limda_1*teta*nu0_1_tY*exp(-1.*nu0_bbb+p_1+y_2)+2.*nu0_1_BK*fi_1_HM*limda_1*exp(nu0_1_R-1.*y_2-1.*p_1)*teta+2.*nu0_1_BK*fi_1_HM*limda_1*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY-2.*nu0_1_BK*fi_1_HM*limda_1*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY*teta+nu0_1_BK*fi_1_HB*limda_1*teta*exp(p_1+y_2)-1.*nu0_1_BK*fi_1_HB*limda_1*exp(-1.*nu0_bbb+p_1+y_2)+2.*nu0_1_BK*fi_1_HM*fi_1_HB*mu_B_1*teta+2.*nu0_1_BK*fi_1_HM*fi_1_HB*mu_B_1*nu0_1_tY-2.*nu0_1_BK*fi_1_HM*fi_1_HB*mu_B_1*nu0_1_tY*teta+2.*nu0_1_BK*fi_1_HM*limda_1*nu0_1_tY*teta+2.*nu0_1_BK*fi_1_HB*b_2_B_1*fi_1_HM*nu0_1_tY*teta*exp(-1.*f_1+nu0_1_R-1.*y_2-1.*p_1)+2.*nu0_1_BK*fi_1_HB*fi_1_HM*teta*exp(p_1+y_2+c_1+nu0_1_MH)+2.*nu0_1_BK*fi_1_HB*nu0_1_tY*teta*limda_1_p*fi_1_HM+2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_BK*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY-2.*nu0_1_BF*fi_1_HM*fi_1_HB*nu0_1_BK*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY*teta-2.*nu0_1_BF*fi_1_HM*fi_1_HB*exp(nu0_1_R-1.*y_2-1.*p_1)*teta-2.*nu0_1_BF*fi_1_HM*fi_1_HB*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY+2.*nu0_1_BF*fi_1_HM*fi_1_HB*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY*teta-2.*exp(-1.*f_1)*nu0_1_BK*fi_1_HB*b_2_B_1*fi_1_HM+2.*nu0_1_BK*fi_1_HM*limda_1*exp(2.*nu0_1_R-2.*y_2-2.*p_1+nu0_bbb)-2.*nu0_1_BK*fi_1_HM*limda_1*teta*exp(2.*nu0_1_R-2.*y_2-2.*p_1+nu0_bbb)+2.*nu0_1_BK*fi_1_HM*limda_1*nu0_1_tY*teta*exp(2.*nu0_1_R-2.*y_2-2.*p_1+nu0_bbb)-2.*nu0_1_BK*fi_1_HM*limda_1*nu0_1_tY*exp(2.*nu0_1_R-2.*y_2-2.*p_1+nu0_bbb))/nu0_1_BK/fi_1_HM/fi_1_HB/(-1.+teta+nu0_1_tY-1.*nu0_1_tY*teta+exp(nu0_1_R-1.*y_2-1.*p_1)-1.*exp(nu0_1_R-1.*y_2-1.*p_1)*teta-1.*exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY+exp(nu0_1_R-1.*y_2-1.*p_1)*nu0_1_tY*teta);
b_2_B_2=1/2*(2*exp(c_2+nu0_2_MH+teta*f_1-y_2_D+y_1_D)*fi_2_HM*omega_2_1+2*nu0_2_tY*fi_2_HM*teta*limda_2_p+2*exp(nu0_2_G+teta*f_1-y_2_D+y_1_D)*fi_2_HM*teta*nu0_2_tY*omega_2_1-2*exp(c_2+nu0_2_MH+teta*f_1-y_2_D+y_1_D)*fi_2_HM*teta*omega_2_1-2*b_1_B_2*exp(f_1+teta*f_1-y_2_D+y_1_D)*fi_2_HM*teta*nu0_2_tY*omega_2_1+2*b_1_B_2*exp(f_1)*fi_2_HM*omega_1_2+limda_2-2*exp(p_2+y_2+c_2+nu0_2_MH+teta*f_1-y_2_D+y_1_D)*fi_2_HM*teta*nu0_2_tY*omega_2_1+2*nu0_2_tY*fi_2_HM*omega_2_1*teta*limda_2_p*exp(teta*f_1-y_2_D+y_1_D)+2*b_1_B_2*exp(f_1+nu0_2_R-y_2-p_2)*fi_2_HM+2*b_1_B_2*exp(f_1+nu0_2_R-y_2-p_2)*fi_2_HM*teta*nu0_2_tY-exp(2*p_2+2*y_2-nu0_bbb+teta*f_1-y_2_D+y_1_D)*limda_2*teta*omega_2_1-2*b_1_B_2*exp(f_1+nu0_2_R-y_2-p_2)*fi_2_HM*teta*nu0_2_tY*omega_1_2-2*b_1_B_2*exp(f_1+teta*f_1-y_2_D+y_1_D+nu0_2_R-y_2-p_2)*fi_2_HM*teta*omega_2_1+2*exp(c_2+nu0_2_MH+teta*f_1-y_2_D+y_1_D)*fi_2_HM*teta*nu0_2_tY*omega_2_1+2*b_1_B_2*exp(f_1+nu0_2_R-y_2-p_2)*fi_2_HM*nu0_2_tY*omega_1_2-exp(2*p_2+2*y_2-nu0_bbb+teta*f_1-y_2_D+y_1_D)*limda_2*nu0_2_tY*omega_2_1+2*b_1_B_2*exp(f_1+nu0_2_R-y_2-p_2)*fi_2_HM*teta*omega_1_2-2*b_1_B_2*exp(f_1+nu0_2_R-y_2-p_2)*fi_2_HM*omega_1_2+2*exp(nu0_2_G)*fi_2_HM*teta*nu0_2_tY+2*exp(c_2+nu0_2_MH)*fi_2_HM*teta*omega_1_2+2*exp(c_2+nu0_2_MH)*fi_2_HM*teta*nu0_2_tY-2*b_1_B_2*exp(f_1)*fi_2_HM-2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*teta*nu0_2_tY*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)-2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)+2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*teta*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)+2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*nu0_2_tY*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)+2*exp(nu0_2_G)*fi_2_HM*teta*omega_1_2+2*exp(nu0_2_G)*fi_2_HM*nu0_2_tY*omega_1_2-exp(p_2+y_2-nu0_bbb+teta*f_1-y_2_D+y_1_D)*limda_2*teta*nu0_2_tY*omega_2_1+limda_2*exp(p_2+y_2+teta*f_1-y_2_D+y_1_D)*nu0_2_tY*omega_2_1+exp(p_2+y_2-nu0_bbb+teta*f_1-y_2_D+y_1_D)*limda_2*nu0_2_tY*omega_2_1+exp(2*p_2+2*y_2-nu0_bbb+teta*f_1-y_2_D+y_1_D)*limda_2*teta*nu0_2_tY*omega_2_1-2*nu0_2_tY*fi_2_HM*omega_1_2*teta*limda_2_p-2*b_1_B_2*exp(f_1)*fi_2_HM*nu0_2_tY*omega_1_2-2*exp(nu0_2_G)*fi_2_HM*teta*nu0_2_tY*omega_1_2-2*b_1_B_2*exp(f_1)*fi_2_HM*teta*omega_1_2-2*b_1_B_2*exp(f_1)*fi_2_HM*teta*nu0_2_tY+2*exp(c_2+nu0_2_MH)*fi_2_HM*nu0_2_tY*omega_1_2+2*b_1_B_2*exp(f_1)*fi_2_HM*teta*nu0_2_tY*omega_1_2-2*exp(c_2+nu0_2_MH)*fi_2_HM*teta*nu0_2_tY*omega_1_2-2*exp(c_2+nu0_2_MH)*fi_2_HM*omega_1_2-2*exp(c_2+nu0_2_MH)*fi_2_HM*teta-2*exp(c_2+nu0_2_MH)*fi_2_HM*nu0_2_tY+exp(p_2+y_2-nu0_bbb)*limda_2*nu0_2_tY-limda_2*teta*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)-limda_2*nu0_2_tY*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)-exp(p_2+y_2-nu0_bbb)*limda_2*teta*omega_1_2-exp(p_2+y_2-nu0_bbb)*limda_2*nu0_2_tY*omega_1_2-exp(p_2+y_2-nu0_bbb)*limda_2*teta*nu0_2_tY+limda_2*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)+2*nu0_2_tr*fi_2_HM+limda_2*teta*nu0_2_tY*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)+exp(p_2+y_2-nu0_bbb)*limda_2*teta+2*nu0_2_tr*fi_2_HM*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)+exp(p_2+y_2-nu0_bbb)*limda_2*omega_1_2+2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*nu0_2_tY-2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*teta*nu0_2_tY+2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*teta*nu0_2_tY*omega_1_2+2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*omega_1_2+2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*teta-2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM+exp(p_2+y_2-nu0_bbb)*limda_2*teta*nu0_2_tY*omega_1_2-2*nu0_2_tr*fi_2_HM*nu0_2_tY*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)+2*nu0_2_tr*fi_2_HM*teta*nu0_2_tY*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)-2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*teta*omega_1_2-2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM*nu0_2_tY*omega_1_2+2*exp(nu0_2_G)*fi_2_HM-exp(p_2+y_2-nu0_bbb+teta*f_1-y_2_D+y_1_D)*limda_2*omega_2_1-2*exp(c_2+nu0_2_MH+teta*f_1-y_2_D+y_1_D)*fi_2_HM*nu0_2_tY*omega_2_1+2*b_1_B_2*exp(f_1+teta*f_1-y_2_D+y_1_D+nu0_2_R-y_2-p_2)*fi_2_HM*omega_2_1-2*exp(p_2+y_2+c_2+nu0_2_MH+teta*f_1-y_2_D+y_1_D)*fi_2_HM*omega_2_1-limda_2*omega_1_2+exp(2*p_2+2*y_2-nu0_bbb+teta*f_1-y_2_D+y_1_D)*limda_2*omega_2_1+limda_2*exp(p_2+y_2+teta*f_1-y_2_D+y_1_D)*teta*omega_2_1-2*b_1_B_2*exp(f_1+nu0_2_R-y_2-p_2)*fi_2_HM*nu0_2_tY+exp(p_2+y_2-nu0_bbb+teta*f_1-y_2_D+y_1_D)*limda_2*teta*omega_2_1-limda_2*exp(p_2+y_2+teta*f_1-y_2_D+y_1_D)*omega_2_1-limda_2*teta*nu0_2_tY*omega_1_2+2*b_1_B_2*exp(f_1+teta*f_1-y_2_D+y_1_D+nu0_2_R-y_2-p_2)*fi_2_HM*teta*nu0_2_tY*omega_2_1+2*b_1_B_2*exp(f_1+teta*f_1-y_2_D+y_1_D)*fi_2_HM*teta*omega_2_1-2*exp(nu0_2_G+teta*f_1-y_2_D+y_1_D)*fi_2_HM*teta*omega_2_1+2*exp(p_2+y_2+c_2+nu0_2_MH+teta*f_1-y_2_D+y_1_D)*fi_2_HM*nu0_2_tY*omega_2_1+2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM*teta+2*exp(nu0_2_G+teta*f_1-y_2_D+y_1_D)*fi_2_HM*omega_2_1-2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM*teta*omega_1_2-2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM*nu0_2_tY*omega_1_2-2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM*teta*nu0_2_tY-limda_2*exp(p_2+y_2)*teta*omega_1_2-limda_2*exp(p_2+y_2)*nu0_2_tY*omega_1_2-limda_2*exp(p_2+y_2)*teta*nu0_2_tY+limda_2*teta*omega_1_2+limda_2*nu0_2_tY*omega_1_2+limda_2*teta*nu0_2_tY-limda_2*teta-limda_2*nu0_2_tY-limda_2*exp(p_2+y_2+teta*f_1-y_2_D+y_1_D)*teta*nu0_2_tY*omega_2_1-2*exp(nu0_2_G+teta*f_1-y_2_D+y_1_D)*fi_2_HM*nu0_2_tY*omega_2_1-2*b_1_B_2*exp(f_1+teta*f_1-y_2_D+y_1_D)*fi_2_HM*omega_2_1+2*b_1_B_2*exp(f_1+teta*f_1-y_2_D+y_1_D)*fi_2_HM*nu0_2_tY*omega_2_1-2*b_1_B_2*exp(f_1+teta*f_1-y_2_D+y_1_D+nu0_2_R-y_2-p_2)*fi_2_HM*nu0_2_tY*omega_2_1+2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM*teta*nu0_2_tY*omega_1_2+2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM*nu0_2_tY-limda_2*exp(p_2+y_2)-2*nu0_2_tr*fi_2_HM*teta*omega_2_1*exp(teta*f_1-y_2_D+y_1_D)+limda_2*exp(p_2+y_2)*omega_1_2-2*b_1_B_2*exp(f_1+nu0_2_R-y_2-p_2)*fi_2_HM*teta+2*nu0_2_tr*fi_2_HM*teta*omega_1_2+2*nu0_2_tr*fi_2_HM*nu0_2_tY*omega_1_2+2*nu0_2_tr*fi_2_HM*teta*nu0_2_tY-2*nu0_2_tr*fi_2_HM*omega_1_2-2*nu0_2_tr*fi_2_HM*teta-2*nu0_2_tr*fi_2_HM*nu0_2_tY+limda_2*exp(p_2+y_2)*teta+limda_2*exp(p_2+y_2)*nu0_2_tY-exp(2*p_2+2*y_2-nu0_bbb)*limda_2*teta-exp(2*p_2+2*y_2-nu0_bbb)*limda_2*omega_1_2-exp(2*p_2+2*y_2-nu0_bbb)*limda_2*nu0_2_tY-2*exp(nu0_2_G)*fi_2_HM*omega_1_2-2*exp(nu0_2_G)*fi_2_HM*teta-2*exp(nu0_2_G)*fi_2_HM*nu0_2_tY-2*nu0_2_tr*fi_2_HM*teta*nu0_2_tY*omega_1_2+2*exp(p_2+y_2+c_2+nu0_2_MH+teta*f_1-y_2_D+y_1_D)*fi_2_HM*teta*omega_2_1+exp(2*p_2+2*y_2-nu0_bbb)*limda_2*teta*omega_1_2+exp(2*p_2+2*y_2-nu0_bbb)*limda_2*nu0_2_tY*omega_1_2+exp(2*p_2+2*y_2-nu0_bbb)*limda_2*teta*nu0_2_tY+2*b_1_B_2*exp(f_1)*fi_2_HM*teta+2*exp(c_2+nu0_2_MH)*fi_2_HM-exp(p_2+y_2-nu0_bbb)*limda_2-2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM+limda_2*exp(p_2+y_2)*teta*nu0_2_tY*omega_1_2+2*b_1_B_2*exp(f_1)*fi_2_HM*nu0_2_tY+exp(2*p_2+2*y_2-nu0_bbb)*limda_2-exp(2*p_2+2*y_2-nu0_bbb)*limda_2*teta*nu0_2_tY*omega_1_2+2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM*omega_1_2)/(-1+nu0_2_tY)/(-1+teta)/(1-exp(nu0_2_R-y_2-p_2)-omega_1_2+omega_1_2*exp(nu0_2_R-y_2-p_2)+omega_2_1*exp(teta*f_1-y_2_D+y_1_D)-omega_2_1*exp(nu0_2_R-y_2-p_2+teta*f_1-y_2_D+y_1_D))/fi_2_HM;
b_2_G=1/2*(2*exp(nu0_2_G)*fi_2_HM+2*exp(c_2+nu0_2_MH)*fi_2_HM+limda_2-exp(p_2+y_2-nu0_bbb)*limda_2+2*nu0_2_tr*fi_2_HM-2*exp(p_2+y_2+c_2+nu0_2_MH)*fi_2_HM-limda_2*exp(p_2+y_2)+exp(2*p_2+2*y_2-nu0_bbb)*limda_2-2*nu0_2_tY*fi_2_HM*exp((y_2_D+log((-1+omega_1_2-omega_2_1*exp(teta*f_1-y_2_D+y_1_D))/teta*(-1+nu0_2_tY)*(-1+teta)/limda_2_p))/(-1+teta)+(-teta*log((-1+omega_1_2-omega_2_1*exp(teta*f_1-y_2_D+y_1_D))/teta*(-1+nu0_2_tY)*(-1+teta)/limda_2_p)-y_2_D)/(-1+teta))+2*nu0_2_tY*fi_2_HM*omega_1_2*exp((y_2_D+log((-1+omega_1_2-omega_2_1*exp(teta*f_1-y_2_D+y_1_D))/teta*(-1+nu0_2_tY)*(-1+teta)/limda_2_p))/(-1+teta)+(-teta*log((-1+omega_1_2-omega_2_1*exp(teta*f_1-y_2_D+y_1_D))/teta*(-1+nu0_2_tY)*(-1+teta)/limda_2_p)-y_2_D)/(-1+teta))-2*nu0_2_tY*fi_2_HM*omega_2_1*exp((y_2_D+log((-1+omega_1_2-omega_2_1*exp(teta*f_1-y_2_D+y_1_D))/teta*(-1+nu0_2_tY)*(-1+teta)/limda_2_p))/(-1+teta)+(-teta*y_2_D-teta*log((-1+omega_1_2-omega_2_1*exp(teta*f_1-y_2_D+y_1_D))/teta*(-1+nu0_2_tY)*(-1+teta)/limda_2_p)-teta*f_1+teta^2*f_1-y_1_D+y_1_D*teta)/(-1+teta))-2*nu0_2_tL*limda_2_p*alfa_2*fi_2_HM)/fi_2_HM/(-1+exp(nu0_2_R-y_2-p_2));
i_2=log(exp(k_2)-exp(k_2-y_2)+exp(k_2-y_2)*nu0_2_SIG);
d_2_B=(fi_2_HB*mu_B_2*exp(-1.*nu0_2_R+r_2_F)+2.*limda_2*exp(nu0_bbb-1.*nu0_2_R-2.*p_2-2.*y_2+3.*r_2_F-1.*k_2)*nu0_2_BF-2.*exp(-3.*y_2-3.*p_2+nu0_bbb+3.*r_2_F-1.*k_2)*limda_2*nu0_2_BF-1.*nu0_2_BF*exp(-1.*p_2-1.*y_2-2.*nu0_2_R+3.*r_2_F)*fi_2_HB*nu0_2_DB+exp(-2.*y_2-2.*p_2+nu0_bbb+2.*r_2_F)*limda_2*nu0_2_DB-1.*nu0_2_BF^2*exp(4.*r_2_F-3.*p_2-3.*y_2-1.*nu0_2_R-1.*k_2)*fi_2_HB-2.*exp(-2.*p_2-2.*y_2-1.*nu0_2_R+3.*r_2_F-1.*k_2)*fi_2_HB*mu_B_2*nu0_2_BF+nu0_2_BF^2*exp(-3.*y_2-3.*p_2+3.*r_2_F-1.*k_2)*fi_2_HB+limda_2*exp(-2.*nu0_2_R+2.*r_2_F)*nu0_2_DB-2.*limda_2*exp(-2.*nu0_2_R+3.*r_2_F-1.*k_2-1.*p_2-1.*y_2)*nu0_2_BF-1.*b_2_B_1*fi_2_HB*exp(-1.*nu0_2_R+r_2_F)+nu0_2_BF*exp(r_2_F-1.*p_2-1.*y_2)*fi_2_HB-1.*b_2_B_2*fi_2_HB*exp(-1.*nu0_2_R+r_2_F)-2.*nu0_2_BF^2*fi_2_HB*exp(-2.*nu0_2_R+3.*r_2_F-1.*k_2-1.*p_2-1.*y_2)-1.*nu0_2_BF*fi_2_HB*exp(-1.*nu0_2_R+r_2_F)-1.*exp(r_2_F-1.*p_2-1.*y_2)*fi_2_HB*mu_B_2+b_2_B_2*exp(r_2_F-1.*p_2-1.*y_2)*fi_2_HB+b_2_B_1*exp(nu0_1_R-1.*y_2-1.*p_1-1.*nu0_2_R+r_2_F)*fi_2_HB+2.*nu0_2_BF^2*exp(-2.*p_2-2.*y_2-2.*nu0_2_R+4.*r_2_F-1.*k_2)*fi_2_HB-1.*limda_2*exp(-1.*nu0_2_R+r_2_F)+exp(r_2_F-1.*p_2-1.*y_2)*limda_2+exp(-1.*p_2-1.*y_2-1.*nu0_2_R+2.*r_2_F)*fi_2_HB*mu_B_2*nu0_2_DB+nu0_2_BF*exp(3.*r_2_F-2.*p_2-2.*y_2-1.*nu0_2_R)*fi_2_HB*nu0_2_DB+2.*b_2_B_1*exp(nu0_1_R-2.*y_2-1.*p_1-2.*nu0_2_R+3.*r_2_F-1.*k_2-1.*p_2)*fi_2_HB*nu0_2_BF-1.*exp(nu0_2_R-2.*y_2-2.*p_2+nu0_bbb+r_2_F)*limda_2-1.*b_2_B_2*exp(-1.*p_2-1.*y_2-1.*nu0_2_R+2.*r_2_F)*fi_2_HB*nu0_2_DB+b_2_B_1*fi_2_HB*exp(-2.*nu0_2_R+2.*r_2_F)*nu0_2_DB-2.*b_2_B_2*fi_2_HB*exp(-2.*nu0_2_R+3.*r_2_F-1.*k_2-1.*p_2-1.*y_2)*nu0_2_BF-1.*nu0_2_BF*exp(-2.*p_2-2.*y_2+2.*r_2_F)*fi_2_HB*nu0_2_DB+b_2_B_2*fi_2_HB*exp(-2.*nu0_2_R+2.*r_2_F)*nu0_2_DB+2.*fi_2_HB*mu_B_2*exp(-2.*nu0_2_R+3.*r_2_F-1.*k_2-1.*p_2-1.*y_2)*nu0_2_BF-1.*b_2_B_1*exp(nu0_1_R-1.*y_2-1.*p_1-2.*nu0_2_R+2.*r_2_F)*fi_2_HB*nu0_2_DB+nu0_2_BF*fi_2_HB*exp(-2.*nu0_2_R+2.*r_2_F)*nu0_2_DB-1.*fi_2_HB*mu_B_2*exp(-2.*nu0_2_R+2.*r_2_F)*nu0_2_DB-2.*b_2_B_1*fi_2_HB*exp(-2.*nu0_2_R+3.*r_2_F-1.*k_2-1.*p_2-1.*y_2)*nu0_2_BF-1.*exp(-1.*p_2-1.*y_2-1.*nu0_2_R+2.*r_2_F)*limda_2*nu0_2_DB+2.*exp(-2.*p_2-2.*y_2-1.*nu0_2_R+3.*r_2_F-1.*k_2)*limda_2*nu0_2_BF-1.*limda_2*exp(nu0_bbb-1.*nu0_2_R-1.*p_2-1.*y_2+2.*r_2_F)*nu0_2_DB+limda_2*exp(nu0_bbb-1.*p_2-1.*y_2+r_2_F)+2.*b_2_B_2*exp(-2.*p_2-2.*y_2-1.*nu0_2_R+3.*r_2_F-1.*k_2)*fi_2_HB*nu0_2_BF)/fi_2_HB/(exp(-1.*nu0_2_R+r_2_F)-1.*exp(-2.*nu0_2_R+2.*r_2_F)*nu0_2_DB+2.*exp(-2.*nu0_2_R+3.*r_2_F-1.*k_2-1.*p_2-1.*y_2)*nu0_2_BF);
d_2_F=(nu0_2_tY*nu0_2_DB/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-2*nu0_2_R+2*r_2_F)-omega_2_1*nu0_2_tY/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(teta*f_1-nu0_2_R+y_1_D-y_2_D+r_2_F)-exp(i_2-nu0_2_R+r_2_F)+1/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-nu0_2_R+r_2_F)+omega_2_1/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(teta*f_1-nu0_2_R+y_1_D-y_2_D+r_2_F)-2*omega_2_1*nu0_2_tY*nu0_2_BF/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(teta*f_1-2*nu0_2_R+y_1_D-k_2+3*r_2_F-y_2-y_2_D-p_2)-nu0_2_DB/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-2*nu0_2_R+2*r_2_F)-nu0_2_tY/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-nu0_2_R+r_2_F)+2*omega_1_2*nu0_2_tY*nu0_2_BF/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)-2*nu0_2_tY*nu0_2_BF/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)+2*omega_2_1*nu0_2_BF/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(teta*f_1-2*nu0_2_R+y_1_D-k_2+3*r_2_F-y_2-y_2_D-p_2)-2*omega_1_2*nu0_2_BF/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)+2*nu0_2_BF/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)-nu0_2_BF*exp(r_2_F-p_2-y_2)+limda_2_p*alfa_2*exp(-2*nu0_2_R+2*r_2_F)*nu0_2_DB-2*limda_2_p*alfa_2*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)*nu0_2_BF+exp(i_2-2*nu0_2_R+2*r_2_F)*nu0_2_DB-nu0_2_BF^2*exp(3*r_2_F-3*p_2-3*y_2-k_2)+nu0_2_BF^2*exp(4*r_2_F-3*p_2-3*y_2-nu0_2_R-k_2)+omega_2_1*nu0_2_tY*nu0_2_DB/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(teta*f_1-2*nu0_2_R+y_1_D-y_2_D+2*r_2_F)-omega_1_2*nu0_2_tY*nu0_2_DB/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-2*nu0_2_R+2*r_2_F)-omega_1_2/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-nu0_2_R+r_2_F)-omega_2_1*nu0_2_DB/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(teta*f_1-2*nu0_2_R+y_1_D-y_2_D+2*r_2_F)+nu0_2_BF*exp(-p_2-y_2-2*nu0_2_R+3*r_2_F)*nu0_2_DB-nu0_2_BF*exp(3*r_2_F-2*p_2-2*y_2-nu0_2_R)*nu0_2_DB+omega_1_2*nu0_2_DB/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-2*nu0_2_R+2*r_2_F)+omega_1_2*nu0_2_tY/(1-nu0_2_tY)/(-1+teta)/(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))*teta*limda_2_p*exp(-nu0_2_R+r_2_F)-2*exp(i_2-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)*nu0_2_BF-2*nu0_2_BF^2*exp(-2*p_2-2*y_2-2*nu0_2_R+4*r_2_F-k_2)-limda_2_p*alfa_2*exp(-nu0_2_R+r_2_F)+nu0_2_BF*exp(2*r_2_F-2*p_2-2*y_2)*nu0_2_DB-nu0_2_BF*exp(-2*nu0_2_R+2*r_2_F)*nu0_2_DB+2*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)*nu0_2_BF^2+nu0_2_BF*exp(-nu0_2_R+r_2_F))/(exp(-nu0_2_R+r_2_F)-exp(-2*nu0_2_R+2*r_2_F)*nu0_2_DB+2*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)*nu0_2_BF);
b_1_B_1=1/2*(-2*fi_1_HM*teta*exp(f_1+p_1+y_2+c_1+nu0_1_MH)+limda_1*teta*exp(f_1-nu0_bbb+2*p_1+2*y_2)+limda_1*nu0_1_tY*exp(f_1-nu0_bbb+2*p_1+2*y_2)+2*fi_1_HM*teta*exp(f_1+nu0_1_G)+2*fi_1_HM*nu0_1_tY*exp(f_1+nu0_1_G)-limda_1*teta*exp(f_1-nu0_bbb+p_1+y_2)-limda_1*nu0_1_tY*exp(f_1-nu0_bbb+p_1+y_2)-limda_1*teta*exp(f_1+p_1+y_2)+exp(f_1)*limda_1*teta-exp(f_1)*limda_1+2*fi_1_HM*exp(f_1+p_1+y_2+c_1+nu0_1_MH)+limda_1*exp(f_1-nu0_bbb+p_1+y_2)-2*fi_1_HM*exp(f_1+nu0_1_G)+limda_1*exp(f_1+p_1+y_2)-limda_1*teta*nu0_1_tY*exp(f_1-nu0_bbb+2*p_1+2*y_2)-2*exp(f_1)*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*nu0_1_tY+2*exp(f_1)*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*teta*nu0_1_tY-2*exp(f_1)*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*teta-limda_1*nu0_1_tY*exp(f_1+p_1+y_2)+2*exp(f_1)*nu0_1_tr*fi_1_HM*nu0_1_tY+2*exp(f_1)*nu0_1_tr*fi_1_HM*teta-exp(f_1)*limda_1*teta*nu0_1_tY+2*b_2_B_1*fi_1_HM*exp(nu0_1_R-y_2-p_1)*teta+2*b_2_B_1*fi_1_HM*exp(nu0_1_R-y_2-p_1)*nu0_1_tY-2*exp(f_1)*nu0_1_tr*fi_1_HM*teta*nu0_1_tY-2*exp(f_1)*nu0_1_tY*teta*limda_1_p*fi_1_HM+2*exp(f_1)*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM-2*exp(f_1)*nu0_1_tr*fi_1_HM-2*b_2_B_1*fi_1_HM*exp(nu0_1_R-y_2-p_1)-2*fi_1_HM*nu0_1_tY*exp(f_1+p_1+y_2+c_1+nu0_1_MH)+2*fi_1_HM*teta*exp(f_1+c_1+nu0_1_MH)+2*fi_1_HM*nu0_1_tY*exp(f_1+c_1+nu0_1_MH)-2*b_2_B_1*fi_1_HM*exp(nu0_1_R-y_2-p_1)*nu0_1_tY*teta+exp(f_1)*limda_1*nu0_1_tY-2*fi_1_HM*exp(f_1+c_1+nu0_1_MH)-2*nu0_1_tY*teta*fi_1_HM*exp(f_1+nu0_1_G)+limda_1*teta*nu0_1_tY*exp(f_1+p_1+y_2)+limda_1*teta*nu0_1_tY*exp(f_1-nu0_bbb+p_1+y_2)+2*nu0_1_tY*teta*fi_1_HM*exp(f_1+p_1+y_2+c_1+nu0_1_MH)-2*nu0_1_tY*teta*fi_1_HM*exp(f_1+c_1+nu0_1_MH)+2*b_2_B_1*fi_1_HM*nu0_1_tY*teta-2*b_2_B_1*fi_1_HM*teta-2*b_2_B_1*fi_1_HM*nu0_1_tY+2*b_2_B_1*fi_1_HM-limda_1*exp(f_1-nu0_bbb+2*p_1+2*y_2))*exp(-f_1)/fi_1_HM/(-1+exp(nu0_1_R-y_2-p_1))/(-1+nu0_1_tY)/(-1+teta);
b_1_G=1/2*(exp(-nu0_bbb+2*p_1+2*y_2)*limda_1+2*exp(p_1+y_2+c_1+nu0_1_MH)*fi_1_HM*teta-2*exp(c_1+nu0_1_MH)*fi_1_HM*teta+exp(-nu0_bbb+p_1+y_2)*limda_1*teta-limda_1*teta-limda_1*nu0_1_tY+2*exp(nu0_1_G)*fi_1_HM*teta*nu0_1_tY+exp(-nu0_bbb+2*p_1+2*y_2)*limda_1*teta*nu0_1_tY-2*exp(p_1+y_2+c_1+nu0_1_MH)*fi_1_HM*teta*nu0_1_tY-limda_1*exp(p_1+y_2)*teta*nu0_1_tY-exp(-nu0_bbb+p_1+y_2)*limda_1*teta*nu0_1_tY+2*exp(c_1+nu0_1_MH)*fi_1_HM*teta*nu0_1_tY-2*exp(c_1+nu0_1_MH)*fi_1_HM*nu0_1_tY-2*exp(nu0_1_G)*fi_1_HM*teta-2*exp(nu0_1_G)*fi_1_HM*nu0_1_tY-exp(-nu0_bbb+2*p_1+2*y_2)*limda_1*teta+limda_1*exp(p_1+y_2)*teta+limda_1*exp(p_1+y_2)*nu0_1_tY+exp(-nu0_bbb+p_1+y_2)*limda_1*nu0_1_tY-exp(-nu0_bbb+2*p_1+2*y_2)*limda_1*nu0_1_tY+limda_1-2*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM-2*nu0_1_tr*fi_1_HM*nu0_1_tY-2*nu0_1_tr*fi_1_HM*teta+2*exp(nu0_1_G)*fi_1_HM+2*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*nu0_1_tY-2*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*teta*nu0_1_tY+limda_1*teta*nu0_1_tY+2*nu0_1_tL*alfa_1*limda_1_p*fi_1_HM*teta+2*nu0_1_tY*teta*limda_1_p*fi_1_HM+2*nu0_1_tr*fi_1_HM*teta*nu0_1_tY+2*exp(c_1+nu0_1_MH)*fi_1_HM-2*exp(p_1+y_2+c_1+nu0_1_MH)*fi_1_HM+2*nu0_1_tr*fi_1_HM-limda_1*exp(p_1+y_2)+2*exp(p_1+y_2+c_1+nu0_1_MH)*fi_1_HM*nu0_1_tY-exp(-nu0_bbb+p_1+y_2)*limda_1)/(-1+exp(nu0_1_R-y_2-p_1))/fi_1_HM/(-1+nu0_1_tY)/(-1+teta);
i_1=log(exp(k_1)-exp(k_1-y_2)+exp(k_1-y_2)*nu0_1_SIG);
d_1_B=(-1.*nu0_1_BF^2*exp(4.*r_1_F-3.*p_1-3.*y_2-1.*nu0_1_R-1.*k_1)*fi_1_HB+2.*nu0_1_BF^2*exp(-2.*nu0_1_R-2.*y_2-2.*p_1+4.*r_1_F-1.*k_1)*fi_1_HB+nu0_1_BF^2*exp(3.*r_1_F-3.*p_1-3.*y_2-1.*k_1)*fi_1_HB-1.*nu0_1_BF*fi_1_HB*exp(-1.*nu0_1_R+r_1_F)+limda_1*exp(-2.*nu0_1_R+2.*r_1_F)*nu0_1_DB+nu0_1_BF*exp(3.*r_1_F-2.*p_1-2.*y_2-1.*nu0_1_R)*fi_1_HB*nu0_1_DB-1.*limda_1*exp(nu0_bbb-1.*nu0_1_R-1.*p_1-1.*y_2+2.*r_1_F)*nu0_1_DB-2.*exp(-1.*nu0_1_R-2.*y_2-2.*p_1+3.*r_1_F-1.*k_1)*fi_1_HB*mu_B_1*nu0_1_BF-1.*limda_1*exp(-1.*nu0_1_R+r_1_F)-1.*exp(-1.*p_1-1.*y_2+r_1_F)*fi_1_HB*mu_B_1-1.*b_1_B_1*fi_1_HB*exp(-1.*nu0_1_R+r_1_F)-1.*b_1_B_2*fi_1_HB*exp(-1.*nu0_1_R+r_1_F)+nu0_1_BF*exp(-1.*p_1-1.*y_2+r_1_F)*fi_1_HB+b_1_B_1*fi_1_HB*exp(-2.*nu0_1_R+2.*r_1_F)*nu0_1_DB-2.*b_1_B_1*fi_1_HB*exp(-2.*nu0_1_R+3.*r_1_F-1.*k_1-1.*p_1-1.*y_2)*nu0_1_BF+b_1_B_2*fi_1_HB*exp(-2.*nu0_1_R+2.*r_1_F)*nu0_1_DB-2.*b_1_B_2*fi_1_HB*exp(-2.*nu0_1_R+3.*r_1_F-1.*k_1-1.*p_1-1.*y_2)*nu0_1_BF+b_1_B_1*exp(-1.*p_1-1.*y_2+r_1_F)*fi_1_HB+2.*fi_1_HB*mu_B_1*exp(-2.*nu0_1_R+3.*r_1_F-1.*k_1-1.*p_1-1.*y_2)*nu0_1_BF+b_1_B_2*exp(nu0_2_R-1.*y_2-1.*p_2-1.*nu0_1_R+r_1_F)*fi_1_HB-1.*b_1_B_2*exp(nu0_2_R-1.*y_2-1.*p_2-2.*nu0_1_R+2.*r_1_F)*fi_1_HB*nu0_1_DB+exp(-2.*y_2-2.*p_1+nu0_bbb+2.*r_1_F)*limda_1*nu0_1_DB+2.*limda_1*exp(nu0_bbb-1.*nu0_1_R-2.*p_1-2.*y_2+3.*r_1_F-1.*k_1)*nu0_1_BF-1.*exp(-1.*nu0_1_R-1.*y_2-1.*p_1+2.*r_1_F)*limda_1*nu0_1_DB-1.*nu0_1_BF*exp(-2.*nu0_1_R-1.*y_2-1.*p_1+3.*r_1_F)*fi_1_HB*nu0_1_DB+nu0_1_BF*fi_1_HB*exp(-2.*nu0_1_R+2.*r_1_F)*nu0_1_DB-1.*fi_1_HB*mu_B_1*exp(-2.*nu0_1_R+2.*r_1_F)*nu0_1_DB-2.*limda_1*exp(-2.*nu0_1_R+3.*r_1_F-1.*k_1-1.*p_1-1.*y_2)*nu0_1_BF-2.*exp(-3.*y_2-3.*p_1+nu0_bbb+3.*r_1_F-1.*k_1)*limda_1*nu0_1_BF-2.*nu0_1_BF^2*fi_1_HB*exp(-2.*nu0_1_R+3.*r_1_F-1.*k_1-1.*p_1-1.*y_2)-1.*b_1_B_1*exp(-1.*nu0_1_R-1.*y_2-1.*p_1+2.*r_1_F)*fi_1_HB*nu0_1_DB-1.*nu0_1_BF*exp(2.*r_1_F-2.*p_1-2.*y_2)*fi_1_HB*nu0_1_DB+limda_1*exp(nu0_bbb-1.*p_1-1.*y_2+r_1_F)+exp(-1.*p_1-1.*y_2+r_1_F)*limda_1-1.*exp(nu0_1_R-2.*y_2-2.*p_1+nu0_bbb+r_1_F)*limda_1+2.*exp(-1.*nu0_1_R-2.*y_2-2.*p_1+3.*r_1_F-1.*k_1)*limda_1*nu0_1_BF+2.*b_1_B_1*exp(-1.*nu0_1_R-2.*y_2-2.*p_1+3.*r_1_F-1.*k_1)*fi_1_HB*nu0_1_BF+exp(-1.*nu0_1_R-1.*y_2-1.*p_1+2.*r_1_F)*fi_1_HB*mu_B_1*nu0_1_DB+2.*b_1_B_2*exp(nu0_2_R-2.*y_2-1.*p_2-2.*nu0_1_R+3.*r_1_F-1.*k_1-1.*p_1)*fi_1_HB*nu0_1_BF+fi_1_HB*mu_B_1*exp(-1.*nu0_1_R+r_1_F))/(exp(-1.*nu0_1_R+r_1_F)-1.*exp(-2.*nu0_1_R+2.*r_1_F)*nu0_1_DB+2.*exp(-2.*nu0_1_R+3.*r_1_F-1.*k_1-1.*p_1-1.*y_2)*nu0_1_BF)/fi_1_HB;
d_1_F=-(exp(i_1-2*nu0_1_R+2*r_1_F+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*omega_2_1+nu0_1_BF*exp(2*r_1_F-2*p_1-2*y_2+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*omega_2_1+nu0_1_BF^2*exp(3*r_1_F-3*p_1-3*y_2-k_1)*omega_1_2+nu0_1_BF*exp(r_1_F-p_1-y_2)*omega_1_2-2*exp(-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*nu0_1_BF^2*omega_1_2+teta*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)-2*alfa_1*limda_1_p*exp(-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*nu0_1_BF*teta*omega_1_2-exp(i_1-nu0_1_R+r_1_F)*teta*omega_1_2+2*exp(i_1-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*nu0_1_BF*omega_1_2+2*nu0_1_BF^2*exp(-2*nu0_1_R-2*y_2-2*p_1+4*r_1_F-k_1)*omega_1_2-nu0_1_BF*exp(-nu0_1_R+r_1_F)*omega_1_2+alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB*omega_2_1-exp(i_1-2*nu0_1_R+2*r_1_F)*nu0_1_DB*omega_1_2-nu0_1_BF^2*exp(4*r_1_F-3*p_1-3*y_2-nu0_1_R-k_1)*omega_1_2+exp(i_1-nu0_1_R+r_1_F)*omega_1_2+nu0_1_BF^2*exp(3*r_1_F-3*p_1-3*y_2-k_1+teta*f_1-y_2_D+y_1_D)+nu0_1_BF*exp(r_1_F-p_1-y_2+teta*f_1-y_2_D+y_1_D)+nu0_1_BF*exp(r_1_F-p_1-y_2+teta*f_1-y_2_D+y_1_D)*teta*omega_2_1-exp(i_1-2*nu0_1_R+2*r_1_F+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*teta*omega_2_1-2*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF^2-alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB+2*nu0_1_BF^2*exp(-2*nu0_1_R-2*y_2-2*p_1+4*r_1_F-k_1+teta*f_1-y_2_D+y_1_D)-nu0_1_BF*exp(r_1_F-p_1-y_2+teta*f_1-y_2_D+y_1_D)*omega_2_1+nu0_1_BF*exp(2*r_1_F-2*p_1-2*y_2)*nu0_1_DB*teta*omega_1_2+2*exp(i_1-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2+teta*f_1-y_2_D+y_1_D)*nu0_1_BF*teta*omega_2_1-exp(i_1-nu0_1_R+r_1_F+teta*f_1-y_2_D+y_1_D)*teta+nu0_1_BF^2*exp(4*r_1_F-3*p_1-3*y_2-nu0_1_R-k_1+teta*f_1-y_2_D+y_1_D)*teta+2*exp(i_1-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2+teta*f_1-y_2_D+y_1_D)*nu0_1_BF-nu0_1_BF^2*exp(4*r_1_F-3*p_1-3*y_2-nu0_1_R-k_1+teta*f_1-y_2_D+y_1_D)*teta*omega_2_1-nu0_1_BF*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)+2*nu0_1_BF^2*exp(-2*nu0_1_R-2*y_2-2*p_1+4*r_1_F-k_1+teta*f_1-y_2_D+y_1_D)*teta*omega_2_1+2*alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF*teta*omega_2_1-nu0_1_BF*exp(-2*nu0_1_R+2*r_1_F)*nu0_1_DB*teta*omega_1_2-nu0_1_BF*exp(2*r_1_F-2*p_1-2*y_2+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*teta*omega_2_1+exp(i_1-2*nu0_1_R+2*r_1_F+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*teta-2*alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF*teta+2*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF^2*teta-nu0_1_BF*exp(3*r_1_F-2*p_1-2*y_2-nu0_1_R+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*teta+nu0_1_BF*exp(3*r_1_F-2*p_1-2*y_2-nu0_1_R+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*teta*omega_2_1+exp(i_1-nu0_1_R+r_1_F+teta*f_1-y_2_D+y_1_D)+alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)*teta*omega_2_1-nu0_1_BF*exp(r_1_F-p_1-y_2+teta*f_1-y_2_D+y_1_D)*teta+alfa_1*limda_1_p*exp(-2*nu0_1_R+2*r_1_F)*nu0_1_DB*teta*omega_1_2+nu0_1_BF*exp(2*r_1_F-2*p_1-2*y_2+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*teta+nu0_1_BF^2*exp(4*r_1_F-3*p_1-3*y_2-nu0_1_R-k_1+teta*f_1-y_2_D+y_1_D)*omega_2_1-nu0_1_BF*exp(-2*nu0_1_R-y_2-p_1+3*r_1_F+teta*f_1-y_2_D+y_1_D)*nu0_1_DB-alfa_1*limda_1_p*exp(-2*nu0_1_R+2*r_1_F)*nu0_1_DB*omega_1_2-exp(i_1-2*nu0_1_R+2*r_1_F+teta*f_1-y_2_D+y_1_D)*nu0_1_DB-2*nu0_1_BF^2*exp(-2*nu0_1_R-2*y_2-2*p_1+4*r_1_F-k_1+teta*f_1-y_2_D+y_1_D)*teta+exp(i_1-nu0_1_R+r_1_F+teta*f_1-y_2_D+y_1_D)*teta*omega_2_1-nu0_1_BF^2*exp(4*r_1_F-3*p_1-3*y_2-nu0_1_R-k_1+teta*f_1-y_2_D+y_1_D)+2*alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF+nu0_1_BF*exp(3*r_1_F-2*p_1-2*y_2-nu0_1_R+teta*f_1-y_2_D+y_1_D)*nu0_1_DB-2*alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF*omega_2_1-nu0_1_BF*exp(-2*nu0_1_R-y_2-p_1+3*r_1_F+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*teta*omega_2_1+2*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF^2*omega_2_1-exp(i_1-nu0_1_R+r_1_F+teta*f_1-y_2_D+y_1_D)*omega_2_1+nu0_1_BF*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB-nu0_1_BF*exp(2*r_1_F-2*p_1-2*y_2+teta*f_1-y_2_D+y_1_D)*nu0_1_DB+nu0_1_BF*exp(-2*nu0_1_R+2*r_1_F)*nu0_1_DB*omega_1_2-nu0_1_BF*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB*teta-alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)*teta-2*exp(i_1-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2+teta*f_1-y_2_D+y_1_D)*nu0_1_BF*omega_2_1+nu0_1_BF*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB*teta*omega_2_1-nu0_1_BF*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB*omega_2_1+nu0_1_BF*exp(-2*nu0_1_R-y_2-p_1+3*r_1_F)*nu0_1_DB*teta*omega_1_2-alfa_1*limda_1_p*exp(-nu0_1_R+r_1_F)*teta*omega_1_2-2*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF^2*teta*omega_2_1-nu0_1_BF*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)*teta*omega_2_1-alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB*teta*omega_2_1-2*nu0_1_BF^2*exp(-2*nu0_1_R-2*y_2-2*p_1+4*r_1_F-k_1+teta*f_1-y_2_D+y_1_D)*omega_2_1+2*teta*limda_1_p*exp(-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*omega_1_2*nu0_1_BF+alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB*teta-teta*limda_1_p*exp(-2*nu0_1_R+2*r_1_F)*omega_1_2*nu0_1_DB+alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)-alfa_1*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)*omega_2_1+teta*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*omega_2_1*nu0_1_DB-2*teta*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*omega_2_1*nu0_1_BF+2*alfa_1*limda_1_p*exp(-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*nu0_1_BF*omega_1_2+nu0_1_BF*exp(-2*nu0_1_R-y_2-p_1+3*r_1_F+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*omega_2_1-nu0_1_BF*exp(3*r_1_F-2*p_1-2*y_2-nu0_1_R+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*omega_2_1+exp(i_1-2*nu0_1_R+2*r_1_F)*nu0_1_DB*teta*omega_1_2-nu0_1_BF*exp(3*r_1_F-2*p_1-2*y_2-nu0_1_R)*nu0_1_DB*teta*omega_1_2-nu0_1_BF*exp(2*r_1_F-2*p_1-2*y_2)*nu0_1_DB*omega_1_2-2*nu0_1_BF^2*exp(-2*nu0_1_R-2*y_2-2*p_1+4*r_1_F-k_1)*teta*omega_1_2-2*exp(i_1-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2+teta*f_1-y_2_D+y_1_D)*nu0_1_BF*teta+nu0_1_BF*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)*omega_2_1-nu0_1_BF*exp(r_1_F-p_1-y_2)*teta*omega_1_2+2*teta*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+3*r_1_F-y_2-k_1-p_1)*nu0_1_BF-2*exp(i_1-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*nu0_1_BF*teta*omega_1_2-nu0_1_BF^2*exp(3*r_1_F-3*p_1-3*y_2-k_1+teta*f_1-y_2_D+y_1_D)*omega_2_1+alfa_1*limda_1_p*exp(-nu0_1_R+r_1_F)*omega_1_2+nu0_1_BF^2*exp(4*r_1_F-3*p_1-3*y_2-nu0_1_R-k_1)*teta*omega_1_2+nu0_1_BF^2*exp(3*r_1_F-3*p_1-3*y_2-k_1+teta*f_1-y_2_D+y_1_D)*teta*omega_2_1-nu0_1_BF^2*exp(3*r_1_F-3*p_1-3*y_2-k_1)*teta*omega_1_2-nu0_1_BF^2*exp(3*r_1_F-3*p_1-3*y_2-k_1+teta*f_1-y_2_D+y_1_D)*teta+nu0_1_BF*exp(-nu0_1_R+r_1_F)*teta*omega_1_2+nu0_1_BF*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)*teta+teta*limda_1_p*exp(-nu0_1_R+r_1_F)*omega_1_2+nu0_1_BF*exp(-2*nu0_1_R-y_2-p_1+3*r_1_F+teta*f_1-y_2_D+y_1_D)*nu0_1_DB*teta+nu0_1_BF*exp(3*r_1_F-2*p_1-2*y_2-nu0_1_R)*nu0_1_DB*omega_1_2-teta*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-nu0_1_R+r_1_F)*omega_2_1-teta*limda_1_p*exp(teta*f_1-y_2_D+y_1_D-2*nu0_1_R+2*r_1_F)*nu0_1_DB-nu0_1_BF*exp(-2*nu0_1_R-y_2-p_1+3*r_1_F)*nu0_1_DB*omega_1_2+2*exp(-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*nu0_1_BF^2*teta*omega_1_2)/(exp(-nu0_1_R+r_1_F)-exp(-2*nu0_1_R+2*r_1_F)*nu0_1_DB+2*exp(-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*nu0_1_BF)/(-exp(teta*f_1-y_2_D+y_1_D)+omega_2_1*exp(teta*f_1-y_2_D+y_1_D)-omega_1_2)/(-1+teta);
m_2_H=log(1/2*(2*fi_2_HM*exp(nu0_bbb+c_2+nu0_2_MH)+exp(nu0_bbb)*limda_2-limda_2*exp(p_2+y_2))/fi_2_HM)+p_2+y_2-nu0_bbb;
l_2=(-nu0_2_L+log(limda_2_p*alfa_2*limda_2*(1-nu0_2_tL)))/(1+omega_2_L);
w_2=log(1/limda_2/(1-nu0_2_tL))+l_2*omega_2_L+nu0_2_L;
b_2_H=(limda_2*exp(nu0_bbb+nu0_2_R-1.*p_2-1.*y_2)+fi_2_HB*mu_B_2-1.*limda_2)/fi_2_HB;
m_1_H=log(1/2*(2*fi_1_HM*exp(nu0_bbb+c_1+nu0_1_MH)+exp(nu0_bbb)*limda_1-limda_1*exp(p_1+y_2))/fi_1_HM)-nu0_bbb+p_1+y_2;
l_1=(-nu0_1_L+log(alfa_1*limda_1_p*limda_1*(1-nu0_1_tL)))/(1+omega_1_L);
w_1=log(1/limda_1/(1-nu0_1_tL))+l_1*omega_1_L+nu0_1_L;
b_1_H=(limda_1*exp(nu0_bbb+nu0_1_R-1.*p_1-1.*y_2)+fi_1_HB*mu_B_1-1.*limda_1)/fi_1_HB;
p_2_P=(y_2_D+log((1-nu0_2_tY)*(-1+teta)*(1-omega_1_2+omega_2_1*exp(teta*f_1-y_2_D+y_1_D))/teta/limda_2_p))/(-1+teta);
limda_2_d=(exp(p_2_P-teta*p_2_P+y_2_D)-omega_1_2*exp(p_2_P-teta*p_2_P+y_2_D)+omega_2_1*exp(p_2_P-teta*p_2_P+teta*f_1+y_1_D)-nu0_2_tY*exp(p_2_P-teta*p_2_P+y_2_D)+nu0_2_tY*omega_1_2*exp(p_2_P-teta*p_2_P+y_2_D)-nu0_2_tY*omega_2_1*exp(p_2_P-teta*p_2_P+teta*f_1+y_1_D)-limda_2_p)/(exp(-teta*p_2_P+y_2_D)-exp(-teta*p_2_P+y_2_D)*omega_1_2+omega_2_1*exp(-teta*p_2_P+teta*f_1+y_1_D));
y_2_P=log(exp(-teta*p_2_P+y_2_D)-exp(-teta*p_2_P+y_2_D)*omega_1_2+omega_2_1*exp(-teta*p_2_P+teta*f_1+y_1_D));
p_1_P=(log((-1+nu0_1_tY)*(-1+teta)*(-exp(teta*f_1-y_2_D+y_1_D)+omega_2_1*exp(teta*f_1-y_2_D+y_1_D)-omega_1_2)/teta/limda_1_p)-teta*f_1+y_2_D)/(-1+teta);
limda_1_d=(-exp(p_1_P-teta*p_1_P+y_1_D)+omega_2_1*exp(p_1_P-teta*p_1_P+y_1_D)-omega_1_2*exp(p_1_P-teta*p_1_P-teta*f_1+y_2_D)+nu0_1_tY*exp(p_1_P-teta*p_1_P+y_1_D)-nu0_1_tY*omega_2_1*exp(p_1_P-teta*p_1_P+y_1_D)+nu0_1_tY*omega_1_2*exp(p_1_P-teta*p_1_P-teta*f_1+y_2_D)+limda_1_p)/(-exp(-teta*p_1_P+y_1_D)+exp(-teta*p_1_P+y_1_D)*omega_2_1-omega_1_2*exp(-teta*p_1_P-teta*f_1+y_2_D));
y_1_P=log(exp(-teta*p_1_P+y_1_D)-exp(-teta*p_1_P+y_1_D)*omega_2_1+omega_1_2*exp(-teta*p_1_P-teta*f_1+y_2_D));
b_2_F=nu0_2_BF;
d_2_FB=(2*fi_2_FB*b_2_F-2*fi_2_FB*nu0_2_BF+exp(-nu0_2_R+r_2_F)-1)/(exp(-nu0_2_R+r_2_F)-exp(-2*nu0_2_R+2*r_2_F)*nu0_2_DB+2*exp(-2*nu0_2_R+3*r_2_F-k_2-p_2-y_2)*b_2_F);
k_2_B=log(b_2_F/nu0_2_BK);
b_1_F=nu0_1_BF;
d_1_FB=(2*fi_1_FB*b_1_F-2*fi_1_FB*nu0_1_BF+exp(-nu0_1_R+r_1_F)-1)/(exp(-nu0_1_R+r_1_F)-exp(-2*nu0_1_R+2*r_1_F)*nu0_1_DB+2*exp(-2*nu0_1_R+3*r_1_F-k_1-p_1-y_2)*b_1_F);
k_1_B=log(b_1_F/nu0_1_BK);
end

p_2=p_2_SS;
y_2_D=y_2_D_SS;
limda_2_p=limda_2_p_SS;
limda_2=limda_2_SS;
y_1_D=y_1_D_SS;
limda_1_p=limda_1_p_SS;
limda_1=limda_1_SS;
b_1_B_2=b_1_B_2_SS;
b_2_B_2=b_2_B_2_SS;
b_2_G=b_2_G_SS;
i_2=i_2_SS;
d_2_B=d_2_B_SS;
d_2_F=d_2_F_SS;
b_1_B_1=b_1_B_1_SS;
b_1_G=b_1_G_SS;
i_1=i_1_SS;
d_1_B=d_1_B_SS;
d_1_F=d_1_F_SS;
m_2_H=m_2_H_SS;
l_2=l_2_SS;
w_2=w_2_SS;
b_2_H=b_2_H_SS;
m_1_H=m_1_H_SS;
l_1=l_1_SS;
w_1=w_1_SS;
b_1_H=b_1_H_SS;
p_2_P=p_2_P_SS;
limda_2_d=limda_2_d_SS;
y_2_P=y_2_P_SS;
p_1_P=p_1_P_SS;
limda_1_d=limda_1_d_SS;
y_1_P=y_1_P_SS;
b_2_F=b_2_F_SS;
d_2_FB=d_2_FB_SS;
k_2_B=k_2_B_SS;
b_1_F=b_1_F_SS;
d_1_FB=d_1_FB_SS;
k_1_B=k_1_B_SS;


r_2_H=nu0_2_R;
r_1_H=nu0_1_R;
limda_2_k=1;
limda_1_k=1;
limda_2_BK=0;
limda_1_BK=0;
limda_2_b=1;
limda_1_b=1;
limda_2_BB=1;
limda_1_BB=1;
tr_2=nu0_2_tr;
g_2=nu0_2_G;
tax_2_Y=nu0_2_tY;
tax_2_L=nu0_2_tL;
r_2=nu0_2_R;
tr_1=nu0_1_tr;
g_1=nu0_1_G;
tax_1_Y=nu0_1_tY;
tax_1_L=nu0_1_tL;
r_1=nu0_1_R;
%y_2=y_1;
%y_1=nu0_Y;
f_2=0;


zzobs_1_C=c_1-y_1_P-p_1_P;
zzobs_1_Ex=log(exp(y_1_P)-omega_1_1*exp(-teta_1*(p_1_P)+y_1_D))-y_1_P;
zzobs_1_fn=-(-p_1+p_2);
zzobs_1_fr=-(0);
zzobs_1_G=g_1-y_1_P-p_1_P;
zzobs_1_I=i_1-y_1_P-p_1_P;
zzobs_1_Im=log(exp(y_1_D)-omega_1_1*exp(-teta_1*(p_1_P)+y_1_D+p_1_P))-y_1_P-p_1_P;   
zzobs_1_M=z_Y+p_1;
zzobs_1_P=p_1;
zzobs_1_R=r_1;
zzobs_1_Y=y_1;

zzobs_1_dR=r_1_F-r_1;
zzobs_1_BK=p_1+z_Y;
zzobs_1_tL=tax_1_L;
zzobs_1_tY=tax_1_Y;
zzobs_1_tr=tr_1/exp(p_2_P+y_2_P);
zzobs_1_WL=w_1+l_1-y_1_P-p_1_P;


qF_z_Y=z_Y;
qF_p_1=p_1;
qF_z_1_DB=z_1_DB;
qF_limda_1_b=limda_1_b;
qF_limda_1_BB=limda_1_BB;
qF_limda_1_BK=limda_1_BK;

qF_p_2=p_2;
qF_z_2_DB=z_2_DB;
qF_limda_2_b=limda_2_b;
qF_limda_2_BB=limda_2_BB;
qF_limda_2_BK=limda_2_BK;

qP_b_1_F=b_1_F;
qP_r_1_F=r_1_F;
qP_b_2_F=b_2_F;
qP_r_2_F=r_2_F;


ys=[qF_z_Y,qF_p_1,qF_z_1_DB,qF_limda_1_b,qF_limda_1_BB,qF_limda_1_BK,qF_p_2,qF_z_2_DB,qF_limda_2_b,qF_limda_2_BB,qF_limda_2_BK,qP_b_1_F,qP_r_1_F,qP_b_2_F,qP_r_2_F,b_1_B_1,b_1_B_2,b_1_F,b_1_G,b_1_H,b_2_B_1,b_2_B_2,b_2_F,b_2_G,b_2_H,c_1,c_2,d_1_B,d_1_F,d_1_FB,d_2_B,d_2_F,d_2_FB,f_1,f_2,g_1,g_2,i_1,i_2,k_1,k_1_B,k_2,k_2_B,l_1,l_2,limda_1,limda_1_b,limda_1_BB,limda_1_BK,limda_1_d,limda_1_k,limda_1_p,limda_2,limda_2_b,limda_2_BB,limda_2_BK,limda_2_d,limda_2_k,limda_2_p,m_1_H,m_2_H,p_1,p_1_P,p_2,p_2_P,r_1,r_1_F,r_1_H,r_2,r_2_F,r_2_H,tax_1_L,tax_1_Y,tax_2_L,tax_2_Y,tr_1,tr_2,w_1,w_2,y_1,y_1_D,y_1_P,y_2,y_2_D,y_2_P,z_1_alfa,z_1_bbb,z_1_BF,z_1_BK,z_1_DB,z_1_DBM,z_1_G,z_1_L,z_1_MH,z_1_R,z_1_SIG,z_1_tL,z_1_tr,z_1_tY,z_1_Y,z_2_alfa,z_2_bbb,z_2_BF,z_2_BK,z_2_DB,z_2_DBM,z_2_G,z_2_L,z_2_MH,z_2_R,z_2_SIG,z_2_tL,z_2_tr,z_2_tY,z_2_Y,z_Y,zzobs_1_C,zzobs_1_Ex,zzobs_1_fn,zzobs_1_fr,zzobs_1_G,zzobs_1_I,zzobs_1_Im,zzobs_1_M,zzobs_1_P,zzobs_1_R,zzobs_1_Y,zzobs_1_dR,zzobs_1_BK,zzobs_1_tL,zzobs_1_tY,zzobs_1_tr,zzobs_1_WL].';
%ZZ=ys([32,45,93,70,38,47,83,33,39,3,50,8,53],1);
%ys=[ys;ZZ];
%ER=feval([M_.fname '_static'],ys,[zeros(27,1)], M_.params);

%keyboard;
if 1==2
    ER=0;
    ER=ER+abs(d_1_FB)*(d_1_FB<0);
    ER=ER+abs(d_2_FB)*(d_2_FB<0);
    ER=ER+abs(b_1_H)*(b_1_H<0);
    ER=ER+abs(b_2_H)*(b_2_H<0);
    ER=ER+abs(b_1_F)*(b_1_F<0);
    ER=ER+abs(b_2_F)*(b_2_F<0);
    ER=ER+abs(b_1_G)*(b_1_G>0);
    ER=ER+abs(b_2_G)*(b_2_G>0);
    
    ys=ys+ER*(-1)^0.5;

end

return



