burn function

data from Section 1.6