function [ys,check] = simple_steadystate(ys,exe)
global M_ lgy_

if isfield(M_,'param_nbr') == 1
NumberOfParameters = M_.param_nbr;
for i = 1:NumberOfParameters
  paramname = deblank(M_.param_names(i,:));
  eval([ paramname ' = M_.params(' int2str(i) ');']);
end
check = 0;
end

%%%% Model equations to be entered here
rstar = (1-bet)/bet;
wstar = (1-alph)*((rstar/alph)^(alph/(alph-1)));
f = @(x) (wstar*(x^(-gamm)))^(1/sigm)-rstar*x*((rstar/alph)^(1/(alph-1)))-wstar*x;
lstar = fzero(f,0.01);
cstar = (wstar*(lstar^(-gamm)))^(1/sigm);
kstar = lstar*((rstar/alph)^(1/(alph-1)));
cs = cstar;
ks = kstar;
rs = rstar;
ls = lstar;
as = 1;
ws = wstar;
c=0;
a=0;
k=0;
l=0;
r=0;
w=0;

for iter = 1:length(M_.params)
    eval([ 'M_.params(' num2str(iter) ') = ' M_.param_names(iter,:) ';' ])
end

if isfield(M_,'param_nbr') == 1

if isfield(M_,'orig_endo_nbr') == 1
NumberOfEndogenousVariables = M_.orig_endo_nbr;
else
NumberOfEndogenousVariables = M_.endo_nbr;
end
ys = zeros(NumberOfEndogenousVariables,1);
for i = 1:NumberOfEndogenousVariables
varname = deblank(M_.endo_names(i,:));
eval(['ys(' int2str(i) ') = ' varname ';']);
end
else
ys=zeros(length(lgy_),1);
for i = 1:length(lgy_)
ys(i) = eval(lgy_(i,:));
end
check = 0;
end
end