#ploting 3 stars phi phi1 = now_arrayphi_star1[i] phi2 = now_arrayphi_star2[i] phi3 = now_arrayphi_star3[i] phiomega = arrayphi_omega[i] plt.figure() plt.plot(-ts,phi1,label='Star1') plt.plot(-ts,phi2,label='Star2') plt.plot(-ts,phi3,label='Star3') plt.plot(-ts,phiomega,label='Omega') plt.xlabel('t (Gyr)') plt.ylabel('phi (deg)') plt.legend() plt.show() star1.animate(d1=['x','R'],d2=['y','z'],width=800) o.plot3d() star1.plot3d() star2.plot3d() star3.plot3d() # In[221]: #animating Omega Cen o.animate(d1=['x','ra'],d2=['y','dec'],width=600) # In[222]: