hodg function

data from Section 1.10