%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) ;
F0 = init ;
options = optimoptions('fsolve','Display','iter','MaxFunEvals',1000000,'MaxIter',100000); % Option to display output

X = fsolve(@myfun,F0,options) ; 
X ;

