利用RK4积分求解Henon Heiles问题并生成Poincare截面
主要函数rk_demo.py绘制了Henon-Heiles轨道问题在两种不同初始能量下的解。该程序调用了rk4_poincare.py库,其中包括基本的四阶Runge-Kutta积分方案,能够实现数值解中的零交叉检测,并生成x=0平面的Poincare映射。以下是该实现的样本输出图像,展示了不同能量下轨道的分布与映射特征。
下载地址
用户评论