%solution
x0=[1 1];
xis=fsolve(@nl_tot,x0);
XI_h=xis(1,1);
XI_f_star=xis(1,2);