text: 
  import  matplotlib.pyplot as plt
  from matplotlib import animation, rc
  from matplotlib.animation import FuncAnimation
  import numpy as np
  from IPython.display import HTML, Image
  
  rc('animation', html='html5')
  
  
  data = np.load('../../data/Orbits.npy')
  Positions = data[:,:,(0,2)]
  
  figure, ax = plt.subplots()
  
  # Setting limits for x and y axis
  r = 2
  ax.set_xlim(-r, r)
  ax.set_ylim(-r, r)
  
  ax2 = ax.twinx()
  ax2.set_xlim(-r, r)
  ax2.set_ylim(-r, r)
  
  ax3 = ax.twinx()
  ax3.set_xlim(-r, r)
  ax3.set_ylim(-r, r)
  
  ax4 = ax.twinx()
  ax4.set_xlim(-r, r)
  ax4.set_ylim(-r, r)
  
  ax5 = ax.twinx()
  ax5.set_xlim(-r, r)
  ax5.set_ylim(-r, r)
  
  # Plotting graph for bodies
  S, = ax.plot(0, 0, 'y.-', label='Sun')
  M, = ax2.plot(0, 0, 'm.-', label='Mercury')
  V, = ax3.plot(0, 0, 'c.-', label='Venus')
  E, = ax4.plot(0, 0, 'b.-', label='Earth') 
  A, = ax5.plot(0, 0, 'r.-', label='Mars')
  
  
  ax.legend([S, M, V, E, A], [S.get_label(), M.get_label(), V.get_label(), E.get_label(), A.get_label()], loc=0)
    
  def animation_function(i):
  
      S.set_data(Positions[:,0,0][i], Positions[:,0,1][i])
      M.set_data(Positions[:,1,0][i], Positions[:,1,1][i])
      V.set_data(Positions[:,2,0][i], Positions[:,2,1][i])
      E.set_data(Positions[:,3,0][i], Positions[:,3,1][i])
      A.set_data(Positions[:,4,0][i], Positions[:,4,1][i])
      return 
  
  
  
  anim = FuncAnimation(figure,func = animation_function,interval = 10)
  plt.show()