function [util_work,util_pens]=util(p,Hp)
%function [util_work,util_pens]=util(p,Hp)
global w beta;
util_work=log(w)-beta*log(p)-(1-beta)*log(Hp);
util_pens=-beta*log(p)-(1-beta)*log(Hp);
return