function parameter=calibration

clear;

%households
BETAH =.99;
R_ss=1/BETAH;

%SB
%S_s=.01/4;
S_s=.016/4;
PHIS_ss=20;
Rs_ss=R_ss+S_s;

Z_ss = Rs_ss - 1;%Target q = 1 and spreads = 0.016
THETA=R_ss/PHIS_ss+S_s;%Taget phis = 20
ETA_S=0.0;

%CB
S_c=.01/4;
%S_c=.008/4;
PHIC_ss=10;
CHI=0.1;

Y_ss=4*Z_ss;
ETA_C = 0.2;

kc = 1;
q  = 1;
nc = q*kc*CHI;
wc = ETA_C*Y_ss;
ALPHAC=((Z_ss+q)/R_ss-q)/kc;
dc=ALPHAC*kc^2/2+q*kc-nc;
SIGMAC=(nc-wc)/((Z_ss+q)*kc-R_ss*dc-wc);
OMEGAK=0.01;
WH=Z_ss;
EPSILON=0.2;
SIGMAS=SIGMAC;


BETAS =.9832;
BETAC =.99;


parameter =[BETAH WH S_s BETAS PHIS_ss SIGMAS ETA_S Z_ss THETA S_c BETAC PHIC_ss CHI ALPHAC ETA_C SIGMAC EPSILON Y_ss OMEGAK];
