Commit 0a3545f2 authored by MichelJuillard's avatar MichelJuillard
Browse files

adding missing test files

parent be97d37c
function [ys_, check_] = expectation_ss_old_steadystate(ys_orig_, exo_)
% Steady state generated by Dynare preprocessor
ys_=zeros(6,1);
global M_
ys_(4)=0;
ys_(6)=0;
ys_(5)=0.3333333333333333;
ys_(3)=((1/M_.params(1)-(1-M_.params(4)))/(M_.params(3)*ys_(5)^(1-M_.params(3))))^(1/(M_.params(3)-1));
ys_(1)=ys_(5)^(1-M_.params(3))*ys_(3)^M_.params(3);
ys_(2)=ys_(1)-M_.params(4)*ys_(3);
M_.params(5)=(1-M_.params(3))*ys_(1)/(ys_(2)*ys_(5)^(1+M_.params(6)));
check_=0;
end
function [ys,check] = wlash1_old_ss_steadystate(ys0,exo)
global M_
check = 0;
params = M_.params;
alpha = params(1);
beta = params(2);
delta = params(3);
gamm = params(4);
phi1 = params(5);
eta = params(6);
a = params(7);
b = params(8);
rho = params(9);
phi2 = params(10);
Psi = params(11);
thetass = params(12);
pi = thetass-1;
en = 1/3;
eR = 1/beta;
y_k = (1/alpha)*(1/beta-1+delta);
ek = en*y_k^(-1/(1-alpha));
ec = ek*(y_k-delta);
em = ec*(a/(1-a))^(-1/b)*((thetass-beta)/thetass)^(-1/b);
ey = ek*y_k;
Xss = a*ec^(1-b)*(1+(a/(1-a))^(-1/b)*((thetass-beta)/thetass)^((b-1)/b));
Psi = (1-alpha)*(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;
ys = [y c k m n R pi z u]';
M_.params(11) = Psi;
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment