updated to-do list to reflect that I added the delta parameter, edited some comments in work.R
made changes so that the updating of some variables and evaluating of P(theta_a + delta < theta_b | data) is in its own seperate function