function z = fixa_micro(x1,delta,delta_e,a,beta,alpha,omega,PHI,mc)
%x1=n
x0(1)=x1(1);
x0(2)=x1(2);
x0(3)=x1(3);

f = @(g)p(g,delta,delta_e,a,beta,alpha,omega,PHI,mc); %The anonomous function
z = fsolve(f,x0,optimset('TolFun',1e-10,'TolX',1e-10,'Display','off'));
% [end_varde z] = fminsearch(f,x0);



function F = p(x,delta,delta_e,a,beta,alpha,omega,PHI,mc)
n=x(1);
h=x(2);
k=x(3);


% F= [k^alpha-(eta*a*(n^(gamma+alpha))*(((k^alpha)*(n^(1-alpha))-delta*k)^sigma_c)/((beta*theta*kappa-(1-eta))*(1-alpha)*(a-1)))];

%More eq: F = [2*x(1) - x(2) - exp(-x(1));
                   % -x(1) + 2*x(2) - exp(-x(2))];
                   
                   


                   
                   
                   
  F = [(((delta_e*a*beta)/((1-beta*(1-delta_e))*(1-alpha)*(1-h)))^(omega/(omega-1)))*((h/(1+(((1-h)/delta_e)^a)*PHI))^(1/(omega-1)))-n;
 (1+(((1-h)/delta_e)^a)*PHI)/(((1-beta*(1-delta))/(beta*alpha*mc))^omega)-k;
 (((1-h)/delta_e)^a)*((alpha*k^((omega-1)/omega)+(1-alpha)*(h*n)^((omega-1)/omega))^(omega/(omega-1))-PHI)-1];
                   
                   




end


end