DAY 3
Lecture Notes

Part 1: Day 3 Intro Slides
Reminder of what we have done/some intro to what we are doing today.

Part 2: Day 3 Intro Slides
Onto better solvers!

Part 3: Day 3 Intro Slides
Multibody systems, Hurray!
Code

Euler's solution with conservation of Energy and Momentum
Some code that does the Euler's method numerical solution to the 2-body problem and looks at the conservation of energy and momentum

Some code that compares Hermite and Euler Integration Schemes with Analytical
Python code comparing different methods of integration for the 2-Body problem.
Choose Your Own Adventure #1: START HERE!

CODE: Make a fake planetary system
Python code for solving for a multiplanet system.
Questions: (1) How can you tell if the system is stable? (2) What effects stability? (3) How hard is it to get a stable system? (4) What determines how long the computation takes?
Choose Your Own Adventure #2

Make a planetary system based on Kepler data
Read in some Kepler data from a confirmed planet system and plot its orbits.
Questions: (1) How can you tell if the system is stable? (2) What effects stability? (3) How long does each system take to compute orbits? Why?
Resources

Indepth Discussion of Several Integrators
Goes over some pros/cons of different sorts of integrators. Hermite integrators are discussed in section 2.3.2.

Super Planet Crash!
Can you make a stable system?

Bonus: Orbits Game
Can you please your alien overlords by building them awesome planetary systems?