var ki k i y c psi;

parameters beta h delta deltai theta phi q;

beta = 0.979;
h = 0.3;
delta = 0.039;
deltai = 0.039;
theta = 0.212;
phi = 0.150;
q = 1;

model; 
(1/c) = beta*(1/c(+1))*(1 - delta + theta*(k^(theta-1)*(ki^phi)*(h(+1)^(1-theta-phi))));
(1/c)*q = beta*(1/c(+1))*(1 - deltai + phi*(k^theta)*(ki^(phi-1))*(h(+1)^(1-theta-phi)));
psi*c/(1-h) = ((1-theta-phi)*(k(-1)^theta)*(ki(-1)^phi))/(h^(theta+phi));
c = y-i;
i = k+ki*q+(1-delta+theta*(k^(theta-1)*(ki^phi)*(h^(1-theta-phi)))*(k(-1))+(1-deltai+(phi*(k^theta)*(ki^(phi-1))*(h^(1-theta-phi))))*(ki(-1)));
y = (k(-1)^theta)*(ki(-1)^phi)*(h^(1-theta-phi));

end;

initval;
ki =steady(beta, h, delta, deltai, theta, phi);
k =steadyk(beta, h, ki, delta, theta, phi);
i =steadyi(beta, h, ki, delta, deltai, theta, phi);
y =steadyy(k, theta, ki, phi, h);
c =steadyc(y, i);
psi =Steadypsi(c, k, ki, h, theta, phi);
end;

resid;