function F=myfunCapitalServices3(x)
% 
cN_c1 = evalin('base', 'cN_c1');
delta_1 = evalin('base', 'delta_1');
beta_1 = evalin('base', 'beta_1');
alphaN_1 = evalin('base', 'alphaN_1');
alphaT_1 = evalin('base', 'alphaT_1');
is_1 = evalin('base', 'is_1');
theta_1 = evalin('base', 'theta_1');
epsilon_1 = evalin('base', 'epsilon_1');

% cN_c1 = 0.53;
% delta_1 = 0.1;
% beta_1 = 1/1.04;
% alphaT_1 = (1-0.67);
% alphaN_1 = alphaT_1;
% is_1 = 0.02;
% theta_1 = 2;
% epsilon_1 = 2;

yT_y1 = x(1);
y1_S1 = x(2);
yT1_ST1=x(3);
mi1_phi1=x(4) ;
omega_1=x(5) ;
pNyN1_SN1=x(6); 
qa1_ss=x(7);

F=[yT_y1 - (1-cN_c1*(1-delta_1/y1_S1)) ; 
   y1_S1 - ((yT_y1/yT1_ST1)^((epsilon_1-1)/epsilon_1) + ((1-yT_y1)/pNyN1_SN1)^((epsilon_1-1)/epsilon_1))^(-epsilon_1/(epsilon_1-1)) ;
   yT1_ST1^((epsilon_1-1)/epsilon_1) - ((1-beta_1*(1-delta_1))/((beta_1*alphaT_1*mi1_phi1/qa1_ss)*(y1_S1*yT_y1)^(1/epsilon_1))) ;
   mi1_phi1 - (omega_1*(yT_y1 - is_1)^(-1/theta_1))*((omega_1*(yT_y1-is_1)^((theta_1-1)/theta_1) + (1-omega_1)*is_1^((theta_1-1)/theta_1))^(1/(theta_1-1))) ;
   omega_1 - ((yT_y1/is_1 - 1)^(1/theta_1))/(1+(yT_y1/is_1 - 1)^(1/theta_1)) ;
   pNyN1_SN1^((epsilon_1-1)/epsilon_1) - ((1-beta_1*(1-delta_1))/(beta_1*alphaN_1*(y1_S1*(1-yT_y1))^(1/epsilon_1))) ;
   qa1_ss - ((omega_1*(yT_y1 - is_1)^((theta_1-1)/theta_1) + (1-omega_1)*is_1^((theta_1-1)/theta_1))^(theta_1/(theta_1-1)))/yT_y1 ;];
  