Numerical differentiation for scores, Fisher informations

This feature has been working well in testing for quite a while now, such that I think it makes sense to merge it in to the main branch.