% -------------------------------------------------------------------------
% computes the steady state analyticaly
% -------------------------------------------------------------------------

function [ys,check] = modelo_cns_steadystate(ys,exe)
  global M_
  
  kap = M_.params(1);
  del = M_.params(2);
  bet = M_.params(3);
  phi = M_.params(4);
  lam = M_.params(5);
  rho = M_.params(6);
  eta = M_.params(7);
  f1  = M_.params(8);
  f2  = M_.params(9);
  f3  = M_.params(10);
  a1  = M_.params(11);
  b1  = M_.params(12);
  
  check = 0;
  
  y  = 0;
  p  = 0;
  r  = 0;
  s  = 0;
  rr = 0;
  dq = 0;
  
  ys =[y p r rr s dq]';
