rk4_tested

A few minor changes in the rk4 and func functions were needed.
Additionally, I've included a 1D harmonic oscillator function to test
the rk4. It seems the rk4 is, finally, correct.