视频来自B站：【【三体模拟】脆弱而优美的稳定三体系统】
视频使用Blender Python建模和渲染
使用RKF45（Runge-Kutta-Fehlberg）方法求运动方程的数值解，计算三体运动轨迹
* 该方法估计每步误差大小，并据此自适应调整积分步长
** 所有计算均取万有引力常数G=1.0
参考资料（三体运动初始条件）
https://observablehq.com/@rreusser/periodic-planar-three-body-orbits
https://numericaltank.sjtu.edu.cn/three-body/three-body.htm
http://three-body.ipb.ac.rs/
https://github.com/mws262/MAE5730_examples/tree/master/3BodySolutions
背景音乐
Tempus Fugit – Inner Music