Added 4th order symplectic integrator

[mckeon-ben]

I adapted this code as part of my master's thesis. Specifically, I added Forest and Ruth's method, a fourth-order symplectic integrator, resulting in substantially reduced render times vs. RK4. The results are explained in further detail in this paper.