More robust rejection calculation

currently a ZeroDivisionError is raised in the rejection calculation here.

In !370 (merged) a new save_divide function is introduced which would allow to return infinity for these cases. However, doing this will cause the plotting to fail, which needs to be adapted accordingly then.