function [ys,check] = ALEstim_slim_Multiplicative_steadystate(ys,junk)

% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This function computes the steady-state values for relevant variables in the model. 
%
% 
%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
global M_ check 


%%
%Read parameter values from M_
Al 	=	M_.params(	1	);
Aq 	=	M_.params(	2	);
b 	=	M_.params(	3	);
beta 	=	M_.params(	4	);
cbar 	=	M_.params(	5	);
cdbar 	=	M_.params(	6	);
cdfbar 	=	M_.params(	7	);
csbar 	=	M_.params(	8	);
D1 	=	M_.params(	9	);
D2 	=	M_.params(	10	);
D3 	=	M_.params(	11	);
D4 	=	M_.params(	12	);
delta 	=	M_.params(	13	);
epspolicy 	=	M_.params(	14	);
epsz 	=	M_.params(	15	);
etac 	=	M_.params(	16	);
etaf 	=	M_.params(	17	);
etai 	=	M_.params(	18	);
etamc 	=	M_.params(	19	);
etami 	=	M_.params(	20	);
gammacd 	=	M_.params(	21	);
gammacmc 	=	M_.params(	22	);
gammadc 	=	M_.params(	23	);
gammadi 	=	M_.params(	24	);
gammaf 	=	M_.params(	25	);
gammaid 	=	M_.params(	26	);
gammaimi 	=	M_.params(	27	);
gammamcd 	=	M_.params(	28	);
gammamid 	=	M_.params(	29	);
gammapi1 	=	M_.params(	30	);
gammapi2 	=	M_.params(	31	);
gammapi3 	=	M_.params(	32	);
gammaq1 	=	M_.params(	33	);
gammaq2 	=	M_.params(	34	);
gammaq3 	=	M_.params(	35	);
gammax 	=	M_.params(	36	);
gammay1 	=	M_.params(	37	);
gammay2 	=	M_.params(	38	);
gammay3 	=	M_.params(	39	);
gbar 	=	M_.params(	40	);
Hbar 	=	M_.params(	41	);
ibar 	=	M_.params(	42	);
idbar 	=	M_.params(	43	);
idfbar 	=	M_.params(	44	);
isbar 	=	M_.params(	45	);
kbar 	=	M_.params(	46	);
koverH 	=	M_.params(	47	);
lambdaf 	=	M_.params(	48	);
lambdaw 	=	M_.params(	49	);
mbar 	=	M_.params(	50	);
mcbar 	=	M_.params(	51	);
mu 	=	M_.params(	52	);
muz 	=	M_.params(	53	);
nubar 	=	M_.params(	54	);
omegac 	=	M_.params(	55	);
omegai 	=	M_.params(	56	);
phi 	=	M_.params(	57	);
phitilde 	=	M_.params(	58	);
phiz 	=	M_.params(	59	);
piebar 	=	M_.params(	60	);
piecbar 	=	M_.params(	61	);
Pii1 	=	M_.params(	62	);
Pii2 	=	M_.params(	63	);
Pii3 	=	M_.params(	64	);
Piim 	=	M_.params(	65	);
Piimc 	=	M_.params(	66	);
PiimcPc 	=	M_.params(	67	);
Piimi 	=	M_.params(	68	);
PiimiPi 	=	M_.params(	69	);
PiiTot 	=	M_.params(	70	);
Pkprime 	=	M_.params(	71	);
qbar 	=	M_.params(	72	);
ratioc 	=	M_.params(	73	);
ratiog 	=	M_.params(	74	);
ratioi 	=	M_.params(	75	);
ratiox 	=	M_.params(	76	);
ratioz 	=	M_.params(	77	);
Rbar 	=	M_.params(	78	);
Rfbar 	=	M_.params(	79	);
rhoeps 	=	M_.params(	80	);
rhoetamc 	=	M_.params(	81	);
rhoetami 	=	M_.params(	82	);
rhoiota 	=	M_.params(	83	);
rholambda 	=	M_.params(	84	);
rhomuz 	=	M_.params(	85	);
rhonu 	=	M_.params(	86	);
rhophitilde 	=	M_.params(	87	);
rhoR 	=	M_.params(	88	);
rhozetac 	=	M_.params(	89	);
rhozetah 	=	M_.params(	90	);
rhozetaq 	=	M_.params(	91	);
rkbar 	=	M_.params(	92	);
sharek 	=	M_.params(	93	);
sigmaa 	=	M_.params(	94	);
sigmal 	=	M_.params(	95	);
sigmaq 	=	M_.params(	96	);
SS 	=	M_.params(	97	);
tauc 	=	M_.params(	98	);
tauk 	=	M_.params(	99	);
tauw 	=	M_.params(	100	);
tauy 	=	M_.params(	101	);
wbar 	=	M_.params(	102	);
xbar 	=	M_.params(	103	);
xid 	=	M_.params(	104	);
ximc 	=	M_.params(	105	);
ximi 	=	M_.params(	106	);
xiw 	=	M_.params(	107	);
ybar 	=	M_.params(	108	);
ybar1 	=	M_.params(	109	);
ybar2 	=	M_.params(	110	);
ybar3 	=	M_.params(	111	);
zbar 	=	M_.params(	112	);
bw 	=	M_.params(	113	);
eta0 	=	M_.params(	114	);
eta1 	=	M_.params(	115	);
eta2 	=	M_.params(	116	);
eta3 	=	M_.params(	117	);
eta4        =	M_.params(	118	);
eta5        =	M_.params(	119	);
eta6        =	M_.params(	120	);
eta7        =	M_.params(	121	);
eta8        =	M_.params(	122	);
eta9        =	M_.params(	123	);
eta10       =	M_.params(	124	);
eta11       =	M_.params(	125	);
eta12       =	M_.params(	126	);
eta13       =	M_.params(	127	);
eta14       =	M_.params(	128	);
eta3bar 	=	M_.params(	129	);
rhotauk 	=	M_.params(	130	);
rhotauy 	=	M_.params(	131	);
rhotauc 	=	M_.params(	132	);
rhotauw 	=	M_.params(	133	);
rhog        =	M_.params(	134	);
rhoepsR 	=	M_.params(	135	);
rhopistar 	=	M_.params(	136	);
rhoystar 	=	M_.params(	137	);
rhoRstar 	=	M_.params(	138	);
cmbar       =	M_.params(	139	);
imbar       =	M_.params(	140	);
rholambdax 	=	M_.params(	141	);
rhozetatilde 	=	M_.params(	142	);
xix         =	M_.params(	143	);
kappaw      =	M_.params(	144	);
kappad      =	M_.params(	145	);
rhopietarget 	=	M_.params(	146	);
kappax      =	M_.params(	147	);
rdeltapi 	=	M_.params(	148	);
rpie        =	M_.params(	149	);
rx          =	M_.params(	150	);
ry          =	M_.params(	151	);
rdeltay 	=	M_.params(	152	);
kappamc 	=	M_.params(	153	);
kappami 	=	M_.params(	154	);
xie         =	M_.params(	155	);
phitildes 	=	M_.params(	156	);
ForLag111 	=	M_.params(	157	);
ForLag112 	=	M_.params(	158	);
ForLag113 	=	M_.params(	159	);
ForLag121 	=	M_.params(	160	);
ForLag122 	=	M_.params(	161	);
ForLag123 	=	M_.params(	162	);
ForLag131 	=	M_.params(	163	);
ForLag132 	=	M_.params(	164	);
ForLag133 	=	M_.params(	165	);
ForLag211 	=	M_.params(	166	);
ForLag212 	=	M_.params(	167	);
ForLag213 	=	M_.params(	168	);
ForLag221 	=	M_.params(	169	);
ForLag222 	=	M_.params(	170	);
ForLag223 	=	M_.params(	171	);
ForLag231 	=	M_.params(	172	);
ForLag232 	=	M_.params(	173	);
ForLag233 	=	M_.params(	174	);
ForLag311 	=	M_.params(	175	);
ForLag312 	=	M_.params(	176	);
ForLag313 	=	M_.params(	177	);
ForLag321 	=	M_.params(	178	);
ForLag322 	=	M_.params(	179	);
ForLag323 	=	M_.params(	180	);
ForLag331 	=	M_.params(	181	);
ForLag332 	=	M_.params(	182	);
ForLag333 	=	M_.params(	183	);
ForLag411 	=	M_.params(	184	);
ForLag412 	=	M_.params(	185	);
ForLag413 	=	M_.params(	186	);
ForLag421 	=	M_.params(	187	);
ForLag422 	=	M_.params(	188	);
ForLag423 	=	M_.params(	189	);
ForLag431 	=	M_.params(	190	);
ForLag432 	=	M_.params(	191	);
ForLag433 	=	M_.params(	192	);
ForShock11 	=	M_.params(	193	);
ForShock12 	=	M_.params(	194	);
ForShock13 	=	M_.params(	195	);
ForShock21 	=	M_.params(	196	);
ForShock22 	=	M_.params(	197	);
ForShock23 	=	M_.params(	198	);
ForShock31 	=	M_.params(	199	);
ForShock32 	=	M_.params(	200	);
ForShock33 	=	M_.params(	201	);
FisLag111 	=	M_.params(	202	);
FisLag112 	=	M_.params(	203	);
FisLag113 	=	M_.params(	204	);
FisLag121 	=	M_.params(	205	);
FisLag122 	=	M_.params(	206	);
FisLag123 	=	M_.params(	207	);
FisLag131 	=	M_.params(	208	);
FisLag132 	=	M_.params(	209	);
FisLag133 	=	M_.params(	210	);
FisLag211 	=	M_.params(	211	);
FisLag212 	=	M_.params(	212	);
FisLag213 	=	M_.params(	213	);
FisLag221 	=	M_.params(	214	);
FisLag222 	=	M_.params(	215	);
FisLag223 	=	M_.params(	216	);
FisLag231 	=	M_.params(	217	);
FisLag232 	=	M_.params(	218	);
FisLag233 	=	M_.params(	219	);
FisShock11 	=	M_.params(	220	);
FisShock12 	=	M_.params(	221	);
FisShock13 	=	M_.params(	222	);
FisShock14 	=	M_.params(	223	);
FisShock15 	=	M_.params(	224	);
FisShock21 	=	M_.params(	225	);
FisShock22 	=	M_.params(	226	);
FisShock23 	=	M_.params(	227	);
FisShock24 	=	M_.params(	228	);
FisShock25 	=	M_.params(	229	);
FisShock31 	=	M_.params(	230	);
FisShock32 	=	M_.params(	231	);
FisShock33 	=	M_.params(	232	);
FisShock34 	=	M_.params(	233	);
FisShock35 	=	M_.params(	234	);
FisShock41 	=	M_.params(	235	);
FisShock42 	=	M_.params(	236	);
FisShock43 	=	M_.params(	237	);
FisShock44 	=	M_.params(	238	);
FisShock45 	=	M_.params(	239	);
FisShock51 	=	M_.params(	240	);
FisShock52 	=	M_.params(	241	);
FisShock53 	=	M_.params(	242	);
FisShock54 	=	M_.params(	243	);
FisShock55 	=	M_.params(	244	);
epstauk 	=	M_.params(	245	);
epstauw	    =	M_.params(	246	);
lambdamc	=	M_.params(	247	);
lambdami	=	M_.params(	248	);


%%


etamc = lambdamc/(lambdamc-1);%transforming markups to demand elasticities used in the programs
etami = lambdami/(lambdami-1);%transforming markups to demand elasticities used in the programs

kappamc=kappad;
kappami=kappad;
kappax=kappad;

%foreign terms of trade  = Nominal exhange rate times foreign price divided by domestic price
gammaf  = 1;

%equation (3.3)
piebar = mu/muz;
piecbar = piebar;  %CPI inflation

%equation (3.2)
Rbar = (piebar*muz-tauk*beta)/((1-tauk)*beta);  %New version 2003-11-28

%equation (3.15)
Rfbar = nubar*Rbar +(1-nubar);

%equation (3.7)
gammacd = ( (1-omegac) + omegac*(etamc/(etamc-1))^(1-etac) )^(1/(1-etac));

%equation (3.8)
gammacmc = ( (1-omegac)*((etamc-1)/etamc)^(1-etac) + omegac )^(1/(1-etac));%New version 2003-11-19

%equation (3.9)
gammaid = ( (1-omegai) + omegai*(etami/(etami-1))^(1-etai) )^(1/(1-etai));

%equation (3.10)
gammaimi = ( (1-omegai)*((etami-1)/etami)^(1-etai) + omegai )^(1/(1-etai)); %New version 2003-11-19

%equation (3.11)
gammamcd = ( ((1-omegac)+omegac*(etamc/(etamc-1))^(1-etac))/((1-omegac)*((etamc-1)/etamc)^(1-etac)+omegac) )^(1/(1-etac));
gammamcd = etamc/(etamc - 1);

%equation (3.12)
gammamid = ( ((1-omegai)+omegai*(etami/(etami-1))^(1-etai))/((1-omegai)*((etami-1)/etami)^(1-etai)+omegai) )^(1/(1-etai));
gammamid = etami/(etami - 1);

%checking that the calculations of the relative prices are consistent
if abs(gammacd/gammacmc-gammamcd) > 1e-9;
   error('The steady state calculations of the consumption relative prices are not correct...');
elseif abs(gammaid/gammaimi-gammamid) > 1e-9;
   error('The steady state calculations of the investment relative prices are not correct...');
end;
   
%equation (3.6)
Pkprime = gammaid;

%equation (3.14)
rkbar = (muz*Pkprime-beta*(1-delta)*Pkprime)/((1-tauk)*beta);

%solution for wbarRfbar from equation 16
wbarRfbar = (lambdaf*( (1/(1-sharek))^(1-sharek) )*( (1/sharek)^sharek )*rkbar^sharek)^(1/(sharek-1));

%solution for wbar
wbar = wbarRfbar/Rfbar;

%solution for k/H from equation 13
koverH = (sharek/(1-sharek))*muz*wbar*Rfbar/rkbar;


%Profit for imported firms
PiimcPc = omegac*(gammamcd-1)*gammacmc^etac;
PiimiPi = omegai*(gammamid-1)*gammaimi^etai;
if PiimcPc < 0 | PiimiPi < 0;
   error('Profits for the importing firm are negative, should not be possible...');
end;

%definition of composite coefficients, equation (3.25)
D1 = (1-omegac)*gammacd^etac+omegac*gammacmc^etac;
D2 = ((1-ratiog)/lambdaf)*(muz^-sharek)*(koverH^sharek) - ...
     ((1-omegai)*(gammaid^etai)+omegai*(gammaimi^etai))*(1-(1-delta)/muz)*koverH;
D3 = ( ((1-tauy)*(1/lambdaw)*(wbar/(1+tauw)))/Al )^(1/sigmal);
D4 = ((muz-beta*b)/((1+tauc)*(muz-b)))*(gammacd^-1);

%equation (3.26)
Hbar = (D3*(D4^(1/sigmal)*((D2/D1)^(-1/sigmal))))^(sigmal/(1+sigmal));
cbar = (D2/D1)*Hbar;
phiz = (1/cbar)*D4;

%equation (3.27)
ybar = (1/lambdaf)*(muz^-sharek)*(koverH^sharek)*Hbar;



%equation (3.18)
phi = (lambdaf-1)*ybar;
if abs(phi-((lambdaf-1)/lambdaf)*(muz^-sharek)*(koverH^sharek)*Hbar) > 1e-6;
   error('Something is wrong with the steady state calcuations of phi');
end;

%solution for kbar
kbar = koverH*Hbar;

%equation (3.20)
ibar = (1-(1-delta)/muz)*kbar;

%equation (3.28)
gbar = ratiog*ybar;

%Import: cons & invest
cmbar = omegac*(gammacmc^etac)*cbar;
imbar = omegai*(gammaimi^etai)*ibar;

%equation (3.24) (xtilde = cxbar + ixbar)
xtilde = cmbar+imbar;
zbar = xtilde; %export equal import in steady state
xbar = xtilde;
% NOTE that xtilde = ystarbar since 
% 1. (cmbar + imbar = cxbar + ixbar)
% 2. zstar = z
% 3. Px = Pstar
% 
% cmbar + imbar = omegac*(gammacmc^etac)*cbar + omegai*(gammaimi^etai)*ibar
% cxbar + ixbar = (gammaxstar^-etaf)*ystar*(zstar/z)
%


%domestic: cons & invest
cdbar = (1-omegac)*(gammacd^etac)*cbar;
idbar = (1-omegai)*(gammaid^etai)*ibar;

%NOTE #1: ChangeHere 2004-03-23 by JL: We have that mbar = nubar*wbar*Hbar+qbar
%since piebar*muz/mu = 1. The calibration of qbar and mbar from eqs 27 and
%28 are problematic since when Aq is fixed, then when sigmaq change in the
%estimation, then the qbar/mbar ratio in the model can be very different from
%0.4, which is roughly the relevant value in the data. Therefore, we change
%the code here to calibrate qbar/mbar rather than letting the ratio of qbar/mbar
%change considerably in the estimation.
%NOTE #2: We change the notation here so that Aq now refer to the desired ratio
%between cash and M3 in the steady state. WARNING FOR THIS CHANGE!!!

%old equations
%equation 27
%qbar = (Aq/((1-tauk)*phiz*(Rbar-1)))^(1/sigmaq);
%equation 28
%mbar = piebar*muz*(nubar*wbar*Hbar+qbar)/mu;%Corrected 2003-12-22 by JL

%new equations
if nubar ~= 0;
   mbar = (nubar*wbar*Hbar)/(1-Aq); 
   qbar = Aq*mbar;
else;
   disp('Warning: nubar is zero, assume that mbar is a small number ...');
   mbar = 1e-6;
   qbar = Aq*mbar;
end;

if (qbar/mbar) > 1;
   error('Steady state: Warning -> MO is bigger than M3');
end;

%computing some alternative measures to verify that the steady-state
%calculations are not necessarily incorrect.

ybar1 = (muz^-sharek)*(kbar^sharek)*Hbar^(1-sharek)-phi;
ybar2 = (cdbar+cmbar) + (idbar+imbar) + gbar + xtilde - zbar;
ybar3 = (wbar*Rfbar*Hbar+rkbar*kbar*(1/muz));

if abs(ybar-ybar1) > 1e-9 | abs(ybar-ybar2) > 1e-9 | abs(ybar-ybar3) > 1e-9;
   error('The steady state calculations are not correct.');
end;

% Get ratios.

% disp('steady state ratios for various aggregates in the model - should match those in the data.');
ratioc  =(cdbar+cmbar)/ybar;
ratioi  =(idbar+imbar)/ybar;
ratiox = xtilde/ybar;
ratioz = zbar/ybar;

if abs(1-(ratiog+ratiox+ratioc+ratioi-ratioz)) > 1e-9;
   ratiog+ratiox+ratioc+ratioi-ratioz
   error('SteadyState: The sum of the GDP component ratios do not sum to unity.');
end;

mcbar= ((1/(1-sharek))^(1-sharek))*((1/sharek)^sharek)*(rkbar^sharek)*((wbar*(nubar*Rbar+1-nubar))^(1-sharek) );

if abs(mcbar-(1/lambdaf)) > 1e-6;
    error('SteadyState: Something is wrong with marginal cost in SS.');
end

profif_bar = (gammamcd - gammaf)*cmbar + (gammamid - gammaf)*imbar; 

%Check the stationarized budget contraint in SS
budgetSS = (1 - Rbar)*mbar + cbar*gammacd + ibar*gammaid - (1 - Rbar)*qbar - rkbar*kbar*(1/muz)  - wbar*Hbar + gbar ... 
           - profif_bar;

if abs(budgetSS) > 1e-6;
   error('SteadyState: Budget restriction does not hold in SS.');
end;


%Parameters to the wage equation
  bw     = (lambdaw*sigmal - (1-lambdaw))/((1-beta*xiw)*(1-xiw)) ;
  eta0   =  bw*xiw;
  eta1   = (sigmal*lambdaw - bw*(1+beta*(xiw^2)));
  eta2   =  bw*beta*xiw;
  eta3   = -bw*xiw;
  eta4   =  bw*beta*xiw;
  eta5   =  bw*kappaw*xiw;
  eta6   = -bw*beta*kappaw*xiw;
  eta7   =  (1-lambdaw);
  eta8   = -(1-lambdaw)*sigmal;
  eta9   = -(1-lambdaw)*(tauy/(1-tauy));
  eta10  = -(1-lambdaw)*(tauw/(1+tauw));
  eta11  = -(1-lambdaw);

  %%
data_pid = 400*(piebar-1)*piebar;
data_wp = 100*log(muz);
data_c = 100*log(muz);
data_i = 100*log(muz);
data_R = 400*(Rbar-1)*Rbar; 
data_y = 100*log(muz);
data_x = 100*log(muz);
data_m = 100*log(muz);
data_pic = 400*(piebar-1)*piebar;
data_pii = 400*(piebar-1)*piebar;
data_ystar = 100*log(muz);
data_pistar = 400*(piebar-1)*piebar;
data_Rstar = 400*(Rbar-1)*Rbar;
%%
ys = zeros( M_.endo_nbr, 1);
ys(55) = data_pid;
ys(56) = data_wp;
ys(57) = data_c;
ys(58) = data_i;
ys(59) = data_R;
ys(60) = data_y;
ys(61) = data_x;
ys(62) = data_m;
ys(63) = data_pic;
ys(64) = data_pii;
ys(65) = data_ystar;
ys(66) = data_pistar;
ys(67) = data_Rstar;
ys(68) = 0;

check = 0;  %Dynare Switch indicating good steady state

M_.params(	1	) =	Al 	;
M_.params(	2	 ) =	Aq 	;
M_.params(	3	 ) =	b 	;
M_.params(	4	 ) =	beta 	;
M_.params(	5	 ) =	cbar 	;
M_.params(	6	 ) =	cdbar 	;
M_.params(	7	 ) =	cdfbar 	;
M_.params(	8	 ) =	csbar 	;
M_.params(	9	 ) =	D1 	;
M_.params(	10	 ) =	D2 	;
M_.params(	11	 ) =	D3 	;
M_.params(	12	 ) =	D4 	;
M_.params(	13	 ) =	delta 	;
M_.params(	14	 ) =	epspolicy 	;
M_.params(	15	 ) =	epsz 	;
M_.params(	16	 ) =	etac 	;
M_.params(	17	 ) =	etaf 	;
M_.params(	18	 ) =	etai 	;
M_.params(	19	 ) =	etamc 	;
M_.params(	20	 ) =	etami 	;
M_.params(	21	 ) =	gammacd 	;
M_.params(	22	 ) =	gammacmc 	;
M_.params(	23	 ) =	gammadc 	;
M_.params(	24	 ) =	gammadi 	;
M_.params(	25	 ) =	gammaf 	;
M_.params(	26	 ) =	gammaid 	;
M_.params(	27	 ) =	gammaimi 	;
M_.params(	28	 ) =	gammamcd 	;
M_.params(	29	 ) =	gammamid 	;
M_.params(	30	 ) =	gammapi1 	;
M_.params(	31	 ) =	gammapi2 	;
M_.params(	32	 ) =	gammapi3 	;
M_.params(	33	 ) =	gammaq1 	;
M_.params(	34	 ) =	gammaq2 	;
M_.params(	35	 ) =	gammaq3 	;
M_.params(	36	 ) =	gammax 	;
M_.params(	37	 ) =	gammay1 	;
M_.params(	38	 ) =	gammay2 	;
M_.params(	39	 ) =	gammay3 	;
M_.params(	40	 ) =	gbar 	;
M_.params(	41	 ) =	Hbar 	;
M_.params(	42	 ) =	ibar 	;
M_.params(	43	 ) =	idbar 	;
M_.params(	44	 ) =	idfbar 	;
M_.params(	45	 ) =	isbar 	;
M_.params(	46	 ) =	kbar 	;
M_.params(	47	 ) =	koverH 	;
M_.params(	48	 ) =	lambdaf 	;
M_.params(	49	 ) =	lambdaw 	;
M_.params(	50	 ) =	mbar 	;
M_.params(	51	 ) =	mcbar 	;
M_.params(	52	 ) =	mu 	;
M_.params(	53	 ) =	muz 	;
M_.params(	54	 ) =	nubar 	;
M_.params(	55	 ) =	omegac 	;
M_.params(	56	 ) =	omegai 	;
M_.params(	57	 ) =	phi 	;
M_.params(	58	 ) =	phitilde 	;
M_.params(	59	 ) =	phiz 	;
M_.params(	60	 ) =	piebar 	;
M_.params(	61	 ) =	piecbar 	;
M_.params(	62	 ) =	Pii1 	;
M_.params(	63	 ) =	Pii2 	;
M_.params(	64	 ) =	Pii3 	;
M_.params(	65	 ) =	Piim 	;
M_.params(	66	 ) =	Piimc 	;
M_.params(	67	 ) =	PiimcPc 	;
M_.params(	68	 ) =	Piimi 	;
M_.params(	69	 ) =	PiimiPi 	;
M_.params(	70	 ) =	PiiTot 	;
M_.params(	71	 ) =	Pkprime 	;
M_.params(	72	 ) =	qbar 	;
M_.params(	73	 ) =	ratioc 	;
M_.params(	74	 ) =	ratiog 	;
M_.params(	75	 ) =	ratioi 	;
M_.params(	76	 ) =	ratiox 	;
M_.params(	77	 ) =	ratioz 	;
M_.params(	78	 ) =	Rbar 	;
M_.params(	79	 ) =	Rfbar 	;
M_.params(	80	 ) =	rhoeps 	;
M_.params(	81	 ) =	rhoetamc 	;
M_.params(	82	 ) =	rhoetami 	;
M_.params(	83	 ) =	rhoiota 	;
M_.params(	84	 ) =	rholambda 	;
M_.params(	85	 ) =	rhomuz 	;
M_.params(	86	 ) =	rhonu 	;
M_.params(	87	 ) =	rhophitilde 	;
M_.params(	88	 ) =	rhoR 	;
M_.params(	89	 ) =	rhozetac 	;
M_.params(	90	 ) =	rhozetah 	;
M_.params(	91	 ) =	rhozetaq 	;
M_.params(	92	 ) =	rkbar 	;
M_.params(	93	 ) =	sharek 	;
M_.params(	94	 ) =	sigmaa 	;
M_.params(	95	 ) =	sigmal 	;
M_.params(	96	 ) =	sigmaq 	;
M_.params(	97	 ) =	SS 	;
M_.params(	98	 ) =	tauc 	;
M_.params(	99	 ) =	tauk 	;
M_.params(	100	 ) =	tauw 	;
M_.params(	101	 ) =	tauy 	;
M_.params(	102	 ) =	wbar 	;
M_.params(	103	 ) =	xbar 	;
M_.params(	104	 ) =	xid 	;
M_.params(	105	 ) =	ximc 	;
M_.params(	106	 ) =	ximi 	;
M_.params(	107	 ) =	xiw 	;
M_.params(	108	 ) =	ybar 	;
M_.params(	109	 ) =	ybar1 	;
M_.params(	110	 ) =	ybar2 	;
M_.params(	111	 ) =	ybar3 	;
M_.params(	112	 ) =	zbar 	;
M_.params(	113	 ) =	bw 	;
M_.params(	114	 ) =	eta0 	;
M_.params(	115	 ) =	eta1 	;
M_.params(	116	 ) =	eta2 	;
M_.params(	117	 ) =	eta3 	;
M_.params(	118	 ) =	eta4 	;
M_.params(	119	 ) =	eta5 	;
M_.params(	120	 ) =	eta6 	;
M_.params(	121	 ) =	eta7 	;
M_.params(	122	 ) =	eta8 	;
M_.params(	123	 ) =	eta9 	;
M_.params(	124	 ) =	eta10 	;
M_.params(	125	 ) =	eta11 	;
M_.params(	126	 ) =	eta12 	;
M_.params(	127	 ) =	eta13 	;
M_.params(	128	 ) =	eta14 	;
M_.params(	129	 ) =	eta3bar 	;
M_.params(	130	 ) =	rhotauk 	;
M_.params(	131	 ) =	rhotauy 	;
M_.params(	132	 ) =	rhotauc 	;
M_.params(	133	 ) =	rhotauw 	;
M_.params(	134	 ) =	rhog 	;
M_.params(	135	 ) =	rhoepsR 	;
M_.params(	136	 ) =	rhopistar 	;
M_.params(	137	 ) =	rhoystar 	;
M_.params(	138	 ) =	rhoRstar 	;
M_.params(	139	 ) =	cmbar 	;
M_.params(	140	 ) =	imbar 	;
M_.params(	141	 ) =	rholambdax 	;
M_.params(	142	 ) =	rhozetatilde 	;
M_.params(	143	 ) =	xix 	;
M_.params(	144	 ) =	kappaw 	;
M_.params(	145	 ) =	kappad 	;
M_.params(	146	 ) =	rhopietarget 	;
M_.params(	147	 ) =	kappax 	;
M_.params(	148	 ) =	rdeltapi 	;
M_.params(	149	 ) =	rpie 	;
M_.params(	150	 ) =	rx 	;
M_.params(	151	 ) =	ry 	;
M_.params(	152	 ) =	rdeltay 	;
M_.params(	153	 ) =	kappamc 	;
M_.params(	154	 ) =	kappami 	;
M_.params(	155	 ) =	xie 	;
M_.params(	156	 ) =	phitildes 	;
M_.params(	157	 ) =	ForLag111 	;
M_.params(	158	 ) =	ForLag112 	;
M_.params(	159	 ) =	ForLag113 	;
M_.params(	160	 ) =	ForLag121 	;
M_.params(	161	 ) =	ForLag122 	;
M_.params(	162	 ) =	ForLag123 	;
M_.params(	163	 ) =	ForLag131 	;
M_.params(	164	 ) =	ForLag132 	;
M_.params(	165	 ) =	ForLag133 	;
M_.params(	166	 ) =	ForLag211 	;
M_.params(	167	 ) =	ForLag212 	;
M_.params(	168	 ) =	ForLag213 	;
M_.params(	169	 ) =	ForLag221 	;
M_.params(	170	 ) =	ForLag222 	;
M_.params(	171	 ) =	ForLag223 	;
M_.params(	172	 ) =	ForLag231 	;
M_.params(	173	 ) =	ForLag232 	;
M_.params(	174	 ) =	ForLag233 	;
M_.params(	175	 ) =	ForLag311 	;
M_.params(	176	 ) =	ForLag312 	;
M_.params(	177	 ) =	ForLag313 	;
M_.params(	178	 ) =	ForLag321 	;
M_.params(	179	 ) =	ForLag322 	;
M_.params(	180	 ) =	ForLag323 	;
M_.params(	181	 ) =	ForLag331 	;
M_.params(	182	 ) =	ForLag332 	;
M_.params(	183	 ) =	ForLag333 	;
M_.params(	184	 ) =	ForLag411 	;
M_.params(	185	 ) =	ForLag412 	;
M_.params(	186	 ) =	ForLag413 	;
M_.params(	187	 ) =	ForLag421 	;
M_.params(	188	 ) =	ForLag422 	;
M_.params(	189	 ) =	ForLag423 	;
M_.params(	190	 ) =	ForLag431 	;
M_.params(	191	 ) =	ForLag432 	;
M_.params(	192	 ) =	ForLag433 	;
M_.params(	193	 ) =	ForShock11 	;
M_.params(	194	 ) =	ForShock12 	;
M_.params(	195	 ) =	ForShock13 	;
M_.params(	196	 ) =	ForShock21 	;
M_.params(	197	 ) =	ForShock22 	;
M_.params(	198	 ) =	ForShock23 	;
M_.params(	199	 ) =	ForShock31 	;
M_.params(	200	 ) =	ForShock32 	;
M_.params(	201	 ) =	ForShock33 	;
M_.params(	202	 ) =	FisLag111 	;
M_.params(	203	 ) =	FisLag112 	;
M_.params(	204	 ) =	FisLag113 	;
M_.params(	205	 ) =	FisLag121 	;
M_.params(	206	 ) =	FisLag122 	;
M_.params(	207	 ) =	FisLag123 	;
M_.params(	208	 ) =	FisLag131 	;
M_.params(	209	 ) =	FisLag132 	;
M_.params(	210	 ) =	FisLag133 	;
M_.params(	211	 ) =	FisLag211 	;
M_.params(	212	 ) =	FisLag212 	;
M_.params(	213	 ) =	FisLag213 	;
M_.params(	214	 ) =	FisLag221 	;
M_.params(	215	 ) =	FisLag222 	;
M_.params(	216	 ) =	FisLag223 	;
M_.params(	217	 ) =	FisLag231 	;
M_.params(	218	 ) =	FisLag232 	;
M_.params(	219	 ) =	FisLag233 	;
M_.params(	220	 ) =	FisShock11 	;
M_.params(	221	 ) =	FisShock12 	;
M_.params(	222	 ) =	FisShock13 	;
M_.params(	223	 ) =	FisShock14 	;
M_.params(	224	 ) =	FisShock15 	;
M_.params(	225	 ) =	FisShock21 	;
M_.params(	226	 ) =	FisShock22 	;
M_.params(	227	 ) =	FisShock23 	;
M_.params(	228	 ) =	FisShock24 	;
M_.params(	229	 ) =	FisShock25 	;
M_.params(	230	 ) =	FisShock31 	;
M_.params(	231	 ) =	FisShock32 	;
M_.params(	232	 ) =	FisShock33 	;
M_.params(	233	 ) =	FisShock34 	;
M_.params(	234	 ) =	FisShock35 	;
M_.params(	235	 ) =	FisShock41 	;
M_.params(	236	 ) =	FisShock42 	;
M_.params(	237	 ) =	FisShock43 	;
M_.params(	238	 ) =	FisShock44 	;
M_.params(	239	 ) =	FisShock45 	;
M_.params(	240	 ) =	FisShock51 	;
M_.params(	241	 ) =	FisShock52 	;
M_.params(	242	 ) =	FisShock53 	;
M_.params(	243	 ) =	FisShock54 	;
M_.params(	244	 ) =	FisShock55 	;
M_.params(	245	 ) =	epstauk 	;
M_.params(	246	 ) =	epstauw	;
M_.params(	247	 ) =    lambdamc;
M_.params(	248	 ) =    lambdami;


end