
function z=example2_fff(y)
z=zeros(6,1);
global ex_ it_ recur_

global alpha beta beta delta psi rho tau theta 
z(1) = exp(y(3))*theta*exp(y(4))^(1+psi) -((1-alpha)*exp(y(6)));
z(2) = exp(y(5)) -(beta*(((exp(y(2))*exp(y(3)))/(exp(y(2))*exp(y(3))))*(exp ...
(y(2))*alpha*exp(y(6))+(1-delta)*exp(y(5)))));
z(3) = exp(y(6)) -(exp(y(1))*(exp(y(5))^alpha)*(exp(y(4))^(1-alpha)));
z(4) = exp(y(5)) -(exp(y(2))*(exp(y(6))-exp(y(3)))+(1-delta)*exp(y(5)));
z(5) = y(1) -(rho*y(1)+tau*y(2)+ex_(it_-1,1));
z(6) = y(2) -(tau*y(1)+rho*y(2)+ex_(it_-1,2));
