function F = Estimation_EtatStationnaire(x,tk,tw,tc,tpi,nivby,mu_A,A0,alpha,b,delta_g,delta_p,phi,rho_h,rho_c,rho_a,rho_g,mu,mu_ikg,rho_ikg,nu,beta,r,theta,zeta,GBS_Y_SS)

    
    y     = x(1); 
    w     = x(2);
    H     = x(3); 
    k     = x(4);
    kp    = x(5);
    c     = x(6);
    ip    = x(7);
    kg    = x(8);
    ig    = x(9);
    pi    = x(10);    
    dette = x(11);
    gbs   = x(12);

    
    F = [ w - (1-alpha)*y/H;
          ((exp(mu_A))^theta/beta - 1 + delta_p) + tk - alpha*(beta^(1+phi))*(y/k)*((k/kp)^(1+phi))*(exp(mu_A));
          zeta*(1+1/nu)*(H)^(1/nu) - (w/c)*((1 - tw)/(1 + tc));
          y/H - ((k/(H*(exp(mu_A))))^alpha)*(A0^(1 - alpha));
          k^(-phi) - ((b^(1 + phi))*kp^(-phi) + ((1-b)^(1 + phi))*kg^(-phi));
          kg  - (1-delta_g)*kg/(exp(mu_A)) - ig;
          kp  - (1-delta_p)*kp/(exp(mu_A)) - ip;
          ig  - mu_ikg*y ;
          dette/y - 4*nivby;
          y - c - ip - ig - gbs;
          pi - (y - w*H - (((exp(mu_A))^theta/beta - 1 + delta_p) + tk)*kp/(exp(mu_A)))*(1 - tpi);
          dette/(exp(mu_A)) + pi + ((exp(mu_A))^theta/beta - 1 + delta_p)*kp/(exp(mu_A)) + H*w*(1 - tw) - c*(1 + tc) - dette/(1 + (((exp(mu_A))^theta/beta - 1 + delta_p) - delta_p)) - ip ];
             
  
 	
  	