copied in actual (modified from previous emails) code

code for integrating the system in time, that is
commented out exploratory prints
don't display the plot of all times, only dat times