Found bug in imf.py when power = -1. Added fix and unit tests.
Lots of miscellaneous plots for different papers