Reproducibility of results

The calc_p_r script is working, but for some reason it doesn't output the numbers I remember (I get other plausible numbers instead). Something is not right here. 

Either I misremember my old experients or I've broken my code at some point.