我对数值模拟（使用我自己的代码，而不是现成的包）磁场中单个带电粒子的运动很感兴趣。该场在时间和空间上会发生变化，但与回旋加速器频率或拉莫尔半径相比可能不会很快。

有许多众所周知的方法可以对场中粒子的运动方程进行积分。但是由于磁场产生的力具有一些特殊的性质，例如不改变粒子的动能。是否有一些数值程序更适合这种特殊情况，例如。这将提供更准确的模拟并更好地保存绝热不变量（何时应该保存）？