#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]:
