diff --git a/tests/steady_state/walsh1_old_ss_steadystate.m b/tests/steady_state/walsh1_old_ss_steadystate.m index 6543d0618..fa6f5d95e 100644 --- a/tests/steady_state/walsh1_old_ss_steadystate.m +++ b/tests/steady_state/walsh1_old_ss_steadystate.m @@ -1,5 +1,7 @@ function [ys,params,check] = walsh1_old_ss_steadystate(ys,exo,M_,options_) +delta = NaN; + % read out parameters to access them with their name NumberOfParameters = M_.param_nbr; for ii = 1:NumberOfParameters @@ -11,35 +13,35 @@ check = 0; %% Enter model equations here - - pi = thetass-1; - en = 1/3; - eR = 1/betta; - y_k = (1/alphha)*(1/betta-1+delta); - ek = en*y_k^(-1/(1-alphha)); - ec = ek*(y_k-delta); - em = ec*(a/(1-a))^(-1/b)*((thetass-betta)/thetass)^(-1/b); - ey = ek*y_k; - Xss = a*ec^(1-b)*(1+(a/(1-a))^(-1/b)*((thetass-betta)/thetass)^((b-1)/b)); - Psi = (1-alphha)*(ey/en)*Xss^((b-phi1)/(1-b))*a*ec^(-b)*(1-en)^eta; - n = log(en); - k = log(ek); - m = log(em); - c = log(ec); - y = log(ey); - R = log(eR); - z = 0; - u = 0; - + +pi = thetass-1; +en = 1/3; +eR = 1/betta; +y_k = (1/alphha)*(1/betta-1+delta); +ek = en*y_k^(-1/(1-alphha)); +ec = ek*(y_k-delta); +em = ec*(a/(1-a))^(-1/b)*((thetass-betta)/thetass)^(-1/b); +ey = ek*y_k; +Xss = a*ec^(1-b)*(1+(a/(1-a))^(-1/b)*((thetass-betta)/thetass)^((b-1)/b)); +Psi = (1-alphha)*(ey/en)*Xss^((b-phi1)/(1-b))*a*ec^(-b)*(1-en)^eta; +n = log(en); +k = log(ek); +m = log(em); +c = log(ec); +y = log(ey); +R = log(eR); +z = 0; +u = 0; + %% end own model equations params=NaN(NumberOfParameters,1); for iter = 1:NumberOfParameters %update parameters set in the file - eval([ 'params(' num2str(iter) ') = ' M_.param_names{iter} ';' ]) + eval([ 'params(' num2str(iter) ') = ' M_.param_names{iter} ';' ]) end NumberOfEndogenousVariables = M_.orig_endo_nbr; %auxiliary variables are set automatically for ii = 1:NumberOfEndogenousVariables - varname = M_.endo_names{ii}; - eval(['ys(' int2str(ii) ') = ' varname ';']); + varname = M_.endo_names{ii}; + eval(['ys(' int2str(ii) ') = ' varname ';']); end