一位理论家认为,太阳系外行星的引力吸引使恒星在椭圆轨道上围绕系统质量中心移动。天文学家测量(通过多普勒位移)相对速度vk = v·n,其中v是恒星的真正速度,n是我们视线范围内的单位向量。我们碰巧看到了这个系统的“边缘”,所以视线矢量在运动平面上。理论家提出了一种供你调查的颂歌体系,以下的图中(r两点代表二阶导数)
这里 r = r (t)和φ=φ(t)。该模型有三个维度参数:p、e,对应椭圆的大小和偏心;和φ0,它们对应于椭圆的方向。通过引入一个新的变量s =\dot{r}(r这里的意思是一阶导数),在初始条件s(0) = 0的情况下,将ode写成3个自控的一阶方程。编写代码数值求解ODE系统,求出{r(t), s(t), φ(t)}在t∈[0,200]上给定参数{p, e, φ0},使用scipy. integration .solve ivp()。