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

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