%PURPOSE: This code computes the analytical Steady State for ict2.mod (Dynare)
%CALLED BY: ict2.mod
%AUTHOR: Benedetto Molinari, 2011


function [ys,check] = ict2_steadystate(junk,ys)


 global bet sig phi csi H alp alp_s alp_i alp_n del_s del_i del_n tg gam_i gam_n
 check=0; %needed for internal Dynare syntax

 
 %1. Normalized SS
  a = 1;
  qi= 1;
  qn= 1;
  h=H;

  %2. Analitical values of SS
  rs=(tg^sig-bet*(1-del_s))/bet;
  ri=(gam_i*tg^sig-bet*(1-del_i))/(gam_i*bet);
  rn=(gam_n*tg^sig-bet*(1-del_n))/(gam_n*bet);
  ksy=alp_s*tg/rs;
  kiy=alp_i*tg*gam_i/ri;
  kny=alp_n*tg*gam_n/rn;
  out=(ksy^alp_s*(kny/gam_n)^alp_n*(kiy/gam_i)^alp_i*tg^(alp-1))^(1/alp)*h;
  ks=ksy*out;
  kn=kny*out;
  ki=kiy*out;
  is=(1-(1-del_s)/tg)*ks;
  iv=(1-(1-del_i)/(tg*gam_i))*ki/gam_i;
  in=(1-(1-del_n)/(tg*gam_n))*kn/gam_n;
  c=out-is-in-iv;
  w=alp*out/h;
  prm=out/h;
  tinv=is+iv+in;
  dum=w+rs;

  %3. Define Steady States using dynare syntax  
  ys =[
      a
      c
      h
      in
      is
      iv
      ki
      kn
      ks
      out
      prm
      qi
      qn
      ri
      rn
      rs
      tinv
      w
      ];