changes to rk45 to showcase two differnet methods

Added a new method to RK4 for calculating the system for apply_derivT the way John did. Put in a switch for going back and forth between the two methods.