Beginnings of code for interpolation model (MacKay and Peto equation 1)