ssigma=2;
ppsi=1;
oomegac = 0.5;
tthetac = 0.5;
ggamma = 0.5;
pphi = 3;
oomegaell = 0.5;
tthetaell = 10;
bbeta = 0.96;
xxib = 1;
eeta = 2;
AH = 1;
aalphah = 0.5;
xxi = 0.01;
rrho = 0.02;
Bbar = 0;
AN = 1;
aalphan = 0.5;
kkappa=0.7;

x0 = zeros(12,1);
x = fsolve(@(x)steadystate(x,ssigma,ppsi,oomegac,tthetac,ggamma,pphi,oomegaell,tthetaell,bbeta,xxib,eeta,AH,aalphah,xxi,rrho,Bbar,AN,aalphan),x0);

Css = exp(x(1));
CTss = exp(x(2));
CNss = exp(x(3));
Lss = exp(x(4));
LHss = exp(x(5));
LNss = exp(x(6));
Bsss = exp(x(7));
WPss = exp(x(8));
PEss = exp(x(9));
PPNss = exp(x(10));
WWHss = exp(x(11));
WWNss = exp(x(12));

CFss = ggamma*CTss;
CHss = (1-ggamma)*CTss;
Rss = 1/bbeta;
isss=1/bbeta-1-rrho;
Bss=Bbar;
PHPFss=1;
