原子尺度简谐振动系统Python欧拉法代码异常问题求助

问题描述 #

我用欧拉法编写代码求解简谐振动系统，将二阶常微分方程拆分为两个一阶方程：
x' = v
v' = -k/m * x
目标是求解x和v随时间变化的函数并绘制x-t曲线。

本次建模针对钠原子，参数如下：

• 质量 m = 3.82×10^-26 kg
• 劲度系数 k = 12.2 N/m
• 初始位置 x₀ = 1.0×10^-10 m
• 初始速度 v₀ = 0 m/s

此前我有一个求解螺旋弹簧简谐振动系统的程序运行完全正常，但该系统参数数值大得多。本次仅修改上述参数后保留原程序其他部分，结果图形未呈现预期的正弦曲线，且持续出现溢出警告。

提问 #

还需要修改哪些内容才能让代码适配当前参数？

（注：原问题附带Python代码截图及结果图，因要求移除外链故未展示）

内容的提问来源于stack exchange，提问作者anon9876