bmt function

data from Section 1.3