Re-add Jacobian adaptive midpoint, generalise interpolators

* Need Jacobian methods to compare with currect methods

* Need to be able to use alternative methods of estimating y_nph etc.