%
% dpgq3F is a matlab function that defines the eulers
% at the 'steady-state'  
%

function [q] = dpgq3F(x,thetabar,alpha,beta,psi,sigma,delta)


% ----- define intra-euler ----- %

ee1 = -psi*(1-x(3))^-sigma + (1-alpha)*(1/x(2))*(x(1)^alpha)*(thetabar^(-alpha))*x(3)^(-alpha)    ;

% ----- define inter-euler ----- %

ee2 = -1/x(2)+beta*(1/thetabar)*(1/x(2))*(alpha*x(1)^(alpha-1)*thetabar^(1-alpha)*x(3)^(1-alpha)+(1-delta));

% ----- define resorce constraint ------ %

ee3 = x(1)^(alpha)*thetabar^(-alpha)*x(3)^(1-alpha)-x(2)-x(1)+(1-delta)*x(1)/thetabar;

% ----- define the equations to be returned ----- %

q = [ee1,ee2,ee3];