ynogkm：一种计算 Kerr-Newman 时空类时测地线的新公开代码

在本文中，我们提出了一种名为$ynogkm$的新公开代码，用于快速计算 Kerr-Newman (K-N) 时空中的类时测地线，这是计算 Kerr 时空中零测地线的$ynogk$代码的直接扩展。遵循$ynogk$中使用的策略，我们也利用 Weierstrass 和 Jacobi 椭圆函数及积分解析及半解析地求解了运动方程，其中 Boyer-Lindquist (B-L) 坐标$r$、$\theta$、$\phi$、$t$和固有时$\sigma$被表示为独立变量$p$（Mino 时间）的函数。所有椭圆积分均采用 Carlson 椭圆积分方法计算，保证了代码的快速运行。最后，将该代码应用于几个示例问题。