code: def accel_ES_fl(q: np.ndarray):
      
      # Number of celestial bodies
      num_bodies: int = 2
      # Number of dimensions in arrays; 3 spatial dimensions times the number of bodies
      dims = 3 * num_bodies
  
      # The force given the positions q of the bodies
      f = force_ES(q)
  
      # The accelerations from this force
      a = np.zeros(dims)
      a[slices[0]] = f[slices[0]] / mass[0]
      a[slices[1]] = f[slices[1]] / mass[1]
      
      return a
docs: Accelaration in the earth-sun system using Fluxion potential energy