I added some python code to check that the derivative matches (more or
less) the slope of the values on a unary function. Currently this is
just applied to harmonic, which is not extremely useful (but I do use
it on other not-yet submitted unary functions).