Is your feature request related to a problem? Please describe.
There exists a finite difference solver in OSTk astro, the boost finite difference solver can be used to reduce code duplication and leverage it's performance
Describe the solution you'd like
https://www.boost.org/doc/libs/1_87_0/boost/math/differentiation/finite_difference.hpp