%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%**计算9个核心变量稳态

function ff=mysecond_steadystate_core(xxx,omegac_p,etac_p,lambdamc_p,omegai_p,etai_p,...
                lambdami_p,omegax_p,etax_p,lambdamx_p,lambdax_p,lambdad_p,muz_p,alpha_p,...
            delta_p,etaf_p,al_p,sigmal_p,lambdaw_p,b_p,beta_p,pitarget_p)
%%
%backing out the guesses for the variables to be solved for
k=xxx(1);
h=xxx(2);
y=xxx(3);
c=xxx(4);
ystar=xxx(5);



%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%辅助简化,或者用到的其他

sbar=1;
gammac=(1-omegac_p+omegac_p*((1-omegac_p)/((sbar*lambdamc_p)^(etac_p-1)-omegac_p)))^(1/(1-etac_p));
gammai=(1-omegai_p+omegai_p*(lambdami_p*sbar*gammac)^(1-etai_p))^(1/(1-etai_p));
gammax=lambdax_p*(1-omegax_p+omegax_p*(lambdamx_p*sbar*gammac)^(1-etai_p))^(1/(1-etai_p))/(sbar*gammac);
rd=pitarget_p*muz_p/beta_p;

rk=(rd/pitarget_p-1+delta_p)*gammai  ;
w=rk*(1-alpha_p)/(alpha_p*muz_p)*1/((lambdad_p*rk*muz_p^(alpha_p-1)/alpha_p)^(1/(1-alpha_p)));

%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%对应的10个方程  
equ1=(lambdad_p*rk*muz_p^(alpha_p-1)/alpha_p)^(1/(1-alpha_p))-h/k;

%equ4,总产出表达
equ2=(k/muz_p)^alpha_p*h^(1-alpha_p)/lambdad_p-y;


%equ5,国内总资源约束,空格为政府支出占比???
equ3=0.2*y+(1-omegac_p)*gammac^etac_p*c+(1-omegai_p)...
*gammai^etai_p*k*(1-(1-delta_p)/muz_p)+(1-omegax_p)...
*(omegax_p*(lambdamx_p*sbar*gammac)^(1-etax_p)+1-omegax_p)^(etax_p/(1-etax_p))*gammax^(-etaf_p)*ystar-y;

%equ6，工资
equ4=al_p*h^sigmal_p*lambdaw_p*gammac/(1/((1-b_p/muz_p)*c)-beta_p*b_p/((muz_p-b_p)*c))-w;


%equ8，nfa稳态,得到稳态进口占产出的稳态,最后空格为稳态的出口占比???
equ5=gammax*gammax^(-etaf_p)*ystar-( omegai_p*(lambdamc_p*sbar)^(-etai_p)*c + omegai_p*(lambdami_p*sbar)^(-etai_p)* (1-(1-delta_p)/muz_p)*k...
+omegax_p*((omegax_p*(lambdamx_p*sbar*gammac)^(1-etax_p)+1-omegax_p)^(1/(1-etax_p))/(lambdamx_p*sbar*gammac))^etax_p*gammax^(-etaf_p)*ystar);





%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%constructing the matrix to be solved with e.g. fsolve
ff=[equ1
    equ2
    equ3
    equ4
    equ5];

