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

global alpha beta beta delta psi rho tau theta 
z(1) = y(3)*theta*y(4)^(1+psi) -((1-alpha)*y(6));
z(2) = y(5) -(beta*(((exp(y(2))*y(3))/(exp(y(2))*y(3)))*(exp(y(2))*alpha* ...
y(6)+(1-delta)*y(5))));
z(3) = y(6) -(exp(y(1))*(y(5)^alpha)*(y(4)^(1-alpha)));
z(4) = y(5) -(exp(y(2))*(y(6)-y(3))+(1-delta)*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));
