decay_* functions should return a function that inherits from "decay_function"

[dhersz]

This way we can check if the functions passed gravity model are in fact decay functions and not, let's say, print().