How do you make a numerical model that shows the motion of the Earth-Sun-moon system? Here, I use vpython to create the calculation.