function [ys,check] = quant_banking_exo_steadystate(ys,exe)
  global M_
  
  alpha = M_.params(1); 
  beta = M_.params(2); 
  delta_k = M_.params(3); 
  psi = M_.params(4); 
  aag = M_.params(5); 
  tau = M_.params(6); 
  eps = M_.params(7); 
  meu=M_.params(8);
  ir=M_.params(9);
  mbbar=M_.params(10);
  gamm=M_.params(11);
  rho_mu=M_.params(12);
  rho_g=M_.params(13);
  rho_mb=M_.params(14);
  

  check = 0;
  
w = ag*(1-alpha)*(kh^alpha)*(lg^(-alpha));
r = ag*alpha*lg^(1-alpha)*kh^(alpha-1);
yh = ag*kh^alpha*lg^(1-alpha);

kh = (lg*w+((1-tau)*r+1-delk)*kh-ch-gh)/(1+gamm);
(psi/(1-lg)) = (w/ch);
gamm = beta*((1-eps)*(1+(1-tau)*r-delk)/(1+beta*eps*(1+ir)/(1+meu)))-1;
bh = kh/(((1-eps)/eps)*mb-1);
gh = tau*yh + (1+meu)*(kh/((1-eps)*mb-1))*((1-beta)/beta) + (meu-ir)*mb;
 
  
  ys =[
w     
r     
yh     
kh
ch
bh     
gh     
          ];