function F = nl_tot(x)
alpha=0.3;
alpha_star=0.15;
theta=0.4;
theta_star=0.4;
eta=1.5;
q=0.61;
q_star=1/q;
F=[(1-alpha)*((1-theta)*x(1))^(1-eta)+alpha*((1-theta)*q*x(2))^(1-eta)-(1+theta*x(1))^(1-eta);
(1-alpha_star)*((1-theta_star)*x(2))^(1-eta)+alpha_star*((1-theta_star)*q_star*x(1))^(1-eta)-(1+theta_star*x(2))^(1-eta)];
end