我正在寻找动画我的图形(下面)和我不确定从何处或如何开始，因为我没有动画经验。我不知道它是如何工作的，也不知道代码的结构应该是什么，所以如果有人能提供一个伪代码或一个算法，我将非常感激。我也提供了我用来绘制下图的代码。在enter code here

from scipy.integrate import odeint

import scipy as sci

import numpy as np

import matplotlib.pyplot as plt

from mpl_toolkits.mplot3d import Axes3D

import matplotlib.animation as ani

# Universal Gravitational Const.

G = 6.674e-11

# Defining Mass

m1 = 2

m2 = 3.5

m3 = 2.3

# Init positions in graph (array)

pos1 = [-5,0,1]

pos2 = [5,0,10]

pos3 = [0,1,3]

p01 = np.array(pos1)

p02 = np.array(pos2)

p03 = np.array(pos3)

# Init velocities (array)

vi1 = [1,0.01,0]

vi2 = [-5,0,1]

vi3 = [0,-1,0]

v01 = np.array(vi1)

v02 = np.array(vi2)

v03 = np.array(vi3)

#Function

def derivs_func(y,t,G,m1,m2,m3):

d1 = y[:3]

d2 = y[3:6]

d3 = y[6:9]

v1 = y[9:12]

v2 = y[12:15]

v3 = y[15:18]

dist12 = np.sqrt((pos2[0]-pos1[0])**2 + (pos2[1]-pos1[1])**2 + (pos2[2]-pos1[2])**2)

dist13 = np.sqrt((pos3[0]-pos1[0])**2 + (pos3[1]-pos1[1])**2 + (pos3[2]-pos1[2])**2)

dist23 = np.sqrt((pos3[0]-pos2[0])**2 + (pos3[1]-pos2[1])**2 + (pos3[2]-pos2[2])**2)

dv1dt = m2 * (d2-d1)/dist12**3 + m3 * (d3-d1)/dist13**3

dv2dt = m1 * (d1-d2)/dist12**3 + m3 * (d3-d2)/dist23**3

dv3dt = m1 * (d1-d3)/dist13**3 + m2 * (d2-d3)/dist23**3

dd1dt = v1

dd2dt = v2

dd3dt = v3

derivs = np.array([dd1dt,dd2dt,dd3dt,dv1dt,dv2dt,dv3dt])

derivs3 = derivs.flatten()

return derivs3

yo = np.array([p01, p02, p03, v01, v02, v03])

y0 = yo.flatten()

time = np.linspace(0,200,500)

sol = odeint(derivs_func, y0, time, args = (G,m1,m2,m3))

x1 = sol[:,:3]

x2 = sol[:,3:6]

x3 = sol[:,6:9]

fig = plt.figure(figsize = (15,15))

ax = fig.add_subplot(111,projection = '3d')

ax.plot(x1[:,0],x1[:,1],x1[:,2],color = 'b')

ax.plot(x2[:,0],x2[:,1],x2[:,2],color = 'm')

ax.plot(x3[:,0],x3[:,1],x3[:,2],color = 'g')

ax.scatter(x1[-1,0],x1[-1,1],x1[-1,2],color = 'b', marker = 'o', s=30, label = 'Mass 1')

ax.scatter(x2[-1,0],x2[-1,1],x2[-1,2],color = 'm', marker = 'o',s=90, label = 'Mass 2')

ax.scatter(x3[-1,0],x3[-1,1],x3[-1,2],color = 'g', marker = 'o',s=60, label = 'Mass 3')

ax.legend()