%find steady state
clear all
clc

Arr={'Y' 'K' 'L' 'VARPHI' 'W' 'RK' 'PI' 'PINT' 'PK' 'DI' ...
    'DK' 'I' 'Q' 'R' 'PSI' 'VO' 'EPSILON' 'THETA' 'CO' 'AO' 'HO' ...
    'SO' 'EO' 'LO' 'TO' 'PTO' 'VY' 'CY' 'AY' 'HY' 'SY' 'OMEGA' 'LY' ... 
    'TY' 'PTY' 'C' 'A' 'LAMBDA' 'V' 'B' 'G' 'E' 'T' 'S' 'NX' ...
    'PHI' 'FSTAR'} ;
init = ones(length(Arr),1) ;
options=optimset('MaxFunEvals',5000000,'MaxIter',500000,'TolFun',1e-30,'TolX',1e-30);


F=fsolve(@(G) steady_state_function(G),init,options);
init=F ;
ss=init.^2 ;

for j=1:length(Arr)
disp([Arr{j},'=', num2str(ss(j),5), ';']) 

end