function F=Fprova(xd, beta, rho, s, eta, omegab, omegaf, omegah, q)

F=[-xd(1)+xd(2)+beta*(rho*xd(1)+(1-rho)*xd(3));

-xd(2)+eta*(1+rho*beta*(xd(5)/q+xd(7)))+(1-eta)*(xd(4)+beta*(s-rho)*(xd(1)-xd(3)));

- xd(3)+xd(4)+beta*(s*xd(1)+(1-s)*xd(3));  

- xd(4)+omegab*xd(2);

- xd(5)+((1-xd(2))/(1-rho*beta)+xd(7))/q;

- xd(6)+omegaf*xd(2);

- xd(7)+omegah*xd(2)];









%- xd(1)+(1-((beta*(1-rho*(1-Gzr))*(beta*s*(1-Gzr)))/((1-rho*beta*(1-Gzr))*(1-beta*(1-s*(1-Gzr))))))^(-1)*xd(5)/(1-rho*beta*(1-Gzr))+beta*(1-rho*(1-Gzr))*xd(2)/((1-rho*beta*(1-Gzr))*(1-beta*(1-s*(1-Gzr))));

%-
%xd(3)+xd(2)/(1-beta*(1-s*(1-Gzr)))+beta*s*(1-Gzr)*xd(1)/(1-beta*(1-s*(1-Gzr)));