%
% Status : main Dynare file 
%
% Warning : this file is generated automatically by Dynare
%           from model file (.mod)

clear all
tic;
global M_ oo_ options_ ys0_ ex0_
options_ = [];
M_.fname = 'ParamValue';
%
% Some global variables initialization
%
global_initialization;
diary off;
logname_ = 'ParamValue.log';
if exist(logname_, 'file')
    delete(logname_)
end
diary(logname_)
M_.param_names = 'bc';
M_.param_names_tex = 'bc';
M_.param_names = strvcat(M_.param_names, 'sigmac');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmac');
M_.param_names = strvcat(M_.param_names, 'betaa');
M_.param_names_tex = strvcat(M_.param_names_tex, 'betaa');
M_.param_names = strvcat(M_.param_names, 'alpha');
M_.param_names_tex = strvcat(M_.param_names_tex, 'alpha');
M_.param_names = strvcat(M_.param_names, 'Omega');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Omega');
M_.param_names = strvcat(M_.param_names, 'alphap');
M_.param_names_tex = strvcat(M_.param_names_tex, 'alphap');
M_.param_names = strvcat(M_.param_names, 'tetap');
M_.param_names_tex = strvcat(M_.param_names_tex, 'tetap');
M_.param_names = strvcat(M_.param_names, 'mup');
M_.param_names_tex = strvcat(M_.param_names_tex, 'mup');
M_.param_names = strvcat(M_.param_names, 'gammap');
M_.param_names_tex = strvcat(M_.param_names_tex, 'gammap');
M_.param_names = strvcat(M_.param_names, 'muW');
M_.param_names_tex = strvcat(M_.param_names_tex, 'muW');
M_.param_names = strvcat(M_.param_names, 'mu');
M_.param_names_tex = strvcat(M_.param_names_tex, 'mu');
M_.param_names = strvcat(M_.param_names, 'varrho');
M_.param_names_tex = strvcat(M_.param_names_tex, 'varrho');
M_.param_names = strvcat(M_.param_names, 'rhoi');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoi');
M_.param_names = strvcat(M_.param_names, 'ap');
M_.param_names_tex = strvcat(M_.param_names_tex, 'ap');
M_.param_names = strvcat(M_.param_names, 'ay');
M_.param_names_tex = strvcat(M_.param_names_tex, 'ay');
M_.param_names = strvcat(M_.param_names, 'eps_tax');
M_.param_names_tex = strvcat(M_.param_names_tex, 'eps\_tax');
M_.param_names = strvcat(M_.param_names, 'delta');
M_.param_names_tex = strvcat(M_.param_names_tex, 'delta');
M_.param_names = strvcat(M_.param_names, 'd');
M_.param_names_tex = strvcat(M_.param_names_tex, 'd');
M_.param_names = strvcat(M_.param_names, 'kapai');
M_.param_names_tex = strvcat(M_.param_names_tex, 'kapai');
M_.param_names = strvcat(M_.param_names, 'sigmal');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmal');
M_.param_names = strvcat(M_.param_names, 'sigmaWbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaWbar');
M_.param_names = strvcat(M_.param_names, 'eps_taxw');
M_.param_names_tex = strvcat(M_.param_names_tex, 'eps\_taxw');
M_.param_names = strvcat(M_.param_names, 'muww');
M_.param_names_tex = strvcat(M_.param_names_tex, 'muww');
M_.param_names = strvcat(M_.param_names, 'alphaw');
M_.param_names_tex = strvcat(M_.param_names_tex, 'alphaw');
M_.param_names = strvcat(M_.param_names, 'tetaw');
M_.param_names_tex = strvcat(M_.param_names_tex, 'tetaw');
M_.param_names = strvcat(M_.param_names, 'gammaw');
M_.param_names_tex = strvcat(M_.param_names_tex, 'gammaw');
M_.param_names = strvcat(M_.param_names, 'Q1wbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Q1wbar');
M_.param_names = strvcat(M_.param_names, 'Q2wbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Q2wbar');
M_.param_names = strvcat(M_.param_names, 'Wstarbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Wstarbar');
M_.param_names = strvcat(M_.param_names, 'lomebar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'lomebar');
M_.param_names = strvcat(M_.param_names, 'muGpbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'muGpbar');
M_.param_names = strvcat(M_.param_names, 'pi');
M_.param_names_tex = strvcat(M_.param_names_tex, 'pi');
M_.param_names = strvcat(M_.param_names, 'SEome');
M_.param_names_tex = strvcat(M_.param_names_tex, 'SEome');
M_.param_names = strvcat(M_.param_names, 'lambdabar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'lambdabar');
M_.param_names = strvcat(M_.param_names, 'cbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'cbar');
M_.param_names = strvcat(M_.param_names, 'Rbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Rbar');
M_.param_names = strvcat(M_.param_names, 'pipbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'pipbar');
M_.param_names = strvcat(M_.param_names, 'rrealbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rrealbar');
M_.param_names = strvcat(M_.param_names, 'lhbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'lhbar');
M_.param_names = strvcat(M_.param_names, 'wbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'wbar');
M_.param_names = strvcat(M_.param_names, 'piwbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'piwbar');
M_.param_names = strvcat(M_.param_names, 'ybar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'ybar');
M_.param_names = strvcat(M_.param_names, 'lbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'lbar');
M_.param_names = strvcat(M_.param_names, 'kbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'kbar');
M_.param_names = strvcat(M_.param_names, 'sbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sbar');
M_.param_names = strvcat(M_.param_names, 'webar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'webar');
M_.param_names = strvcat(M_.param_names, 'zbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'zbar');
M_.param_names = strvcat(M_.param_names, 'Q1bar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Q1bar');
M_.param_names = strvcat(M_.param_names, 'Q2bar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Q2bar');
M_.param_names = strvcat(M_.param_names, 'Pstarbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Pstarbar');
M_.param_names = strvcat(M_.param_names, 'Fomebar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Fomebar');
M_.param_names = strvcat(M_.param_names, 'abar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'abar');
M_.param_names = strvcat(M_.param_names, 'rhoA');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoA');
M_.param_names = strvcat(M_.param_names, 'sigmaA');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaA');
M_.param_names = strvcat(M_.param_names, 'Fnomebar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'Fnomebar');
M_.param_names = strvcat(M_.param_names, 'omebar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'omebar');
M_.param_names = strvcat(M_.param_names, 'SIGbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'SIGbar');
M_.param_names = strvcat(M_.param_names, 'muGbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'muGbar');
M_.param_names = strvcat(M_.param_names, 'SIGpbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'SIGpbar');
M_.param_names = strvcat(M_.param_names, 'ffomebar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'ffomebar');
M_.param_names = strvcat(M_.param_names, 'deltaubar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'deltaubar');
M_.param_names = strvcat(M_.param_names, 'deltauubar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'deltauubar');
M_.param_names = strvcat(M_.param_names, 'adjibar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'adjibar');
M_.param_names = strvcat(M_.param_names, 'adjiibar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'adjiibar');
M_.param_names = strvcat(M_.param_names, 'nwbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'nwbar');
M_.param_names = strvcat(M_.param_names, 'ibar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'ibar');
M_.param_names = strvcat(M_.param_names, 'ubar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'ubar');
M_.param_names = strvcat(M_.param_names, 'gbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'gbar');
M_.param_names = strvcat(M_.param_names, 'chilbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'chilbar');
M_.param_names = strvcat(M_.param_names, 'rhoChil');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoChil');
M_.param_names = strvcat(M_.param_names, 'rhoG');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoG');
M_.param_names = strvcat(M_.param_names, 'ktildebar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'ktildebar');
M_.param_names = strvcat(M_.param_names, 'bbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'bbar');
M_.param_names = strvcat(M_.param_names, 'xxbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'xxbar');
M_.param_names = strvcat(M_.param_names, 'qbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'qbar');
M_.param_names = strvcat(M_.param_names, 'rkbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rkbar');
M_.param_names = strvcat(M_.param_names, 'rrbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rrbar');
M_.param_names = strvcat(M_.param_names, 'LAMbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'LAMbar');
M_.param_names = strvcat(M_.param_names, 'vbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'vbar');
M_.param_names = strvcat(M_.param_names, 'cebar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'cebar');
M_.param_names = strvcat(M_.param_names, 'gammabar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'gammabar');
M_.param_names = strvcat(M_.param_names, 'deltatild');
M_.param_names_tex = strvcat(M_.param_names_tex, 'deltatild');
M_.param_names = strvcat(M_.param_names, 'rhoGamma');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoGamma');
M_.param_names = strvcat(M_.param_names, 'rhoPref');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoPref');
M_.param_names = strvcat(M_.param_names, 'rhoEm');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoEm');
M_.param_names = strvcat(M_.param_names, 'rhoTechno');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoTechno');
M_.param_names = strvcat(M_.param_names, 'prefbar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'prefbar');
M_.param_names = strvcat(M_.param_names, 'embar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'embar');
M_.param_names = strvcat(M_.param_names, 'technobar');
M_.param_names_tex = strvcat(M_.param_names_tex, 'technobar');
M_.param_names = strvcat(M_.param_names, 'sigmaGamma');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaGamma');
M_.param_names = strvcat(M_.param_names, 'sigmaPref');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaPref');
M_.param_names = strvcat(M_.param_names, 'sigmaTechno');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaTechno');
M_.param_names = strvcat(M_.param_names, 'sigmaEm');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaEm');
M_.param_names = strvcat(M_.param_names, 'sigmaChil');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaChil');
M_.param_names = strvcat(M_.param_names, 'sigmaG');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaG');
M_.param_names = strvcat(M_.param_names, 'sigmaSIG');
M_.param_names_tex = strvcat(M_.param_names_tex, 'sigmaSIG');
M_.param_names = strvcat(M_.param_names, 'rhoSIG');
M_.param_names_tex = strvcat(M_.param_names_tex, 'rhoSIG');
M_.exo_det_nbr = 0;
M_.exo_nbr = 0;
M_.endo_nbr = 0;
M_.param_nbr = 97;
M_.Sigma_e = zeros(0, 0);
M_.orig_endo_nbr = 0;
M_.aux_vars = [];
options_.block=0;
options_.bytecode=0;
options_.use_dll=0;
erase_compiled_function('ParamValue_dynamic');
M_.params( 3 ) = 0.99;
betaa = M_.params( 3 );
M_.params( 1 ) = 0;
bc = M_.params( 1 );
M_.params( 2 ) = 1;
sigmac = M_.params( 2 );
M_.params( 20 ) = 1;
sigmal = M_.params( 20 );
M_.params( 4 ) = 0.35;
alpha = M_.params( 4 );
M_.params( 5 ) = 0.9846;
Omega = M_.params( 5 );
M_.params( 17 ) = 0.025;
delta = M_.params( 17 );
tau     = 0;        % tax (does not appear in .mod)
M_.params( 19 ) = 5.86;
kapai = M_.params( 19 );
M_.params( 6 ) = 0.67;
alphap = M_.params( 6 );
M_.params( 9 ) = 0.75;
gammap = M_.params( 9 );
M_.params( 7 ) = 11;
tetap = M_.params( 7 );
M_.params( 8 ) = M_.params(7)/(M_.params(7)-1);
mup = M_.params( 8 );
SNbar    = 1/mup;
M_.params( 24 ) = 0.75;
alphaw = M_.params( 24 );
M_.params( 26 ) = 0.75;
gammaw = M_.params( 26 );
M_.params( 25 ) = 11;
tetaw = M_.params( 25 );
M_.params( 23 ) = M_.params(25)/(M_.params(25)-1);
muww = M_.params( 23 );
kappap  = (1-alphap)*(1-alphap*betaa)/(alphap);                          
kappaw  = ((1-alphaw)*(1-betaa*alphaw)/(alphaw*(1+sigmal*tetaw)));
xbar    = 2;            % = ktildebar/nbar
M_.params( 77 ) = 1.004962931573204;
rrbar = M_.params( 77 );
M_.params( 36 ) = 1/M_.params(3);
Rbar = M_.params( 36 );
M_.params( 76 ) = M_.params(77)*M_.params(36);
rkbar = M_.params( 76 );
M_.params( 81 ) = 0.9728;
gammabar = M_.params( 81 );
M_.params( 12 ) = 0.01;
varrho = M_.params( 12 );
Eome     = 1;       % mean of omegabar -- BGG
omebar0  = 1;       % Initial value of omegabar 
SEome0   = 0.28;    % s.e of omegabar -- BGG
gamma0   = 0.9;
mu0      = 0.1;
paramm0  = [omebar0 SEome0 gamma0 mu0];
paramc   = [tau betaa delta alpha Omega tetap xbar rrbar Eome];
options  = optimset('Display','iter');
[x,Fval,exitflag] = fsolve('funx2',paramm0,options,paramc);      % what value of omebar such that: logncdf(omebar,muW,sigmaW) = PHI;
M_.params( 56 ) = x(1);
omebar = M_.params( 56 );
M_.params( 33 ) = x(2);
SEome = M_.params( 33 );
M_.params( 81 ) = x(3);
gammabar = M_.params( 81 );
M_.params( 11 ) = x(4);
mu = M_.params( 11 );
M_.params( 21 ) = sqrt(log(1+M_.params(33)^2/Eome^2));
sigmaWbar = M_.params( 21 );
M_.params( 10 ) = log(Eome)-0.5*M_.params(21)^2;
muW = M_.params( 10 );
tbar       = (muW - log(omebar))/sigmaWbar;
M_.params( 30 ) = log(M_.params(56));
lomebar = M_.params( 30 );
M_.params( 51 ) = normcdf(M_.params(30),M_.params(10),M_.params(21));
Fomebar = M_.params( 51 );
aa = exp(-0.5 * ((log(omebar) - muW)/sigmaWbar)^2) / (omebar * sqrt(2*pi) * sigmaWbar);
M_.params( 60 ) = aa;
ffomebar = M_.params( 60 );
ffomep    = -(ffomebar/omebar)*(1 - (muW-log(omebar))/(sigmaWbar*sigmaWbar));  
M_.params( 55 ) = Eome-exp(M_.params(10)+M_.params(21)*M_.params(21)/2)*normcdf(M_.params(21)+tbar,0,1);
Fnomebar = M_.params( 55 );
M_.params( 57 ) = M_.params(56)*(1-M_.params(51))+M_.params(55);
SIGbar = M_.params( 57 );
M_.params( 59 ) = 1-M_.params(51);
SIGpbar = M_.params( 59 );
M_.params( 58 ) = M_.params(55)*M_.params(11);
muGbar = M_.params( 58 );
M_.params( 31 ) = M_.params(56)*M_.params(11)*M_.params(60);
muGpbar = M_.params( 31 );
SIGppbar  = - ffomebar;                       % SIGMA''(omebar)
M_.params( 58 ) = M_.params(55)*M_.params(11);
muGbar = M_.params( 58 );
M_.params( 31 ) = M_.params(56)*M_.params(11)*M_.params(60);
muGpbar = M_.params( 31 );
muGppbar  = mu*(ffomebar + omebar*ffomep);     % mu*G''(omebar)
vareps   = 1e-15;
SEome1   = SEome + vareps;
sigmaW1  = sqrt(log(1 + (SEome1^2)/(Eome^2)));  
muW1     = log(Eome) - 0.5*sigmaWbar^2;
sigmaW   = sigmaW1;
M_.params( 10 ) = muW1;
muW = M_.params( 10 );
tbar1    = (muW1 - log(omebar))/sigmaW1;
Fome1    = logncdf(omebar,muW1,sigmaW1);   % F(omebar)
fome1    = lognpdf(omebar,muW1,sigmaW1);   % f(omebar)
fomep1   = -(fome1/omebar)*(1 - (muW1-log(omebar))/(sigmaW1*sigmaW1));   % f'(omebar)
Fnome1   = Eome - exp(muW1+sigmaW1*sigmaW1/2)*normcdf(tbar1+sigmaW1,0,1); % Fnormal(omebar) % Gamma
SIG1     = omebar*(1-Fome1) + Fnome1;      % SIGMA(omebar)
SIGp1    = 1 - Fome1;                      % SIGMA'(omebar)
SIGpp1   = - fome1;                        % SIGMA''(omebar)
muG1     = mu*Fnome1;                      % mu*G(omebar)
muGp1    = mu*omebar*fome1;                % mu*G'(omebar)
muGpp1   = mu*(fome1 + omebar*fomep1);     % mu*G''(omebar)
SIGsbar     = (SIG1-SIGbar)/vareps;
SIGpsbar    = (SIGp1-SIGpbar)/vareps;
muGsbar     = (muG1-muGbar)/vareps;
muGpsbar    = (muGp1-muGpbar)/vareps;
GbarY      = 0.19;         
M_.params( 16 ) = 1/(M_.params(7)-1);
eps_tax = M_.params( 16 );
M_.params( 22 ) = 1/(M_.params(25)-1);
eps_taxw = M_.params( 22 );
M_.params( 50 ) = 1^(1/(1-M_.params(7)));
Pstarbar = M_.params( 50 );
M_.params( 45 ) = (1+M_.params(16))*M_.params(50)/M_.params(8);
sbar = M_.params( 45 );
M_.params( 47 ) = M_.params(76)-(1-M_.params(17));
zbar = M_.params( 47 );
YbarK   = zbar/(sbar*alpha);
IbarY   = delta*(1/YbarK);
CebarK  = varrho*(1-gammabar)*(rkbar - Rbar*(1-1/xbar) - muGbar*rkbar);
CebarY  = CebarK/YbarK;
CbarY   = 1 - IbarY - CebarY - GbarY  - muGbar*rkbar*(1/YbarK) ;
WebarK  = sbar*(1-Omega)*(1-alpha)*YbarK;
VbarN   = rkbar*xbar - Rbar*(xbar - 1) - muGbar*rkbar*xbar;
CebarN  = (1-gammabar)*varrho*VbarN;
LbarK   = (YbarK)^(1/(1-alpha));
M_.params( 42 ) = 1;
ybar = M_.params( 42 );
M_.params( 41 ) = 1;
piwbar = M_.params( 41 );
M_.params( 37 ) = 1;
pipbar = M_.params( 37 );
M_.params( 38 ) = M_.params(36);
rrealbar = M_.params( 38 );
M_.params( 67 ) = 0;
ubar = M_.params( 67 );
M_.params( 75 ) = 1;
qbar = M_.params( 75 );
M_.params( 44 ) = 1/(YbarK*1/M_.params(42));
kbar = M_.params( 44 );
M_.params( 72 ) = M_.params(44);
ktildebar = M_.params( 72 );
M_.params( 66 ) = M_.params(42)*IbarY;
ibar = M_.params( 66 );
M_.params( 80 ) = M_.params(42)*CebarY;
cebar = M_.params( 80 );
M_.params( 35 ) = M_.params(42)*CbarY;
cbar = M_.params( 35 );
M_.params( 68 ) = M_.params(42)*GbarY;
gbar = M_.params( 68 );
M_.params( 65 ) = M_.params(44)/xbar;
nwbar = M_.params( 65 );
M_.params( 79 ) = VbarN*M_.params(65);
vbar = M_.params( 79 );
M_.params( 34 ) = (1-M_.params(3)*M_.params(1))*(M_.params(35)*(1-M_.params(1)))^((-1)/M_.params(2));
lambdabar = M_.params( 34 );
M_.params( 43 ) = M_.params(44)*LbarK;
lbar = M_.params( 43 );
M_.params( 39 ) = M_.params(43)^(1/M_.params(5));
lhbar = M_.params( 39 );
M_.params( 46 ) = M_.params(44)*WebarK;
webar = M_.params( 46 );
M_.params( 40 ) = M_.params(42)*M_.params(5)*M_.params(45)*(1-M_.params(4))/M_.params(39);
wbar = M_.params( 40 );
M_.params( 29 ) = M_.params(40);
Wstarbar = M_.params( 29 );
M_.params( 69 ) = (1+M_.params(22))*M_.params(29)*M_.params(34)*M_.params(39)^(-M_.params(20))/M_.params(23);
chilbar = M_.params( 69 );
M_.params( 73 ) = M_.params(72)-M_.params(65);
bbar = M_.params( 73 );
M_.params( 74 ) = M_.params(72)/M_.params(65);
xxbar = M_.params( 74 );
M_.params( 78 ) = M_.params(59)/(M_.params(59)-M_.params(31));
LAMbar = M_.params( 78 );
M_.params( 52 ) = M_.params(76)/M_.params(38)*M_.params(74)*(M_.params(57)-M_.params(58))/(M_.params(74)-1);
abar = M_.params( 52 );
M_.params( 82 ) = M_.params(47);
deltatild = M_.params( 82 );
M_.params( 18 ) = M_.params(82)/M_.params(17);
d = M_.params( 18 );
M_.params( 48 ) = M_.params(42)*M_.params(45)/(1-M_.params(3)*M_.params(6));
Q1bar = M_.params( 48 );
M_.params( 49 ) = (1+M_.params(16))*M_.params(42)/(1-M_.params(3)*M_.params(6));
Q2bar = M_.params( 49 );
M_.params( 27 ) = M_.params(69)*M_.params(39)^M_.params(20)/(1-M_.params(3)*M_.params(24));
Q1wbar = M_.params( 27 );
M_.params( 28 ) = (1+M_.params(22))*M_.params(34)/(1-M_.params(3)*M_.params(24));
Q2wbar = M_.params( 28 );
M_.params( 61 ) = M_.params(17);
deltaubar = M_.params( 61 );
M_.params( 62 ) = 0;
deltauubar = M_.params( 62 );
M_.params( 63 ) = 0;
adjibar = M_.params( 63 );
M_.params( 64 ) = 0;
adjiibar = M_.params( 64 );
M_.params( 87 ) = 1;
prefbar = M_.params( 87 );
M_.params( 88 ) = 1;
embar = M_.params( 88 );
M_.params( 89 ) = 1;
technobar = M_.params( 89 );
M_.params( 13 ) = 0.878;
rhoi = M_.params( 13 );
M_.params( 14 ) = 1.852;
ap = M_.params( 14 );
M_.params( 15 ) = 0.07825;
ay = M_.params( 15 );
aq       = 0;
M_.params( 83 ) = 0.561;
rhoGamma = M_.params( 83 );
M_.params( 84 ) = 0.903;
rhoPref = M_.params( 84 );
M_.params( 86 ) = 0.945;
rhoTechno = M_.params( 86 );
M_.params( 85 ) = 0;
rhoEm = M_.params( 85 );
M_.params( 70 ) = 0.2;
rhoChil = M_.params( 70 );
M_.params( 71 ) = 0.6;
rhoG = M_.params( 71 );
M_.params( 97 ) = 0.98;
rhoSIG = M_.params( 97 );
M_.params( 53 ) = 0.5;
rhoA = M_.params( 53 );
M_.params( 90 ) = 0.0001;
sigmaGamma = M_.params( 90 );
M_.params( 91 ) = 0.0001;
sigmaPref = M_.params( 91 );
M_.params( 92 ) = 0.0001;
sigmaTechno = M_.params( 92 );
M_.params( 93 ) = 0.0001;
sigmaEm = M_.params( 93 );
M_.params( 94 ) = 0.0001;
sigmaChil = M_.params( 94 );
M_.params( 95 ) = 0.0001;
sigmaG = M_.params( 95 );
M_.params( 96 ) = 0.0001;
sigmaSIG = M_.params( 96 );
M_.params( 54 ) = 0.0001;
sigmaA = M_.params( 54 );
save_params_and_steady_state('ParamMatrix');
save('ParamValue_results.mat', 'oo_', 'M_', 'options_');
diary off

disp(['Total computing time : ' dynsec2hms(toc) ]);
