Finish implementing the abstract Integrator class and a subclass for implementing RK4 approximations. The Integrator class will need a DiffEQ definition that will be fed into different algorithms in the subclasses.
Perhaps the issue with the integration blowing up is that there are discontinuities created by using a local fit rather than using a global fit for the function. Perhaps a global fit will fix this issue.
Finish implementing the abstract Integrator class and a subclass for implementing RK4 approximations. The Integrator class will need a DiffEQ definition that will be fed into different algorithms in the subclasses.