我试图复制柯蒂斯在轨道力学上的插曲，但无法完全理解。 但是，我通过从np.arctan切换到np.arctan2取得了进展。
我可能没有正确实现arctan2吗？ 导入实验室
导入编号为NP
e=NP.arange (0.0，1.0，0.15 ) (.reshape (-1，1 ) ) ) ) ) ) ) ) ) ) ) )。
nu=NP.linspace (0.001，2 * NP.pi-0.001，50000 ) ) ) ) ) ) ) nu=NP.linspace (0.001，2 * NP.pi-0.001，50000 ) )
m2 evals=(2* NP.arctan2(1，1/) ) (1-e ) ) 0.5*NP.tan ) nu/2 ) -
e*(1-e**2) *0.5*NP.sin(nu )1e*NP.cos ) nu ) )
fig2=pylab.figure (
ax2=fig2.add_subplot(111 ) )
for Me2，_einzip(m2evals，e.ravel () ) :
ax2.plot(nu.ravel )、Me2、label=str ) _e ) )
pylab.legend (
pylab.xlim ()，7.75 ) )
py lab.ylim (((0，2 * NP.pi ) ) ) ) ) ) ) ) ) ) ) ) ) ) ) ) )。
pylab.show () )。
下图中出现了中断。 此函数是平滑的，应该在y范围(0，2pi )内连接0和2pi，而不是使0和2pi接触。
教科书的故事和方程式：
cqdyf的要求，我说：
问题可能出在以“原理值”为输出的arctan函数上。
因此，如果x是第二象限或第三象限的角，Arctan(tan(x ) )不会生成x。 如果将Arctan(tan ) x )从x=0绘制到x=Pi，则可以看到在x=Pi/2时存在不连续的跳跃。
在你的情况下，我想写arctan (1，1/arg )而不是arctan(arg )。 其中arg是Arctan函数的参数。 因此，如果arg为负值，arctan2将角度转向第二象限而不是第四象限。 ”