% % 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) ]);