function F = myfun_GK_macro(c,sigmab,ksi,betta,theta,omega,chim,phi_1)
F(1)=c(2) - (1-sigmab/betta)*(1+omega*chim)/(c(1)*(sigmab+ksi)-sigmab/betta + omega*(sigmab+ksi)*chim*(1/betta-c(1))); %Equation 1
F(2)=c(2) - ((1-sigmab)+sigmab*c(2)*theta+phi_1*omega*(1/betta-c(1)))/(theta-((1-sigmab)+sigmab*c(2)*theta)*(betta*(c(1)-1)));%Equation 2
end
