问题如下：我在磁场中有一个粒子，该粒子在所有轴上都有分量，并且与所有坐标有关
B =（B_x（x，y，z），B_y（x，y，z），B_z（x，y，z））
我求解了运动方程，发现了粒子在每个时间步的位置（x，y，z）和速度（v_x，v_y，v_z）。
我想做的是计算第二个绝热不变量J，定义为： Second adiabatic invariant，其中a和b是磁场中粒子的转折点，v_parallel是粒子速度与磁场速度的平行部分。另外，我试图找到第三个绝热不变量phi，即磁通量。我想这样做是为了查看它们是否（在时间上）被保存。
我已经获得了速度的平行和垂直分量，但是我不确定如何从数值上获得第二和第三不变性。 我很乐意为此提供一些帮助。
谢谢。
