It’s about the composite parameters set before the model
-command :
RbSS=(BETAS*(1-gama*(1-RRRss))+BETAF*gama)/(BETAS*BETAF);
z1=J/(1-BETAS); %qHstoCs
z2=J*RbSS/((1-BETAB)*RbSS-Kss*(1-BETAB*RbSS)); %qHbtoCb
z3=(J*Kss)/(RbSS*(1-BETAB)-Kss*(1-BETAB*RbSS)); %bSStoCbSS
z4=1-(1-RbSS)*(z3); %NbwbtoCb
z5=BETAS/(BETAS*gama*(1-RRRss)+BETAS*RbSS-BETAS-gama); %bSStoCfSS
z6=((1-ALPHA)/(XSS*z4)); %CbtoY
z7=(z3/z5)*z6; %CftoY
z8=1-z6-z7; %CstoY
z9=(z1*z8)/(z2*z6); %HsSS/HbSS
z90=z3*z6; %btoY
z91=((1/BETAS)-1)*gama*z90+((ALPHA-1)/XSS)+1; %CstoY
z92=z1*z91; %qHstoY
z93=ALPHA/XSS; %wsNstoY
See