function [Hp]=house(H) % global beta expenditures; Hp=(1-beta)*expenditures./H; return