Experiment: Demonstrate Kozai Mechanism
Initial conditions: The kozai.xml file provided by Dan.
Parameters: RK4 accuracy, Timestep 30 minutes (except at the very end, where it was reduced to more accurately check for collision).
Observations: Orbit of moon oscillates around inclination 90, tending more towards positive growth rather than negative. During periods of greatest positive deviation from 90, eccentricity grows markedly, while it comes to a standstill at periods of greatest negative deviation. With each such cycle, the effect becomes more pronounced, and the deviation in inclination greater. Eccentricity never reduces; it only grows exponentially, and as it does so, deviation from inclination 90 becomes increasingly pronounced.
Collision took place near the end of the 5th year, with inclination exceeding 110 degrees, and eccentricity somewhere between 0.99 and 1. Timestep was reduced to 5 minutes for the last few cycles.
Regrettably, I have no visual data for that final year which exhibited the greatest changes - I should not have stuck so doggedly to the policy of taking one picture at the end of each year.
Below is an orbital plot covering the time from the beginning of the simulation to the end of the 4th. year. Each year is color coded to make it easier to see the age of the data.
Conclusion: The Kozai Mechanism can be demonstrated with US2