function F = Equilibrio(x)
global M_ %  beta alfa delta_k delta gamma lambda rho_z mi d omega H eta;


%  uso gli stessi valori dal file .mod (paper)
beta1        =M_.params(1);  %  beta     = 0.99;
alfa1        =M_.params(2);  %  alfa     = 0.69;
delta_k1     =M_.params(3);  %  delta_k  = 0.018;
delta1       =M_.params(4);  %  delta    = 0.095;
gamma1       =M_.params(5);  %  gamma    = 0.06;
lambda1      =M_.params(6);  %  lambda   = 0.7;
rho_z1       =M_.params(7);  %  rho_z    = 0.95;
mi1          =M_.params(8);  %  mi       = 0.5;
d1           =M_.params(9);  %  d        = 0.4;
omega1       =M_.params(10); %  omega      = 0.29;
H1           =M_.params(11); %  H        = 1.05;
eta1         =M_.params(12); %  eta      = 0.8;



%%%%%%%%%%%%%%% EQUAZIONI ORIGINALI %%%%%%%%%%%%%%

%  USO IL MODELLO NON LINEARIZZATO

%  cambio n con N  in tutte le formnule
%  N=(1-delta)*(N(-1)+m(-1));
%  c=beta*(c(+1)*(r(+1)+(1-delta_k)));
%  k=i+(1-delta_k)*k(-1);
%  m=eta*(v^lambda)*u^(1-lambda);
%  q=m/v;
%  p=m/u;
%  //u*p=v*q;
%  theta=v/u;
%  y=((N^alfa)*(k(-1)^(1-alfa)))*z;
%  y=C+i+gamma*v;
%  gamma/q=beta*(1-delta)*(c(+1)/c)*(alfa*(y(+1)/N(+1))-w(+1)+gamma/q(+1));
%  ((l-c*d)/p)=beta*(1-delta)*(c(+1)*w(+1)-l(+1)+((l(+1)-c(+1)*d)/p(+1)));
%  w=mi*(alfa*(y/N))+(1-mi)*(l/c);
%  r=(1-alfa)*(y/k(-1));
%  c=1/C;
%  L=1-u-N;
%  l=H*(L^(-omega));
%  z=rho_z*z(-1)+e;


%  USO IL MODELLO NON LINEARIZZATO SENZA INDICI TEMPORALI

F = [-x(2)+(1-delta1)*(x(2)+x(8));
-x(1)+beta1*(x(1)*(x(17)+(1-delta_k1)));
-x(4)+x(6)+(1-delta_k1)*x(4);
-x(8)+eta1*(x(7)^lambda1)*(x(9)^(1-lambda1));
-x(10)+x(8)/x(7);
-x(11)+x(8)/x(9);
%  //u*p=v*q;
-x(12)+x(7)/x(9);
-x(16)+((x(2)^alfa1)*(x(4)^(1-alfa1)))*x(18);
-x(16)+x(5)+x(6)+gamma1*x(7);
-(gamma1/x(10))+beta1*(1-delta)*(x(1)/x(1))*(alfa1*(x(16)/x(2))-x(13)+gamma1/x(10));
-((x(15)-x(1)*d1)/x(11))+beta1*(1-delta1)*(x(1)*x(13)-x(15)+((x(15)-x(1)*d1)/x(11)));
-x(13)+mi1*(alfa1*(x(16)/x(2)))+(1-mi1)*(x(15)/x(1));
-x(17)+(1-alfa1)*(x(16)/x(4));
-x(1)+(1/x(5));
-x(14)+1-x(9)-x(2);
-x(15)+H1*(x(14)^(-omega1));
-x(18)+rho_z1*x(18)+x(19)];
   