
function [ys,check] = model_steadystate(ys,exe)
  global M_
  
  beta = M_.params(1); 
  theta = M_.params(2); 
  gamma = M_.params(3); 
  alpha = M_.params(4); 
  rho = M_.params(5); 
  phi = M_.params(6); 
  rho1 = M_.params(7); 
  omega=M_.params(8);
  
  check = 0;
  
eta=-1/theta;
kappa=1;
a=1;
i=(1-beta)/beta;
pa=1;
k=((1+eta)*(theta/(theta-1))*kappa)^((gamma-1)/(gamma+phi))*(1/(1-alpha*beta))*(a/(kappa)^(1/(1+phi)))^((1-gamma)*(1+phi)/(gamma+phi));
f=((1+eta)*(theta/(theta-1))*kappa)^((gamma-1)/(gamma+phi))*(1/(1-alpha*beta))*(a/(kappa)^(1/(1+phi)))^((1-gamma)*(1+phi)/(gamma+phi));
l=((1+eta)*(theta/(theta-1))*kappa)^(-1/(gamma+phi))*(a)^((1-gamma)/(phi+gamma));
c=((1+eta)*(theta/(theta-1))*kappa)^(-1/(gamma+phi))*(a)^((phi+1)/(phi+gamma));
pi=1;
alphat=alpha;
  
  
  ys =[
a   
alphat     
c    
eta     
f  
i     
k     
kappa     
l     
pa
pi        ];
