Strange changes to internal imports that seem to be needed.
Extend BDF formula generator to work for equating to a derivative at any point

See Hairer et. al. for what I mean...
Attempt to make adaptive midpoint with eBDF and averages of derivatives
Add assert_sym_eq for comparing symbolic expressions