Code running into problem

[jpfeifer]

1. How do the mode_check-plots look like?
2. For the unbounded density, see Unbounded density