function uss = utilfun(uvar)

% ****************************************************
% ****************************************************

%parameter values
beta = .96;
alpha = .29;
gamma = 2;
theta = .6;
omega= 1.42;


k = uvar(1);
x = uvar(2);
l = uvar(3);
h = uvar(4);
kprime = uvar(5);

% Objective function:
uss = (( ((k*h)^alpha)*(l^(1-alpha)) + (1 - (h^omega)/omega) *exp(-x) *k -exp(-x)* kprime - (l^(1+theta))/(1+theta))^(1-gamma)-1)/(1-gamma);