function  u = qr2u(x,p);
%QR2U  is the utility function called by qr2lq.m
%       u = qr2u(x,p)  takes inputs x=[k,1,ome,k'-k] and p=[alpha,tau]
%       and computes:
%
%                 u = (exp(ome)*k^alpha-k'+k)^(1-tau)/(1-tau)
%

%       Ellen R. McGrattan, 3-23-89

alpha = p(1);
tau   = p(2);

u = (exp(x(3))*x(1)^alpha-x(4))^(1-tau)/(1-tau);

